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ABSTRACT 



Abstract: This paper presents the first cosmological results based on Planck measurements of the cosmic microwave background (CMB) temper- 
ature and lensing-potential power spectra. We find that the Planck spectra at high multipoles (t > 40) are extremely well described by the standard 
£f*l spatially-flat six-parameter ACDM cosmology with a power-law spectrum of adiabatic scalar perturbations. Within the context of this cosmology, 
("I") the Planck data determine the cosmological parameters to high precision: the angular size of the sound horizon at recombination, the physical den- 
fV~) sities of baryons and cold dark matter, and the scalar spectral index are estimated to be ft, = (1.04147±0.00062)x 10~ 2 , Sl b h 2 = 0.02205 ±0.00028, 
H Si^h 2 = 0.1199 ± 0.0027, and n s = 0.9603 ± 0.0073, respectively (68% errors). For this cosmology, we find a low value of the Hubble constant, 
!* H = 67.3 ± 1.2 km s 1 Mpc~', and a high value of the matter density parameter, fi m = 0.3 15 ±0.017. These values are in tension with recent direct 
. £h measuremen ts of H and the magnitude-redshift relation for Type la supernovae, but are in excellent agreement with geometrical constraints from 
k> baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent level preci- 
sion using Planck CMB data alone. We use high-resolution CMB data together with Planck to provide greater control on extragalactic foreground 
components in an investigation of extensions to the six-parameter ACDM model. We present selected results from a large grid of cosmological 
models, using a range of additional astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured 
over the standard six-parameter ACDM cosmology. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor 
modes and to changes in the matter content of the Universe. We find a 95% upper limit of r .oo2 < 0.1 1 on the tensor-to-scalar ratio. There is no 
evidence for additional neutrino-like relativistic particles beyond the three families of neutrinos in the standard model. Using BAO and CMB data, 
we find /V cff = 3.30 ± 0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the sum of neutrino masses. 
Our results are in excellent agreement with big bang nucleosynthesis and the standard value of iV e fr = 3.046. We find no evidence for dynamical 
dark energy; using BAO and CMB data, the dark energy equation of state parameter is constrained to be w = -1.13*qJq. We also use the Planck 
data to set limits on a possible variation of the fine-structure constant, dark matter annihilation and primordial magnetic fields. Despite the success 
of the six-parameter ACDM model in describing the Planck data at high multipoles, we note that this cosmology does not provide a good fit to the 
temperature power spectrum at low multipoles. The unusual shape of the spectrum in the multipole range 20 < t < 40 was seen previously in the 
WMAP data and is a real feature of the primordial CMB anisotropies. The poor fit to the spectrum at low multipoles is not of decisive significance, 
but is an "anomaly" in an otherwise self-consistent analysis of the Planck temperature data. 

Key words. Cosmology: observations - Cosmology: theory - cosmic microwave background - cosmological parameters 
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1. Introduction 

The discovery of the cosmic microwave background (CMB) by 
Penzias & Wilson (1965) established the modem paradigm of 
the hot big bang cosmology. Almost immediately after this sem- 
inal discovery, searches began for anisotropics in the CMB - 
the primordial signatures of the fluctuations that grew to form 
the structure that we see today 1 . Convincing evidence for a 
dipole anisotropy was reported by Smoot et al. (1977), but de- 
spite many attempts, the detection of higher-order anisotropies 
proved elusive until the first results from the Cosmic Background 
Explorer (COBE; Smoot et al. 1992). The COBE results estab- 
lished the existence of a nearly scale-invariant spectrum of pri- 
mordial fluctuations on angular scales larger than 7°, consistent 
with the predictions of inflationary cosmology, and stimulated a 
new generation of precision measurements of the CMB of which 
this set of papers forms a part. 

CMB anisotropies are widely recognized as one of the 
most powerful probes of cosmology and early-Universe physics. 
Given a set of initial conditions and assumptions concerning 
the background cosmology, the angular power spectrum of the 
CMB anisotropies can be computed numerically to high preci- 
sion using linear perturbation theory (see Sect. 2). The combi- 
nation of precise experimental measurements and accurate the- 
oretical predictions can be used to set tight constraints on cos- 
mological parameters. The influential results from the Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite (Bennett et al. 
2003; Spergel et al. 2003), following on from earlier ground- 
based and sub-orbital experiments 2 , demonstrated the power of 
this approach, which has been followed by all subsequent CMB 
experiments. 

Planck 3 is the third-generation space mission, following 
COBE and WMAP, dedicated to measurements of the CMB 
anistropies. The primary aim of Planck (Planck Collaboration 
2005) is to measure the temperature and polarization 
anisotropies with micro-Kelvin sensitivity per resolution ele- 
ment over the entire sky. The wide frequency coverage of Planck 
(30-857 GHz) was chosen to provide accurate discrimination 
of Galactic emission from the primordial anisotropies and to 
enable a broad range of ancilliary science, such as detections 
of galaxy clusters, extragalactic point sources and the prop- 
erties of Galactic dust emission. This paper, one of a set as- 
sociated with the 2013 release of data from the Planck mis- 
sion (Planck Collaboration 1 2013), describes the first cosmolog- 
ical parameter results from the Planck temperature power spec- 
trum. 

The results from WMAP (see Bennett et al. 2012 
and Hinshaw et al. 2012 for the final nine-year WMAP re- 
sults) together with those from high-resolution ground-based 



'Corresponding author: G. Efstathiou <gpe@ast . cam. ac .uk> 
'For a good review of the early history of CMB studies see 
Peebles et al. (2009). 

2 It is worth highlighting here the pre- WMAP constraints on 
the geometry of the Universe by the BOOMERang (Balloon 
Observations of Millimetric Extragalactic Radiation and 
Geomagnetics; de Bernardis et al. 2000) and MAXIMA (Millimeter- 
wave Anisotropy Experiment Imaging Array; Balbi et al. 2000) 
experiments. 

^Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



CMB experiments (e.g., Reichardt et al. 2012; Story et al. 
2012; Sievers et al. 2013) are remarkably consistent with the 
predictions of a "standard" cosmological model. This model is 
based upon a spatially-flat, expanding Universe whose dynamics 
are governed by General Relativity and whose constituents are 
dominated by cold dark matter (CDM) and a cosmological 
constant (A) at late times. The primordial seeds of structure 
formation are Gaussian-distributed adiabatic fluctuations with 
an almost scale-invariant spectrum. This model (which will be 
referred to as the base ACDM model in this paper) is described 
by only six key parameters. Despite its simplicity, the base 
ACDM model has proven to be successful in describing a wide 
range of cosmological data in addition to the CMB, including 
the Type la supernovae magnitude-distance relation, baryon 
acoustic oscillation measurements, the large-scale clustering of 
galaxies and cosmic shear (as reviewed in Sect. 5). 

Nevertheless, there have been some suggestions of new 
physics beyond that assumed in the base ACDM model. 
Examples include various large-angle "anomalies" in the CMB 
(as reviewed by the WMAP team in Bennett et al. 2011) and 
hints of new physics, such as additional relativistic particles, that 
might steepen the high multipole "damping tail" of the CMB 
temperature power spectrum (Dunkley et al. 2011; Hou et al. 
2012). Furthermore, developments in early-Universe cosmology 
over the last 20 years or so have led to a rich phenomenology (see 
e.g., Baumann 2009 for a review). It is easy to construct models 
that preserve the main features of simple single-field inflationary 
models, but lead to distinctive observational signatures such as 
non-Gaussianity, isocurvature modes or topological defects. 

A major goal of the Planck experiment is to test the ACDM 
model to high precision and identify areas of tension. From 
previous CMB experiments and other cosmological probes, 
we know that any departures from the standard six-parameter 
ACDM cosmology are likely to be small and challenging to de- 
tect. Planck, with its combination of high sensitivity, wide fre- 
quency range and all-sky coverage, is uniquely well-suited to 
this challenge. 

The focus of this paper is to investigate cosmological con- 
straints from the temperature power spectrum measured by 
Planck. Figure 1 summarizes some important aspects of the 
Planck temperature power spectrum; we plot this as = 
1(1 + l)C(/2n (a notation we will use throughout this pa- 
per) versus multipole I. The temperature likelihood used in 
this paper is a hybrid: over the multipole range I = 2-A9, 
the likelihood is based on a component-separation algorithm 
applied to 91% of the sky (Planck Collaboration XII 2013; 
Planck Collaboration XV 2013). The likelihood at higher multi- 
poles is constructed from cross-spectra over the frequency range 
100-217 GHz, as discussed in Planck Collaboration XV (2013). 
It is important to recognize that unresolved foregrounds (and 
other factors such as beams and calibration uncertainties) need 
to be modelled to high precision to achieve the science goals 
of this paper. There is therefore no unique ''Planck primordial 
temperature spectrum". Figure 1 is based on a full likelihood so- 
lution for foreground and other "nuisance" parameters assuming 
a cosmological model. A change in the cosmology will lead to 
small changes in the Planck primordial CMB power spectrum 
because of differences in the foreground solution. Neverthess, 
Fig. 1 provides a good illustration of the precision achieved by 
Planck. The precision is so high that conventional power spec- 
trum plots (shown in the upper panel of Fig. 1) are usually un- 
informative. We therefore place high weight in this paper on 
plots of residuals with respect to the best-fit model (shown in the 
lower panel). Figure 1 also serves to illustrate the highly inter- 
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Fig. 1. Planck foreground-subtracted temperature power spectrum (with foreground and other "nuisance" parameters fixed to their 
best-fit values for the base ACDM model). The power spectrum at low multipoles (I = 2-49, plotted on a logarithmic multi- 
pole scale) is determined by the Commander algorithm applied to the Planck maps in the frequency range 30-353 GHz over 
91% of the sky. This is used to construct a low-multipole temperature likelihood using a Blackwell-Rao estimator, as described 
in Planck Collaboration XV (2013). The asymmetric error bars show 68% confidence limits and include the contribution from un- 
certainties in foreground subtraction. At multipoles 50 < I < 2500 (plotted on a linear multipole scale) we show the best-fit CMB 
spectrum computed from the CamSpec likelihood (see Planck Collaboration XV 2013) after removal of unresolved foreground com- 
ponents. The light grey points show the power spectrum multipole-by-multipole. The blue points show averages in bands of width 
A{ * 31 together with lcr errors computed from the diagonal components of the band-averaged covariance matrix (which includes 
contributions from beam and foreground uncertainties). The red line shows the temperature spectrum for the best-fit base ACDM 
cosmology. The lower panel shows the power spectrum residuals with respect to this theoretical model. The green lines show the 
±lcr errors on the individual power spectrum estimates at high multipoles computed from the CamSpec covariance matrix. Note the 
change in vertical scale in the lower panel at t — 50. 
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connected nature of this series of papers. The temperature like- 
lihood used in this paper utilizes data from both the Planck Low 
Frequency Instrument (LFI) and High Frequency Instrument 
(HFI). The data-processing chains for these two instruments 
and beam calibrations are described in Planck Collaboration II 
(2013), Planck Collaboration VI (2013), and associated papers. 
Component separation is described in Planck Collaboration XII 
(2013) and the temperature power spectrum and likelihood, as 
used in this paper, are described in Planck Collaboration XV 
(2013). Planck Collaboration XV (2013) also presents a de- 
tailed analysis of the robustness of the likelihood to various 
choices, such as frequency ranges and sky masks (and also com- 
pares the likelihood to results from an independent likelihood 
code based on different assumptions, see also Appendix B). 
Consistency of the Planck maps across frequencies is demon- 
strated in Planck Collaboration XI (2013), and the level of con- 
sistency with WMAP is also assessed. 

This paper is closely linked to other papers reporting 
cosmological results in this series. We make heavy use of 
the gravitational lensing power spectrum and likelihood es- 
timated from an analysis of the 4-point function of the 
Planck maps (Planck Collaboration XVII 2013). The present 
paper concentrates on simple parameterizations of the spec- 
trum of primordial fluctuations. Tests of specific models of 
inflation, isocurvature modes, broken scale-invariance etc. are 
discussed in Planck Collaboration XXII (2013). Here, we as- 
sume throughout that the initial fluctuations are Gaussian and 
statistically isotropic. Precision tests of non-Gaussianity, from 
Planck estimates of the 3- and 4-point functions of the temper- 
ature anisotropies, are presented in Planck Collaboration XXIV 
(2013). Tests of isotropy and additional tests of non-Gaussianity 
using Planck data are discussed in Planck Collaboration XXIII 
(2013). 

The outline of the paper is as follows. In Sect. 2 we define our 
notation and cosmological parameter choices. This section also 
summarizes aspects of the Markov Chain Monte Carlo (MCMC) 
sampler used in this paper and of the CMB Boltzmann code 
used to predict theoretical temperature power spectra. Section 3 
presents results on cosmological parameters using Planck data 
alone. For this data release we do not use Planck polarization 
data in the likelihood, and we therefore rely on WMAP polar- 
ization data at low multipoles to constrain the optical depth, r, 
from reionization. An interesting aspect of Sect. 3 is to assess 
whether CMB gravitational lensing measurements from Planck 
can be used to constrain the optical depth without the use of 
WMAP polarization measurements. 

Section 4 introduces additional CMB temperature data from 
high-resolution experiments. This section presents a detailed 
description of how we have modified the Planck model for 
unresolved foreground and "nuisance" parameters introduced 
in Planck Collaboration XV (2013) to enable the Planck spec- 
tra to be used together with those from other CMB experiments. 
Combining high-resolution CMB experiments with Planck mit- 
igates the effects of unresolved foregrounds which, as we will 
show, can affect cosmological parameters (particularly for ex- 
tensions to the base ACDM model) if the foreground parameters 
are allowed too much freedom. Section 4 ends with a detailed 
analysis of whether the base ACDM model provides an accept- 
able fit to the CMB temperature power spectra from Planck and 
other experiments. 

It is well known that certain cosmological parameter combi- 
nations are highly degenerate using CMB power spectrum mea- 
surements alone (Zaldarriaga et al. 1997; Efstathiou & Bond 
1999; Howlett et al. 2012). These degeneracies can be broken 



by combining with other cosmological data (though the Planck 
lensing analysis does help to break the principal "geometrical" 
degeneracy, as discussed in Sect. 5.1). Section 5 discusses ad- 
ditional "astrophysical" data that are used in combination with 
Planck. Since the Planck temperature data are so precise, we 
have been selective in the additional data sets that we have cho- 
sen to use. Section 5 discusses our rationale for making these 
choices. 

Having made a thorough investigation of the base ACDM 
model, Sect. 6 describes extended models, including models 
with non-power-law spectral indices, tensor modes, curvature, 
additional relativistic species, neutrino masses and dynamical 
dark energy. This section also discusses constraints on models 
with annihilating dark matter, primordial magnetic fields and a 
time-variable fine-structure constant. 

Finally, we present our conclusions in Sect. 7. Appendix A 
contrasts the Planck best-fit ACDM cosmology with that de- 
termined recently by combining data from the South Pole 
Telescope with WMAP (Story et al. 2012). Appendix B dis- 
cusses the dependence of our results for extended models 
on foreground modelling and likelihood choices, building on 
the discussion in Planck Collaboration XV (2013) for the base 
ACDM model. 



2. Model, parameters, and methodology 

2. 1 . Theoretical model 

We shall treat anisotropies in the CMB as small fluctua- 
tions about a Friedmann-Robertson- Walker metric whose evo- 
lution is described by General Relativity. We shall not con- 
sider modified gravity scenarios or "active" sources of fluc- 
tuations such as cosmic defects. The latter are discussed in 
Planck Collaboration XXV (2013). Under our assumptions, the 
evolution of the perturbations can be computed accurately us- 
ing a CMB Boltzmann code once the initial conditions, ion- 
ization history and constituents of the Universe are specified. 
We discuss each of these in this section, establishing our nota- 
tion. Our conventions are consistent with those most commonly 
adopted in the field and in particular with those used in the camb 4 
Boltzmann code (Lewis et al. 2000), which is the default code 
used in this paper. 

2.1 .1 . Matter and radiation content 

We adopt the usual convention of writing the Hubble constant 
at the present day as Hq = 100/ikms -1 Mpc -1 . For our base- 
line model, we assume that the cold dark matter is pressureless, 
stable and non-interacting, with a physical density u> c = Q. c h 2 . 
The baryons, with density o>b = Qb/1 2 , are assumed to consist 
almost entirely of hydrogen and helium; we parameterize the 
mass fraction in helium by Yp. The process of standard big bang 
nucleosynthesis (BBN) can be accurately modelled, and gives a 
predicted relation between Yp, the photon-baryon ratio, and the 
expansion rate (which depends on the number of relativistic de- 
grees of freedom). By default we use interpolated results from 
the PArthENoPE BBN code (Pisanti et al. 2008) to set Yp fol- 
lowing Hamann et al. (2011), which for the Planck best-fitting 
base model (assuming no additional relativistic components and 
negligible neutrino degeneracy) gives Yp = 0.2477. We shall 
compare our results with the predictions of BBN in Sect. 6.4. 
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Table 1. Cosmological parameters used in our analysis. For each, we give the symbol, prior range, value taken in the base ACDM 
cosmology (where appropriate), and summary definition (see text for details). The top block contains parameters with uniform 
priors that are varied in the MCMC chains. The ranges of these priors are listed in square brackets. The lower blocks define various 
derived parameters. 



Parameter Prior range Baseline Definition 



(Ob = Sl b h 2 [0.005, 0. 1] ... Baryon density today 

w c = Cl c h 2 [0.001,0.99] ... Cold dark matter density today 

100e M c [0.5,10.0] ... 100 x approximation to r./D A (CosmoMC) 

t [0.01,0.8] ... Thomson scattering optical depth due to reionization 

S1 K [-0.3,0.3] Curvature parameter today with fi tot = 1 - Q. K 

2 m v [0, 5] 0.06 The sum of neutrino masses in eV 

" 1 v ff s teriic [0, 3] Effective mass of sterile neutrino in eV 

wo [-3.0,-0.3] -1 Dark energy equation of state", w(a) = wq + (1 - a)w a 

w a [-2, 2] As above (perturbations modelled using PPF) 

N eS [0.05, 10.0] 3.046 Effective number of neutrino-like relativistic degrees of freedom (see text) 

K P [0.1,0.5] BBN Fraction of baryonic mass in helium 

A L [0, 10] 1 Amplitude of the lensing power relative to the physical value 

n s [0.9,1.1] ... Scalar spectrum power-law index (k = 0.05Mpc~') 

n t «t = -'"0.05/8 Inflation Tensor spectrum power-law index (kg = 0.05Mpc~') 

dn s /dlnk [-1,1] Running of the spectral index 

ln(10 10 A s ) [2.7,4.0] ... Log power of the primordial curvature perturbations (kg = 0.05 Mpc~') 

r .05 [0, 2] Ratio of tensor primordial power to curvature power at k = 0.05 Mpc~' 

fl A ... Dark energy density divided by the critical density today 

to ... Age of the Universe today (in Gyr) 

Q m ... Matter density (inc. massive neutrinos) today divided by the critical density 

<Tg ... RMS matter fluctuations today in linear theory 

z m ... Redshift at which Universe is half reionized 

Hq [20,100] ... Current expansion rate in km s~'Mpc~' 

ro.002 Ratio of tensor primordial power to curvature power at k = 0.002 Mpc~' 

10 9 A S ... 10 5 x dimensionless curvature power spectrum at ko = 0.05 Mpc~' 

co m = Cl m h 2 ... Total matter density today (inc. massive neutrinos) 

z, ... Redshift for which the optical depth equals unity (see text) 

r, = r s (z*) ... Comoving size of the sound horizon at z = z, 

1006, ... 100 x angular size of sound horizon at z = z, (?"»/Da) 

z drag ... Redshift at which baryon-drag optical depth equals unity (see text) 

'"drag = ''s(Zdiag) • • • • • • • Comoving size of the sound horizon at z = Zdrag 

ko ... Characteristic damping comoving wavenumber (Mpc~') 

1OO0D • • • 100 x angular extent of photon diffusion at last scattering (see text) 

z cq ... Redshift of matter-radiation equality (massless neutrinos) 

1OO0 CIJ ... 100 x angular size of the comoving horizon at matter-radiation equality 

?"drag/-Dv(0.57) .... ... BAO distance ratio at z = 0.57 (see Sect. 5.2) 



For dynamical dark energy models with constant equation of state, we denote the equation of state by w and adopt the same prior as for wo- 



The photon temperature today is well measured to be Tq = 
2.7255 ± 0.0006 K (Fixsen 2009); we adopt T = 2.7255 K as 
our fiducial value. We assume full thermal equilibrium prior to 
neutrino decoupling. The decoupling of the neutrinos is nearly, 
but not entirely, complete by the time of electron-positron anni- 
hilation. This leads to a slight heating of the neutrinos in addition 
to that expected for the photons and hence to a small departure 
from the thermal equilibrium prediction T y = (ll/4)^ 3 r v be- 
tween the photon temperature T y and the neutrino temperature 
T v . We account for the additional energy density in neutrinos by 
assuming that they have a thermal distribution with an effective 
energy density 



7 / 4 \ 4/3 

*-8 n p y> 



(i) 



with N e ff = 3.046 in the baseline model (Mangano et al. 2002, 
2005). This density is divided equally between three neutrino 
species while they remain relativistic. 

In our baseline model we assume a minimal-mass normal 
hierarchy for the neutrino masses, accurately approximated for 



current cosmological data as a single massive eigenstate with 
m v = 0.06 eV (Q. v h 2 * £ m v /93.04 eV * 0.0006; corrections 
and uncertainties at the meV level are well below the accuracy 
required here). This is consistent with global fits to recent os- 
cillation and other data (Forero et al. 2012), but is not the only 
possibility. We discuss more general neutrino mass constraints 
in Sect. 6.3. 



We shall also consider the possibility of extra radiation, 
beyond that included in the Standard Model. We model this 
as additional massless neutrinos contributing to the total A^ e ff 
determining the radiation density as in Eq. (1). We keep the 
mass model and heating consistent with the baseline model at 



(massive) 
eff 



A'eif = 3.046, so there is one massive neutrino with N[ 
3.046/3 ~ 1.015, and massless neutrinos with // m<issless ) = 
/V e ff - 1.015. In the case where /V e ff < 1.015 we use one mas- 
sive eigenstate with reduced temperature. 
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2.1.2. Ionization history 

To make accurate predictions for the CMB power spectra, the 
background ionization history has to be calculated to high ac- 
curacy. Although the main processes that lead to recombina- 
tion at z ~ 1090 are well understood, cosmological param- 
eters from Planck can be sensitive to sub-percent differences 
in the ionization fraction x e (Hu et al. 1995; Lewis et al. 2006; 
Rubino-Martin et al. 2009; Shaw & Chluba 2011). The process 
of recombination takes the Universe from a state of fully ion- 
ized hydrogen and helium in the early Universe, through to the 
completion of recombination with residual fraction x e ~ 10 4 . 
Sensitivity of the CMB power spectrum to x e enters through 
changes to the sound horizon at recombination, from changes 
in the timing of recombination, and to the detailed shape of the 
recombination transition, which affects the thickness of the last- 
scattering surface and hence the amount of small-scale diffusion 
(Silk) damping, polarization, and line-of-sight averaging of the 
perturbations. 

Since the pioneering work of Peebles (1968) and 
Zeldovich et al. (1969), which identified the main physical 
processes involved in recombination, there has been signif- 
icant progress in numerically modelling the many relevant 
atomic transitions and processes that can affect the details of 
the recombination process (Hu et al. 1995; Seageretal. 2000; 
Wongetal. 2008; Hirata & Switzer 2008; Switzer & Hirata 
2008; Rubino-Martin et al. 2009; Chluba & Thomas 2011; 
Ali-Haimoud & Hirata 2011). In recent years a consen- 
sus has emerged between the results of two multi-level 
atom codes HyRec 5 (Switzer & Hirata 2008; Hirata 2008; 
Ali-Haimoud & Hirata 2011), and CosmoRec 6 (Chluba et al. 
2010; Chluba & Thomas 2011), demonstrating agreement at a 
level better than that required for Planck (differences less that 
4 x 10~ 4 in the predicted temperature power spectra on small 
scales). 

These recombination codes are remarkably fast, given the 
complexity of the calculation. However, the recombination his- 
tory can be computed even more rapidly by using the sim- 
ple effective three-level atom model developed by Seager et al. 
(2000) and implemented in the recfast code 7 , with appropri- 
ately chosen small correction functions calibrated to the full 
numerical results (Wong et al. 2008; Rubino-Martin et al. 2009; 
Shaw & Chluba 201 1). We use recfast in our baseline param- 
eter analysis, with correction functions adjusted so that the pre- 
dicted power spectra C( agree with those from the latest ver- 
sions of HyRec (January 2012) and CosmoRec (v2) to better than 
0.05 % 8 . We have confirmed, using importance sampling, that 
cosmological parameter constraints using recfast are consis- 
tent with those using CosmoRec at the 0.05cr level. Since the re- 
sults of the Planck parameter analysis are crucially dependent on 
the accuracy of the recombination history, we have also checked, 
following Lewis et al. (2006), that there is no strong evidence for 
simple deviations from the assumed history. However, we note 
that any deviation from the assumed history could significantly 
shift parameters compared to the results presented here and we 
have not performed a detailed sensitivity analysis. 

The background recombination model should accurately 
capture the ionization history until the Universe is reionized 
at late times via ultra-violet photons from stars and/or active 

5 http : //www . sns . ias . edu/~yacine/hyrec/hyrec . html 
6 http: //www. chluba.de/CosmoRec/ 

7 http : //www . astro . ubc . ca/people/ scott/recf ast . html 
8 The updated recfast used here in the baseline model is publicly 
available as version 1.5.2 and is the default in camb as of October 2012. 



galactic nuclei. We approximate reionization as being relatively 
sharp, with the mid-point parameterized by a redshift z re (where 
x e = f /2) and width parameter Az le = 0.5. Hydrogen reion- 
ization and the first reionization of helium are assumed to oc- 
cur simultaneously, so that when reionization is complete x e = 
/ = 1 + fm ~ 1-08 (Lewis 2008), where fu e is the helium- 
to-hydrogen ratio by number. In this parameterization, the opti- 
cal depth is almost independent of Az re and the only impact of 
the specific functional form on cosmological parameters comes 
from very small changes to the shape of the polarization power 
spectrum on large angular scales. The second reionization of he- 
lium (i.e., He + — > He ++ ) produces very small changes to the 
power spectra (At ~ 0.001, where t is the optical depth to 
Thomson scattering) and does not need to be modelled in detail. 
We include the second reionization of helium at a fixed redshift 
of z = 3.5 (consistent with observations of Lyman-o' forest lines 
in quasar spectra, e.g., Becker et al. 201 1), which is sufficiently 
accurate for the parameter analyses described in this paper. 

2.1 .3. Initial conditions 

In our baseline model we assume purely adiabatic scalar per- 
turbations at very early times, with a (dimensionless) curvature 
power spectrum parameterized by 

P K (k)=A s {-\ , (2) 

with n s and dn^/dlnk taken to be constant. For most of this 
paper we shall assume no "running", i.e., a power-law spec- 
trum with dn s /dlnk = 0. The pivot scale, ko, is chosen to be 
ko = 0.05 Mpc~', roughly in the middle of the logarithmic range 
of scales probed by Planck. With this choice, n s is not strongly 
degenerate with the amplitude parameter A s . 

The amplitude of the small-scale linear CMB power spec- 
trum is proportional to e~ 2r A s . Because Planck measures this 
amplitude very accurately there is a tight linear constraint be- 
tween t and In A s (see Sect. 3.4). For this reason we usually use 
lnA s as a base parameter with a flat prior, which has a signifi- 
cantly more Gaussian posterior than A s . A linear parameter re- 
definition then also allows the degeneracy between t and A s to be 
explored efficiently. (The degeneracy between r and A s is broken 
by the relative amplitudes of large-scale temperature and polar- 
ization CMB anisotropics and by the non-linear effect of CMB 
lensing.) 

We shall also consider extended models with a significant 
amplitude of primordial gravitational waves (tensor modes). 
Throughout this paper, the (dimensionless) tensor mode spec- 
trum is parameterized as a power-law with 9 

(3) 

We define ro.05 = A t /A s , the primordial tensor-to-scalar ratio at 
k = ko. Our constraints are only weakly sensitive to the tensor 
spectral index, n t (which is assumed to be close to zero), and 
we adopt the theoretically motivated single-field inflation con- 
sistency relation n t = -ro.os/8, rather than varying n t indepen- 
dently. We put a flat prior on ro.05, t> ut a l so report the constraint 
at k = 0.002 Mpc~' (denoted ro.002), which is closer to the scale 

9 For a transverse-traceless spatial tensor H-,j, the tensor part of the 
metric is ds 2 = a 2 [dn 2 - (Sjj + 2Hjj)dx'dx j ], and P, is denned so that 
P t (k) = Sw t <2Hi ; 2H«>. 
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at which there is some sensitivity to tensor modes in the large- 
angle temperature power spectrum. Most previous CMB experi- 
ments have reported constraints on ro.rxn- For further discussion 
of the tensor-to-scalar ratio and its implications for inflationary 
models see Planck Collaboration XXII (2013). 

2.1 .4. Dark energy 

In our baseline model we assume that the dark energy is a cos- 
mological constant with current density parameter Q.^. When 
considering a dynamical dark energy component, we parame- 
terize the equation of state either as a constant w or as a function 
of the cosmological scale factor, a, with 

p 

w(a) = - = w + (1 - a)w a , (4) 
P 

and assume that the dark energy does not interact with other con- 
stituents other than through gravity. Since this model allows the 
equation of state to cross below -1, a single-fluid model can- 
not be used self-consistently. We therefore use the parameterized 
post-Friedmann (PPF) model of Fang et al. (2008a). For models 
with w > -1, the PPF model agrees with fluid models to signif- 
icantly better accuracy than required for the results reported in 
this paper. 

2.1.5. Power spectra 

Over the last decades there has been significant progress in 
improving the accuracy, speed and generality of the numer- 
ical calculation of the CMB power spectra given an ion- 
ization history and set of cosmological parameters (see e.g., 
Bond & Efstathiou 1987; Sugiyama 1995; Ma & Bertschinger 
1995; Seljak & Zaldarriaga 1996; Zaldarriaga et al. 1998; 
Lewis et al. 2000; Lesgourgues & Tram 201 1). Our baseline nu- 
merical Boltzmann code is camb 10 (March 2013; Lewis et al. 
2000), a parallelized line-of-sight code developed from cmbf ast 
(Seljak & Zaldarriaga 1996) and Cosmics (Bertschinger 1995; 
Ma & Bertschinger 1995), which calculates the lensed CMB 
temperature and polarization power spectra. The code has been 
publicly available for over a decade and has been very well tested 
(and improved) by the community. Numerical stability and ac- 
curacy of the calculation at the sensitivity of Planck has been 
explored in detail (Hamann et al. 2009; Lesgourgues 2011b; 
Howlett et al. 2012), demonstrating that the raw numerical pre- 
cision is sufficient for numerical errors on parameter constraints 
from Planck to be less than 10% of the statistical error around 
the assumed cosmological model. (For the high multipole CMB 
data at ( > 2000 introduced in Sect. 4, the default camb settings 
are adequate because the power spectra of these experiments are 
dominated by unresolved foregrounds and have large errors at 
high multipoles.) To test the potential impact of camb errors, we 
importance-sample a subset of samples from the posterior pa- 
rameter space using higher accuracy settings. This confirms that 
differences purely due to numerical error in the theory prediction 
are less than 10% of the statistical error for all parameters, both 
with and without inclusion of CMB data at high multipoles. We 
also performed additional tests of the robustness and accuracy 
of our results by reproducing a fraction of them with the inde- 
pendent Boltzmann code class (Lesgourgues 2011a; Bias et al. 
2011). 

In the parameter analysis, information from CMB lensing 
enters in two ways. Firstly, all the CMB power spectra are 

10 http : //camb . info 



modelled using the lensed spectra, which includes the approx- 
imately 5% smoothing effect on the acoustic peaks due to lens- 
ing. Secondly, for some results we include the Planck lensing 
likelihood, which encapsulates the lensing information in the 
(mostly squeezed-shape) CMB trispectrum via a lensing po- 
tential power spectrum (Planck Collaboration XVII 2013). The 
theoretical predictions for the lensing potential power spectrum 
are calculated by camb, optionally with corrections for the non- 
linear matter power spectrum, along with the (non-linear) lensed 
CMB power spectra. For the Planck temperature power spec- 
trum, corrections to the lensing effect due to non-linear struc- 
ture growth can be neglected, however the impact on the lens- 
ing potential reconstruction is important. We use the halofit 
model (Smith et al. 2003) as updated by Takahashi et al. (2012) 
to model the impact of non-linear growth on the theoretical pre- 
diction for the lensing potential power. 

2.2. Parameter choices 

2.2.1. Base parameters 

The first section of Table 1 lists our base parameters that have 
flat priors when they are varied, along with their default values 
in the baseline model. When parameters are varied, unless oth- 
erwise stated, prior ranges are chosen to be much larger than the 
posterior, and hence do not affect the results of parameter esti- 
mation. In addition to these priors, we impose a "hard" prior on 
the Hubble constant of [20, 100] kms" 1 Mpc" 1 . 

2.2.2. Derived parameters 

Matter-radiation equality z eq is defined as the redshift at which 
p 7 + p v — p c + pb (where p v approximates massive neutrinos as 
massless). 

The redshift of last-scattering, z», is defined so that the op- 
tical depth to Thomson scattering from z — (conformal time 
T] = T]o) to z = Zm is unity, assuming no reionization. The optical 
depth is given by 

r(77)= ffdrf, (5) 

•J '10 

where f = -an e cr T (and n e is the density of free electrons and 
<x T is the Thomson cross section). We define the angular scale of 
the sound horizon at last-scattering, 6* - rs(z*)/D\(z<,), where r s 
is the sound horizon 

n iz> drf 

r s (z)= ; (6) 

Jo V3(l +R) 

withfl = 3p b /(4p r ). 

Baryon velocities decouple from the photon dipole when 
Compton drag balances the gravitational force, which happens 
at Td ~ 1, where (Hu & Sugiyama 1996) 

T d (i/)s ffdrf/R. (7) 
J no 

Here, again, t is from recombination only, without reioniza- 
tion contributions. We define a drag redshift Zdmg, so that 
T d(?7(zdiag)) = 1 ■ The sound horizon at the drag epoch is an im- 
portant scale that is often used in studies of baryon acoustic os- 
cillations; we denote this as r,j ra g = r s (zd ra g)- We compute Zdrag 
and rjrag numerically from camb (see Sect. 5.2 for details of ap- 
plication to BAO data). 



7 



Planck Collaboration: Cosmological parameters 



The characteristic wavenumber for damping, kp, is given by 
1 



. r> , 1 /? 2 + i6(i+/?)/i5 

^^J, rf "f (1 + *) 2 



(8) 



We define the angular damping scale, 6>d = tt/ (k D Dpj, where Da 
is the comoving angular diameter distance to z». 

For our purposes, the normalization of the power spectrum 
is most conveniently given by A s . However, the alternative mea- 
sure erg is often used in the literature, particularly in studies of 
large-scale structure. By definition, erg is the rms fluctuation in 
total matter (baryons + CDM + massive neutrinos) in 8 hr x Mpc 
spheres at z = 0, computed in linear theory. It is related to the 
dimensionless matter power spectrum, P m , by 



dk 



r m (k) 



3MkR) 



kR 



(9) 



where R = 8/i~'Mpc and j\ is the spherical Bessel function of 
order 1. 

In addition, we compute £2 m /i 3 (a well-determined combina- 
tion orthogonal to the acoustic scale degeneracy in flat models; 
see e.g., Percival et al. 2002 and Howlett et al. 2012), 10 9 A s e ~ 2t 
(which determines the small-scale linear CMB anisotropy 
power), J"o.oo2 (the ratio of the tensor to primordial curvature 
power at k — 0.002 Mpc -1 ), Q. v h 2 (the physical density in mas- 
sive neutrinos), and the value of Yp from the BBN consistency 
condition. 

2.3. Likelihood 

Planck Collaboration XV (2013) describes the Planck tempera- 
ture likelihood in detail. Briefly, at high multipoles (C > 50) we 
use the 100, 143 and 217 GHz temperature maps to form a high 
multipole likelihood following the CamSpec methodology de- 
scribed in Planck Collaboration XV (2013). Apodized Galactic 
masks, including an apodized point source mask, are applied 
to individual detector/detector-set maps at each frequency. The 
masks are carefully chosen to limit contamination from diffuse 
Galactic emission to low levels (less than 20/xK 2 at all multi- 
poles). Thus we retain 57.8% of the sky at 100 GHz and 37.3% 
of the sky at 143 and 217 GHz. Mask deconvolved and beam 
corrected cross-spectra are computed for all detector/detector- 
set combinations and compressed to form averaged 100 x 100, 
143 x 143, 143 x 217 and 217 x 217 pseudo-spectra (note 
that we do not retain the 100 x 143 and 100 x 217 cross- 
spectra in the likelihood). Semi-analytic covariance matrices for 
these pseudo-spectra (Efstathiou 2004) are used to form a high- 
multipole likelihood in a fiducial Gaussian likelihood approxi- 
mation (Hamimeche & Lewis 2008). 

At low multipoles (2 < I < 49) the temperature likeli- 
hood is based on a Blackwell-Rao estimator applied to Gibbs 
samples computed by the Commander algorithm (Eriksen et al. 
2008) from Planck maps in the frequency range 30-353 GHz 
over 9 1 % of the sky. The likelihood at low multipoles therefore 
accounts for errors in foreground cleaning. 

Detailed consistency tests of both the high- and low- 
multipole components of the temperature likelihood are pre- 
sented in Planck Collaboration XV (2013). The high-multipole 
Planck likelihood requires a number of additional parame- 
ters to describe unresolved foreground components and other 
"nuisance" parameters (such as beam eigenmodes). The model 
adopted for Planck is described in Planck Collaboration XV 
(2013). A self-contained account is given in Sect. 4 which gen- 
eralizes the model to allow matching of the Planck likelihood to 



the likelihoods from high-resolution CMB experiments. A com- 
plete list of the foreground and nuisance parameters is given in 
Table 4. 



2.4. Sampling and confidence intervals 

We sample from the space of possible cosmological parame- 
ters with Markov Chain Monte Carlo (MCMC) exploration us- 
ing CosmoMC (Lewis & Bridle 2002). This uses a Metropolis- 
Hastings algorithm to generate chains of samples for a set of cos- 
mological parameters, and also allows for importance sampling 
of results to explore the impact of small changes in the analy- 
sis. The set of parameters is internally orthogonalized to allow 
efficient exploration of parameter degeneracies, and the baseline 
cosmological parameters are chosen following Kosowsky et al. 
(2002), so that the linear orthogonalisation allows efficient ex- 
ploration of the main geometric degeneracy (Bond et al. 1997). 
The code has been publicly available for a decade, has been 
thoroughly tested by the community, and has recently been ex- 
tended to sample efficiently large numbers of "fast" parameters 
by use of a speed-ordered Cholesky parameter rotation and a 
fast-parameter "dragging" scheme described by Neal (2005) and 
Lewis (2013). 

For our main cosmological parameter runs we execute eight 
chains until they are converged, and the tails of the distribu- 
tion are well enough explored for the confidence intervals for 
each parameter to be evaluated consistently in the last half of 
each chain. We check that the spread in the means between 
chains is small compared to the standard deviation, using the 
standard Gelman and Rubin (Gelman & Rubin 1992) criterion 
R - 1 < 0.01 in the least-converged orthogonalized parame- 
ter. This is sufficient for reliable importance sampling in most 
cases. We perform separate runs when the posterior volumes dif- 
fer enough that importance sampling is unreliable. Importance- 
sampled and extended data-combination chains used for this pa- 
per satisfy R— 1 < 0.1, and in almost all cases are closer to 0.01. 
We discard the first 30% of each chain as burn in, where the 
chains may be still converging and the sampling may be signif- 
icantly non-Markovian. This is due to the way CosmoMC learns 
an accurate orthogonalisation and proposal distribution for the 
parameters from the sample covariance of previous samples. 

From the samples, we generate estimates of the posterior 
mean of each parameter of interest, along with a confidence in- 
terval. We generally quote 68% limits in the case of two-tail 
limits, so that 32% of samples are outside the limit range, and 
there are 16% of samples in each tail. For parameters where the 
tails are significantly different shapes, we instead quote the inter- 
val between extremal points with approximately equal marginal- 
ized probability density. For parameters with prior bounds we ei- 
ther quote one-tail limits or no constraint, depending on whether 
the posterior is significantly non-zero at the prior boundary. Our 
one-tail limits are always 95% limits. For parameters with nearly 
symmetric distribution we sometimes quote the mean and stan- 
dard deviation (±l<x). The samples can also be used to estimate 
one, two and three-dimensional marginalized parameter posteri- 
ors. We use variable-width Gaussian kernel density estimates in 
all cases. 

We have also performed an alternative analysis to the one 
described above, using an independent statistical method based 
on frequentist profile likelihoods (Wilks 1938). This gives fits 
and error bars for the baseline cosmological parameters in ex- 
cellent agreement for both Planck and Planck combined with 
high-resolution CMB experiments, consistent with the Gaussian 
form of the posteriors found from full parameter space sampling. 
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In addition to posterior means, we also quote maximum- 
likelihood parameter values. These are generated using the 
BOBYQA bounded minimization routine 1 1 . Precision is limited by 
stability of the convergence, and values quoted are typically re- 
liable to within Ax 2 ~ 0.6, which is the same order as differ- 
ences arising from numerical errors in the theory calculation. 
For poorly constrained parameters the actual value of the best- 
fit parameters is not very numerically stable and should not be 
over-interpreted; in particular, highly degenerate parameters in 
extended models and the foreground model can give many ap- 
parently significantly different solutions within this level of ac- 
curacy. The best-fit values should be interpreted as giving typ- 
ical theory and foreground power spectra that fit the data well, 
but are generally non-unique at the numerical precision used; 
they are however generally significantly better fits than any of 
the samples in the parameter chains. Best-fit values are useful 
for assessing residuals, and differences between the best-fit and 
posterior means also help to give an indication of the effect of 
asymmetries, parameter-volume and prior-range effects on the 
posterior samples. We have cross-checked a small subset of the 
best-fits with the widely used MINUIT software (James 2004), 
which can give somewhat more stable results. 



3. Constraints on the parameters of the base 
ACDM model from Planck 

In this section we discuss parameter constraints from Planck 
alone in the ACDM model. Planck provides a precision mea- 
surement of seven acoustic peaks in the CMB temperature power 
spectrum. The range of scales probed by Planck is sufficiently 
large that many parameters can be determined accurately with- 
out using low-/' polarization information to constrain the optical 
depth, or indeed without using any other astrophysical data. 

However, because the data are reaching the limit of as- 
trophysical confusion, interpretation of the peaks at higher 
multipoles requires a reliable model for unresolved fore- 
grounds. We model these here parametrically, as described in 
Planck Collaboration XV (2013), and marginalize over the pa- 
rameters with wide priors. We give a detailed discussion of con- 
sistency of the foreground model in Sect. 4, making use of other 
high-/" CMB observations, although as we shall see the param- 
eters of the base ACDM model have a weak sensitivity to fore- 
grounds. 

As foreground modelling is not especially critical for the 
base ACDM model, we have decided to present the Planck con- 
straints early in this paper, ahead of the detailed descriptions of 
the foreground model, supplementary high-resolution CMB data 
sets, and additional astrophysical data sets. The reader can there- 
fore gain a feel for some of the key Planck results before being 
exposed to the lengthier discussions of Sects. 4 and 5, which are 
essential for the analysis of extensions to the base ACDM cos- 
mology presented in Sect. 6. 

In addition to the temperature power spectrum measurement, 
the Planck lensing reconstruction (discussed in more detail in 
Sect. 5.1 and Planck Collaboration XVII 2013) provides a dif- 
ferent probe of the perturbation amplitudes and geometry at late 
times. CMB lensing can break degeneracies inherent in the tem- 
perature data alone, especially the geometric degeneracy in non- 
flat models, providing a strong constraint on spatial curvature 
using only CMB data. The lensing reconstruction constrains the 



1 'http : // www . damtp . cam . ac . uk/user/na/NA_paper s/ 
NA2089_86.pdf 



matter fluctuation amplitude, and hence the accurate measure- 
ment of the temperature anisotropy power can be used together 
with the lensing reconstruction to infer the relative suppression 
of the temperature anisotropies due to the finite optical depth 
to reionization. The large-scale polarization from nine years of 
WMAP observations (Bennett et al. 2012) gives a constraint on 
the optical depth consistent with the Planck temperature and 
lensing spectra. Nevertheless, the WMAP polarization constraint 
is somewhat tighter, so by including it we can further improve 
constraints on some parameters. 

We therefore also consider the combination of the Planck 
temperature power spectrum with a WMAP polarization low- 
multipole likelihood (Bennett et al. 2012) at I < 23 (denoted 
WP), as discussed in Planck Collaboration XV (2013) 12 . We re- 
fer to this CMB data combination as Planck+WP. 

Table 2 summarizes our constraints on cosmological pa- 
rameters from the Planck temperature power spectrum alone 
(labelled "Planck"), from Planck in combination with Planck 
lensing (Planck+lensing) and with WMAP low-{ polariza- 
tion (Planck+WP). Figure 2 shows a selection of correspond- 
ing constraints on pairs of parameters and fully marginalized 
one-parameter constraints compared to the final results from 
WMAP (Bennett et al. 2012). 

3. 1 . Acoustic scale 

The characteristic angular size of the fluctuations in the CMB is 
called the acoustic scale. It is determined by the comoving size 
of the sound horizon at the time of last-scattering, r s (z*), and the 
angular diameter distance at which we are observing the fluc- 
tuations, D A (z*). With accurate measurement of seven acoustic 
peaks, Planck determines the observed angular size 9* = r s /D A 
to better than 0.1% precision at 1 o~. 

G m = (1.04148 ± 0.00066) x 10~ 2 = 0.596724° + 0.00038°. (10) 

Since this parameter is constrained by the positions of the peaks 
but not their amplitudes, it is quite robust; the measurement is 
very stable to changes in data combinations and the assumed 
cosmology. Foregrounds, beam uncertainties, or any system- 
atic effects which only contribute a smooth component to the 
observed spectrum will not substantially affect the frequency 
of the oscillations, and hence this determination is likely to 
be Planck's most robust precision measurement. The situation 
is analogous to baryon acoustic oscillations measurements in 
large-scale structure surveys (see Sect. 5.2), but the CMB acous- 
tic measurement has the advantage that it is based on observa- 
tions of the Universe when the fluctuations were very accurately 
linear, so second and higher-order effects are expected to be neg- 
ligible 13 . 

The tight constraint on 9* also implies tight constraints on 
some combinations of the cosmological parameters that deter- 
mine Da and r s . The sound horizon r s depends on the physical 

12 The WP likelihood is based on the WMAP likelihood module as 
distributed at http : / /lambda . gsf c . nasa . gov . 

13 Note, however, that Planck's measurement of 9, is now so accu- 
rate that 0(1CT 3 ) effects from aberration due to the relative motion be- 
tween our frame and the CMB rest-frame are becoming non-negligible; 
see Planck Collaboration XXVII (2013). The statistical anisotropy in- 
duced would lead to dipolar variations at the 10~ 3 level in ft, determined 
locally on small regions of the sky. For Planck, we average over many 
such regions and we expect that the residual effect (due to asymmetry 
in the Galactic mask) on the marginalised values of other parameters is 
negligible. 
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matter density parameters, and Da depends on the late-time evo- 
lution and geometry. Parameter combinations that fit the Planck 
data must be constrained to be close to a surface of constant 6*. 
This surface depends on the model that is assumed. For the base 
ACDM model, the main parameter dependence is approximately 
described by a 0.3% constraint in the three-dimensional Q m -/z- 
Cl\,h 2 subspace: 

Q m /! 3 - 2 (CVi 2 r°- 54 = 0.695 ± 0.002 (68%; Planck). (1 1) 

Reducing further to a two-dimensional subspace gives a 0.6% 
constraint on the combination 

£l m h 3 = 0.0959 + 0.0006 (68%; Planck). (12) 

The principle component analysis direction is actually Q m h 2 - 93 
but this is conveniently close to Q m /? 3 and gives a similar con- 
straint. The simple form is a coincidence of the ACDM cos- 
mology, error model, and particular parameter values of the 



model (Percival et al. 2002; Howie tt et al. 2012). The degener- 
acy between Ho and Q m is illustrated in Fig. 3: parameters are 
constrained to lie in a narrow strip where Q. m h 3 is nearly con- 
stant, but the orthogonal direction is much more poorly con- 
strained. The degeneracy direction involves consistent changes 
in the Hq, Q m , and Q.bh 2 parameters, so that the ratio of the sound 
horizon and angular diameter distance remains nearly constant. 
Changes in the density parameters, however, also have other 
effects on the power spectrum and the spectral index n s also 
changes to compensate. The degeneracy is not exact; its extent 
is much more sensitive to other details of the power spectrum 
shape. Additional data can help further to restrict the degeneracy. 
Figure 3 shows that adding WMAP polarization has almost no ef- 
fect on the Q.yjr' measurement, but shrinks the orthogonal direc- 
tion slightly from Q m /T 3 = 1 .03 ± 0. 1 3 to fi m /T 3 = 1.04 ± 0.11. 



10 



Planck Collaboration: Cosmological parameters 







Planck 


P/fl«eA:+lensing 


Planck+W 


Parameter 
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Best fit 68% limits 
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fin, 
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Ip 

Age/Gyr . 

z« 

r* 
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2drag - ■ ■ ■ 
''drag - ■ ■ ■ 

*D 

1OO0TJ . . . 



*--eq 

lOOSeq 

rd ral! /Ov(0.57) 



0.6825 
0.3175 
0.8344 
11.35 
67.11 
2.215 
0.14300 
0.09597 
0.247710 
13.819 
1090.43 
144.58 
1.04139 
1059.32 
147.34 
0.14026 
0.161332 
3402 
0.8128 
0.07130 



0.02207 ± 0.00033 
0.1 196 ±0.0031 
1.04132 ±0.00068 
0.097 ± 0.038 
0.9616 ± 0.0094 
3.103 ±0.072 

0.686 ± 0.020 
0.314 ±0.020 
0.834 ± 0.027 
1 1 4 +40 
67.4 ± 1.4 
2.23 ±0.16 
0.1423 ±0.0029 
0.09590 ± 0.00059 
0.24771 ±0.00014 
13.813 ±0.058 
1090.37 ± 0.65 
144.75 ± 0.66 
1.04148 ±0.00066 
1059.29 ± 0.65 
147.53 ±0.64 
0.14007 ±0.00064 
0.16137 ±0.00037 
3386 ± 69 
0.816 ±0.013 
0.0716 ±0.0011 



0.022242 
0.11805 
1.04150 
0.0949 
0.9675 
3.098 

0.6964 
0.3036 
0.8285 
11.45 
68.14 
2.215 
0.14094 
0.09603 
0.247785 
13.784 
1090.01 
145.02 
1.04164 
1059.59 
147.74 
0.13998 
0.161196 
3352 
0.8224 
0.07207 



0.02217 ±0.00033 
0.1 186 ±0.0031 
1.04141 ±0.00067 

0.089 ± 0.032 
0.9635 ± 0.0094 
3.085 ± 0.057 

0.693 ±0.019 
0.307 ± 0.019 
0.823 ±0.018 
10.81I;' 

67.9 ± 1.5 

2 19+ 012 

0.1414 ± 0.0029 
0.09593 ± 0.00058 
0.24775 ±0.00014 
13.796 ±0.058 
1090.16 ±0.65 
144.96 ± 0.66 
1.04156 ±0.00066 
1059.43 ± 0.64 
147.70 ± 0.63 
0.13996 ±0.00062 
0.16129 ±0.00036 
3362 ± 69 
0.821 ±0.013 
0.0719 ±0.0011 



0.022032 
0.12038 
1.04119 
0.0925 
0.9619 
3.0980 

0.6817 
0.3183 
0.8347 
11.37 
67.04 
2.215 
0.14305 
0.09591 
0.247695 
13.8242 
1090.48 
144.58 
1.04136 
1059.25 
147.36 
0.14022 
0.161375 
3403 
0.8125 
0.07126 



0.02205 ± 0.00028 
0.1 199 ±0.0027 
1.04131 ±0.00063 

0.9603 ± 0.0073 



0.829 ±0.012 
11.1 ± 1.1 
67.3 ± 1.2 
2 196+ 0051 
0.1426 ±0.0025 
0.09589 ± 0.00057 
0.24770 ±0.00012 
13.817 ±0.048 
1090.43 ± 0.54 
144.71 ±0.60 
1.04147 ±0.00062 
1059.25 ± 0.58 
147.49 ± 0.59 
0.14009 ±0.00063 
0.16140 ±0.00034 
3391 ± 60 
0.815 ±0.011 
0.07147 ±0.00091 



Table 2. Cosmological parameter values for the six-parameter base ACDM model. Columns 2 and 3 give results for the Planck 
temperature power spectrum data alone. Columns 4 and 5 combine the Planck temperature data with Planck lensing, and columns 
6 and 7 include WMAP polarization at low multipoles. We give best fit parameters as well as 68% confidence limits for constrained 
parameters. The first six parameters have flat priors. The remainder are derived parameters as discussed in Sect. 2. Beam, calibration 
parameters, and foreground parameters (see Sect. 4) are not listed for brevity. Constraints on foreground parameters for Planck+WP 
are given later in Table 5. 



3.2. Hubble parameter and dark energy density 

The Hubble constant, Hq, and matter density parameter, D. m , 
are only tightly constrained in the combination Cl m h 3 discussed 
above, but the extent of the degeneracy is limited by the effect 
of Q. m h 2 on the relative heights of the acoustic peaks. The pro- 
jection of the constraint ellipse shown in Fig. 3 onto the axes 
therefore yields useful marginalized constraints on H and Q m 
(or equivalently Qa) separately. We find the 2% constraint on 

H = (67.4 ± 14)kms"' Mpc" 1 (68%; Planck). (13) 

The corresponding constraint on the dark energy density param- 
eter is 

Q A = 0.686 + 0.020 (68%; Planck), (14) 

and for the physical matter density we find 

D. m h 2 =0.1423 + 0.0029 (68%; Planck). (15) 

Note that these indirect constraints are highly model depen- 
dent. The data only measure accurately the acoustic scale, and 



the relation to underlying expansion parameters (e.g., via the 
angular-diameter distance) depends on the assumed cosmology, 
including the shape of the primordial fluctuation spectrum. Even 
small changes in model assumptions can change Ho noticeably; 
for example, if we neglect the 0.06 eV neutrino mass expected 
in the minimal hierarchy, and instead take £ m v = 0, the Hubble 
parameter constraint shifts to 

H = (68.0 + 1.4) kms" 1 Mpc" 1 (68%; Planck, £ m y = 0). (16) 
3.3. Matter densities 

Planck can measure the matter densities in baryons and dark 
matter from the relative heights of the acoustic peaks. However, 
as discussed above, there is a partial degeneracy with the spec- 
tral index and other parameters that limits the precision of the 
determination. With Planck there are now enough well measured 
peaks that the extent of the degeneracy is limited, giving Q.\,h 2 to 
an accuracy of 1.5% without any additional data: 

Q h h 2 = 0.02207 + 0.00033 (68%; Planck). (17) 
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Fig. 3. Constraints in the Q m -Z/o plane. Points show samples 
from the Planck-only posterior, coloured by the corresponding 
value of the spectral index n s . The contours (68% and 95%) 
show the improved constraint from PZandr+lensing+WP The 
degeneracy direction is significantly shortened by including WP, 
but the well-constrained direction of constant Q. m h 3 (set by the 
acoustic scale), is determined almost equally accurately from 
Planck alone. 



large-scale polarization (Bennett et al. 2012). The large-scale E- 
mode polarization measurement is very challenging because it 
is a small signal relative to polarized Galactic emission on large 
scales, so this Planck polarization-free result is a valuable cross- 
check. The posterior for the Planck temperature power spectrum 
measurement alone also consistently peaks at t ~ 0.1, where the 
constraint on the optical depth is coming from the amplitude of 
the lensing smoothing effect and (to a lesser extent) the relative 
power between small and large scales. 

Since lensing constrains the underlying fluctuation ampli- 
tude, the matter density perturbation power is also well deter- 
mined: 



o-g = 0.823 + 0.018 (68%; Planck+lensing). 



(20) 



Much of the residual uncertainty is caused by the degeneracy 
with the optical depth. Since the small-scale temperature power 
spectrum more directly fixes <T%e~ T , this combination is tightly 
constrained: 



o- 8 e~ T = 0.753 ± 0.01 1 (68%; P/anc/t+lensing). 



(21) 



The estimate of <r% is significantly improved to <x 8 = 0.829 ± 
0.012 by using the WMAP polarization data to constrain the op- 
tical depth, and is not strongly degenerate with Q. m . (We shall 
see in Sect. 5.5 that the Planck results are discrepant with re- 
cent estimates of combinations of <x 8 and Q m from cosmic shear 
measurements and counts of rich clusters of galaxies.) 



Adding WMAP polarization information shrinks the errors by 
only 10%. 

The dark matter density is slightly less accurately measured 
at around 3%: 

Q c /i 2 = 0.1196 ± 0.0031 (68%; Planck). (18) 
3.4. Optical depth 

Small-scale fluctuations in the CMB are damped by Thomson 
scattering from free electrons produced at reionization. This 
scattering suppresses the amplitude of the acoustic peaks by e~ 2r 
on scales that correspond to perturbation modes with wavelength 
smaller than the Hubble radius at reionization. Planck measures 
the small-scale power spectrum with high precision, and hence 
accurately constrains the damped amplitude <?~ 2t A s . With only 
unlensed temperature power spectrum data, there is a large de- 
generacy between t and A s , which is weakly broken only by the 
power in large-scale modes that were still super-Hubble scale 
at reionization. However, lensing depends on the actual ampli- 
tude of the matter fluctuations along the line of sight. Planck 
accurately measures many acoustic peaks in the lensed tempera- 
ture power spectrum, where the amount of lensing smoothing de- 
pends on the fluctuation amplitude. Furthermore Planck's lens- 
ing potential reconstruction provides a more direct measurement 
of the amplitude, independently of the optical depth. The combi- 
nation of the temperature data and Planck's, lensing reconstruc- 
tion can therefore determine the optical depth t relatively well. 
The combination gives 

t = 0.089 + 0.032 (68%; P/ancfc+lensing). (19) 

As shown in Fig. 4 this provides marginal confirmation (just un- 
der 2cr) that the total optical depth is significantly higher than 
would be obtained from sudden reionization at z ~ 6, and is con- 
sistent with the WMAP-9 constraint, t = 0.089 + 0.014, from 



3.5. Spectral index 

The scalar spectral index defined in Eq. (2) is measured by 
Planck data alone to 1% accuracy: 

n s = 0.9616 ± 0.0094 (68%; Planck). (22) 

Since the optical depth t affects the relative power between large 
scales (that are unaffected by scattering at reionization) and in- 
termediate and small scales (that have their power suppressed 
by e~ 2r ), there is a partial degeneracy with « s . Breaking the de- 
generacy between t and n s using WMAP polarization leads to a 
small improvement in the constraint: 

n s = 0.9603 ± 0.0073 (68%; Planck+WP). (23) 

Comparing Eqs. (22) and (23), it is evident that the Planck tem- 
perature spectrum spans a wide enough range of multipoles to 
give a highly significant detection of a deviation of the scalar 
spectral index from exact scale invariance (at least in the base 
ACDM cosmology) independent of WMAP polarization infor- 
mation. 

One might worry that the spectral index parameter is degen- 
erate with foreground parameters, since these act to increase 
smoothly the amplitudes of the temperature power spectra at 
high multipoles. The spectral index is therefore liable to po- 
tential systematic errors if the foreground model is poorly con- 
strained. Figure 4 shows the marginalized constraints on the 
ACDM parameters for various combinations of data, includ- 
ing adding high-resolution CMB measurements. As will be dis- 
cussed in Sect. 4, the use of high-resolution CMB provides 
tighter constraints on the foreground parameters (particularly 
"minor" foreground components) than from Planck data alone. 
However, the small shifts in the means and widths of the distri- 
butions shown in Fig. 4 indicate that, for the base ACDM cos- 
mology, the errors on the cosmological parameters are not lim- 
ited by foreground uncertainties when considering Planck alone. 
The effects of foreground modelling assumptions and likelihood 
choices on constraints on n s are discussed in Appendix B. 
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Fig. 4. Marginalized constraints on parameters of the base ACDM model for various data combinations. 
Table 3. Summary of the CMB temperature data sets used in this analysis. 



Frequency Area Cin ttm S 

cut ^CMB ^tSZ ^Radio ^IR 

Experiment [GHz] [deg 2 ] [mjy] [GHz] [GHz] [GHz] [GHz] 

Planck 100 23846 50 1200 ... 100.0 103.1 

Planck 143 15378 50 2000 ... 143.0 145.1 ... 146.3 

Planck 217 15378 500 2500 ... 217.0 ... ... 225.7 

ACT(D13) 148 600 540 9440 15.0 148.4 146.9 147.6 149.7 

ACT(D13) 218 600 1540 9440 15.0 218.3 220.2 217.6 219.6 

SPT-high(R12) 95 800 2000 10000 6.4 95.0 97.6 95.3 97.9 

SPT-high(R12) 150 800 2000 10000 6.4 150.0 152.9 150.2 153.8 

SPT-high(R12) 220 800 2000 10000 6.4 220.0 218.1 214.1 219.6 

a Flux-density cut applied to the map by the point-source mask. For Planck the point-source mask is based on a composite of sources identified in 
the 100-353 GHz maps, so there is no simple flux cut. 



4. Planck combined with high-resolution CMB 
experiments: the base ACDM model 

The previous section adopted a foreground model with rela- 
tively loose priors on its parameters. As discussed there and 
in Planck Collaboration XV (2013), for the base ACDM model, 
the cosmological parameters are relatively weakly correlated 
with the parameters of the foreground model and so we ex- 
pect that the cosmological results reported in Sect. 3 are ro- 
bust. Fortunately, we can get an additional handle on unresolved 
foregrounds, particularly "minor" components such as the ki- 
netic SZ effect, by combining the Planck data with data from 
high-resolution CMB experiments. The consistency of results 
obtained with Planck data alone and Planck data combined with 
high-resolution CMB data gives added confidence to our cosmo- 
logical results, particularly when we come to investigate exten- 
sions to the base ACDM cosmology (Sect. 6). In this section, 
we review the high-resolution CMB data (hereafter, usually de- 
noted highL) that we combine with Planck and then discuss how 
the foreground model is adapted (with additional "nuisance" pa- 
rameters) to handle multiple CMB data sets. We then discuss the 
results of an MCMC analysis of the base ACDM model combin- 
ing Planck data with the high-^ data. 



4.1. Overview of the high-E CMB data sets 

The Atacama Cosmology Telescope (ACT) mapped the sky from 
2007 to 2010 in two distinct regions, the equatorial stripe (ACTe) 
along the celestial equator, and the southern stripe (ACTs) along 
declination -55°, observing in total about 600 deg 2 . The ACT 
data sets at 148 and 218 GHz are presented in Das et al. (2013, 
hereafter D13) and cover the angular scales 540 < I < 9440 at 
148 GHz and 1540 < I < 9440 at 218 GHz. Beam errors are 
included in the released covariance matrix. We include the ACT 
148 x 148 spectra for t > 1000, and the ACT 148 x218 and 21 8 x 
218 spectra for i > 1500. The inclusion of ACT spectra to I — 
1000 improves the accuracy of the inter-calibration parameters 
between the high-£ experiments and Planck. 

The South Pole Telescope observed a region of sky over 
the period 2007-10. Spectra are reported in Keisler et al. (201 1, 
hereafter Kll) and Story et al. (2012, hereafter S12) for angu- 
lar scales 650 < t < 3000 at 150 GHz, and in Reichardtet al. 
(2012, hereafter R12) for angular scales 2000 < € < 10000 at 
95, 150 and 220 GHz. Beam errors are included in the released 
covariance matrices used to form the SPT likelihood. The param- 
eters of the base ACDM cosmology derived from the WMAP- 
7+S12 data and (to a lesser extent) from Kll are in tension with 
Planck. Since the S12 spectra have provided the strongest CMB 
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constraints on cosmological parameters prior to Planck, this dis- 
crepancy merits a detailed analysis, which is presented in ap- 
pendix A. The S12 and Kll data are not used in combination 
with Planck in this paper. Since the primary purpose of includ- 
ing high-£ CMB data is to provide stronger constraints on fore- 
grounds, we use the R12 SPT data at i > 2000 only in combina- 
tion with Planck. We ignore any correlations between ACT/SPT 
and Planck spectra over the overlapping multipole ranges. 

Table 3 summarizes some key features of the CMB data sets 
used in this paper. 



4.2. Model of unresolved foregrounds and "nuisance" 
parameters 

The model for unresolved foregrounds used in the Planck 
likelihood is described in detail in Planck Collaboration XV 
(2013). Briefly, the model includes power spectrum templates 
for clustered extragalactic point sources (the cosmic infra-red 
background, hereafter CIB), thermal (tSZ) and kinetic (kSZ) 
Sunyaev-Zeldovich contributions, and the cross-correlation 
(tSZxCIB) between infra-red galaxies and the thermal Sunyaev- 
Zeldovich effect. The model also includes amplitudes for the 
Poisson contributions from radio and infra-red galaxies. The 
templates are described in Planck Collaboration XV (2013) and 
are kept fixed here. (Appendix B discusses briefly a few tests 
showing the impact of varying some aspects of the foreground 
model.) The model for unresolved foregrounds is similar to the 
models developed by the ACT and the SPT teams (e.g., R12; 
Dunkley et al. 2013). The main difference is in the treatment of 
the Poisson contribution from radio and infra-red galaxies. In 
the ACT and SPT analyses, spectral models are assumed for ra- 
dio and infra-red galaxies. The Poisson point source contribu- 
tions can then be described by an amplitude for each population, 
assuming either fixed spectral parameters or solving for them. 
In addition, one can add additional parameters to describe the 
decorrelation of the point source amplitudes with frequency (see 
e.g., Millea et al. 2012). The Planck model assumes free am- 
plitudes for the point sources at each frequency, together with 
appropriate correlation coefficients between frequencies. The 
model is adapted to handle the ACT and SPT data as discussed 
later in this section. 

Figure 5 illustrates the importance of unresolved foregrounds 
in interpreting the power spectra of the three CMB data sets. 
The upper panel of Fig. 5 shows the Planck temperature spec- 
tra at 100, 143, and 217 GHz, without corrections for unre- 
solved foregrounds (to avoid overcrowding, we have not plot- 
ted the 143 x 217 spectrum). The solid (red) lines show the 
best-fit base ACDM CMB spectrum corresponding to the com- 
bined Planck+ACT+SPT+WMAP polarization likelihood anal- 
ysis, with parameters listed in Table 5. The middle panel shows 
the SPT spectra at 95, 150 and 220 GHz from S12 and R12. 
In this figure, we have recalibrated the R12 power spectra to 
match Planck using calibration parameters derived from a full 
likelihood analysis of the base ACDM model. The S12 spec- 
trum plotted is exactly as tabulated in S12, i.e., we have not re- 
calibrated this spectrum to Planck. (The consistency of the S12 
spectrum with the theoretical model is discussed in further detail 
in Appendix A.) The lower panel of Fig. 5 shows the ACT spec- 
tra from D 13, recalibrated to Planck with calibration coefficients 
determined from a joint likelihood analysis. The power spectra 
plotted are an average of the ACTe and ACTs spectra, and in- 
clude the small Galactic dust corrections described in Das et al. 
(2013). 
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Fig. 5. Top: Planck spectra at 100, 143 and 217 GHz without 
subtraction of foregrounds. Middle: SPT spectra from R12 at 
95, 150 and 220 GHz, recalibrated to Planck using the best- 
fit calibration, as discussed in the text. The S12 SPT spec- 
trum at 150 GHz is also shown, but without any calibration cor- 
rection. This spectrum is discussed in detail in Appendix A, 
but is not used elsewhere in this paper. Bottom: ACT spectra 
(weighted averages of the equatorial and southern fields) from 
D13 at 148 and 220 GHz, and the 148x220 GHz cross-spectrum, 
with no extragalactic foreground corrections, recalibrated to the 
Planck spectra as discussed in the text. The solid line in each 
panel shows the best-fit base ACDM model from the combined 
Planck+WP+highL fits listed in Table 5. 
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The small-scale SPT (R12) and ACT (D13) data are domi- 
nated by the extragalactic foregrounds and hence are highly ef- 
fective in constraining the multi-parameter foreground model. 
In contrast, Planck has limited angular resolution and therefore 
limited ability to constrain unresolved foregrounds. Planck is 
sensitive to the Poisson point source contribution at each fre- 
quency and to the CIB contribution at 217 GHz. Planck has 
some limited sensitivity to the tSZ amplitude from the 100 GHz 
channel (and almost no sensitivity at 143 GHz). The remaining 
foreground contributions are poorly constrained by Planck and 
highly degenerate with each other in a Planck-alone analysis. 
The main gain in combining Planck with the high-resolution 
ACT and SPT data is in breaking some of the degeneracies be- 
tween foreground parameters which are poorly determined from 
Planck data alone. 

An important extension of the foreground parameterization 
described here over that developed in Planck Collaboration XV 
(2013) concerns the use of effective frequencies. Different exper- 
iments (and different detectors within a frequency band) have 
non-identical bandpasses (Planck Collaboration IX 2013) and 
this needs to be taken into account in the foreground mod- 
elling. Consider, for example, the amplitude of the CIB tem- 
plate at 217 GHz, A^ IB , introduced in Planck Collaboration XV 
(2013). The effective frequency for a dust-like component for 
the averaged 217 GHz spectrum used in the Planck likelihood 
is 225.7 GHz. To avoid cumbersome notation, we solve for the 
CIB amplitude A?J? at the CMB effective frequency of 217 GHz. 
The actual amplitude measured in the Planck 217 GHz band is 
1.33AS5, reflecting the different effective frequencies of a dust- 
like component compared to the blackbody primordial CMB 
(see Eq. 30 below). With appropriate effective frequencies, the 
single amplitude A?J? can be used to parameterize the CIB con- 
tributions to the ACT and SPT power spectra in their respective 
218 and 220 GHz bands. A similar methodology is applied to 
match the tSZ amplitudes for each experiment. 

The relevant effective frequencies for the foreground pa- 
rameterization discussed below are listed in Table 3. For the 
high resolution experiments, these are as quoted in R12 and 
Dunkley et al. (2013). For Planck these effective frequencies 
were computed from the individual HFI bandpass measure- 
ments (Planck Collaboration IX 2013), and vary by a few per- 
cent from detector to detector. The numbers quoted in Table 3 
are based on an approximate average of the individual de- 
tector bandpasses using the weighting scheme for individual 
detectors/detector-sets applied in the CamSpec likelihood. (The 
resulting bandpass correction factors for the tSZ and CIB ampli- 
tudes should be accurate to better than 5%.) 

The ingredients of the foreground model and associated 
"nuisance" parameters are summarized in the following para- 
graphs. 

Calibration factors To combine the Planck, ACT and SPT like- 
lihoods it is important to incorporate relative calibration factors, 
since the absolute calibrations of ACT and SPT have large errors 
(e.g., around 3.5% for the SPT 150 GHz channel). We introduce 
three map calibration parameters y^, y^I and y^o t0 rescale 
the R12 SPT spectra. These factors rescale the cross-spectra at 
frequencies v,- and vj as 



r v,xv, SPT SPTWfXvj 



(24) 



In the analysis of ACT, we solve for different map calibration 
factors for the ACTe and ACTs spectra, yf^ e , yf£? s , y^i^' m ^ 
y A f 8 Ts - In addition, we solve for the 100 x 100 and 217 x 217 



Planck power-spectrum calibration factors cioq and czn, with 
priors as described in Planck Collaboration XV (2013); see also 
Table 4. (The use of map calibration factors for ACT and SPT 
follows the conventions adopted by the ACT and SPT teams, 
while for the Planck power spectrum analysis we have consis- 
tently used power-spectrum calibration factors.) 

In a joint parameter analysis of P/a«cA:+ACT+SPT, the in- 
clusion of these calibration parameters leads to recalibrations 
that match the ACT, SPT and Planck 100 GHz and 217 GHz 
channels to the calibration of the Planck 143 x 143 spectrum 
(which, in turn, is linked to the calibration of the HFI 143-5 
detector, as described in Planck Collaboration XV 2013). It is 
worth mentioning here that the Planck 143 x 143 GHz spec- 
trum is 2.5% lower than the WMAP-9 combined V+W power 
spectrum (Hinshaw et al. 2012). This calibration offset between 
Planck HFI channels and WMAP is discussed in more detail in 
Planck Collaboration XI (2013) and in Appendix A. 

Poisson point source amplitudes To avoid any possible biases 
in modelling a mixed population of sources (synchrotron+dusty 
galaxies) with differing spectra, we solve for each of the Poisson 
point source amplitudes as free parameters. Thus, for Planck 
we solve for A™, A^ 3 , and A^ 7 , giving the amplitude of the 
Poisson point source contributions to ©3000 for the 100 x 100, 
143 x 143, and 217 x 217 spectra. The units of A™ are there- 
fore pK 2 . The Poisson point source contribution to the 143 x 217 
spectrum is expressed as a correlation coefficient, r^ 3x211 ; 



, n 143x217 _ PS [aPS a PS 

^3000 _ '143x217 \ /1 143 /1 217- 



(25) 



Note that we do not use the Planck 100 x 143 and 100 x 217 
spectra in the likelihood, and so we do not include correlation 
coefficients r^ jQxU3 or ^fo 0x2 i7- (These spectra carry little addi- 
tional information on the primordial CMB, but would require 
additional foreground parameters had we included them in the 
likelihood.) 

In an analogous way, the point source amplitudes for ACT 
and SPT are characterized by the amplitudes A?^ ACT , A^ S ' ACT 



.PS.SPT A PS.SPT , A PS.SPT ^ „ ;„ „„ ;to nfllV2 
A 95 ' A 150 • dnQ/1 220 



148 ' "217 

(all in units of pK z ) and three corre- 



lation coefficients r£f xl50 , r g^ x2 2^ an ^ r i50x220- T ne l ast °f these 
correlation coefficients is common to ACT and SPT. 



Kinetic SZ The kSZ template used here is from Trac et al. 
(201 1). We solve for the amplitude A ksz (in units of yuK 2 ): 



D) sz = A 



kSZ 



kSZ template 



kSZ template " 



(26) 



3000 



Thermal SZ We use the e = 0.5 tSZ template from 
Efstathiou & Migliaccio (2012) normalized to a frequency of 
143 GHz. 

For cross-spectra between frequencies v, and vj, the tSZ tem- 
plate is normalized as 

■*tSZ template 



tSZ.. 



= A 



.fivdf(vj) 'J>; 



143 f 2 (Vn) WSZ template ' 
-^3000 



(27) 



where vo is the reference frequency of 143 GHz, D ( 
the template spectrum at 143 GHz, and 



tSZ template 



/(V) = U 



+ 1 



with x 



hv 



(28) 
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The tSZ contribution is therefore characterized by the amplitude 
AfZ in units of fjK 2 . 

We neglect the tSZ contribution for any spectra involving the 
Planck 217 GHz, ACT 218 GHz, and SPT 220 GHz channels, 
since the tSZ effect has a null point at v = 217 GHz. (For Planck 
the bandpasses of the 217 GHz detectors see less than 0.1% of 
the 143 GHz tSZ power.) 



Cosmic infrared background The CIB contributions are ne- 
glected in the Planck 100 GHz and SPT 95 GHz bands and in 
any cross-spectra involving these frequencies. The CIB power 
spectra at higher frequencies are characterized by three ampli- 
tude parameters and a spectral index, 



= A 



^jCIB217x217 _ -JCIB 

-^7 _ ^217 



143 



CIB 



3000 



3000 



CIB 
'143x21 



I A CIB a CIB / j V 
7 V A 143 A 217^ 000 j 



(29a) 



(29b) 



(29c) 



where ArJ? and A™ are expressed in /xK 2 . As explained above, 
we define these amplitudes at the Planck CMB frequencies of 
143 and 217 GHz and compute scalings to adjust these ampli- 
tudes to the effective frequencies for a dust-like spectrum for 
each experiment. The scalings are 



cm Vj x 



CIB,,„ x „ n 

_ (£) y i0 x V<) 



/ g(vdg(vj) W vtVj r B Vi (T A ) B V] (T A ) 
3000 \g(vio)g(vjo)l\vioVjo) B Ya (T d ) B Vj0 (T d y 

(30) 

where ByiJi) is the Planck function at a frequency v, 

g(v) = [dB v (T)/dTT l \ TcuB (31) 

converts antenna temperature to thermodynamic temperature, v, 
and Vi refer to the Planck/ ACT/SPT dust effective frequencies, 
and v,o and vjq refer to the corresponding reference CMB Planck 
frequencies. In the analysis presented here, the parameters of 
the CIB spectrum are fixed to /3 A = 2.20 and T d = 9.7 K, as 
discussed in Addison et al. (2012a). The model of Eq. (30) then 
relates the Planck reference amplitudes of Eqs. (29b-29c) to 
the neighbouring Planck, ACT, and SPT effective frequencies, 
assuming that the CIB is perfectly correlated over these small 
frequency ranges. 

It has been common practice in recent CMB parameter 
studies to fix the slope of the CIB spectrum to y CIB = 0.8 
(e.g., Story et al. 2012; Dunkley et al. 2013). In fact, the shape 
of the CIB spectrum is poorly constrained at frequencies be- 
low 353 GHz and we have decided to reflect this uncertainty 
by allowing the slope y CIB to vary. We adopt a Gaussian 
prior on y CIB with a mean of 0.7 and a dispersion of 0.2. 
In reality, the CIB spectrum is likely to have some degree 
of curvature reflecting the transition between linear (two-halo) 
and non-linear (one-halo) clustering (see e.g., Cooray & Sheth 
2002; Planck Collaboration XVIII 2011; Amblard et al. 2011; 
Thackeretal. 2012). However, a single power law is an ad- 
equate approximation within the restricted multipole range 
(500 < I < 3000) over which the CIB contributes sig- 
nificantly to the P/a/icAyACT/SPT high-frequency spectra (as 
judged by the foreground-corrected power spectrum residuals 
shown in Figs. 7, 8 and 9 below). The prior on y CIB is moti- 
vated, in part, by the map-based Planck CIB analysis discussed 



in Planck Collaboration (2013). Appendix B explores different 
parameterizations of the CIB power spectrum. 

Thermal-SZ/CIB cross-correlation The cross-correlation be- 
tween dust emission from CIB galaxies and SZ emission from 
clusters (tSZxCIB) is expected to be non-zero. Because of un- 
certainties in the modelling of the CIB, it is difficult to compute 
this correlation with a high degree of precision. Addison et al. 
(2012b) present a halo-model approach to model this term and 
conclude that anti-correlations of around 10-20% are plausible 
between the clustered CIB components and the SZ at 150 GHz. 
The tSZxCIB correlation is therefore expected to make a mi- 
nor contribution to the unresolved foreground emission, but it is 
nevertheless worth including to determine how it might interact 
with other sub-dominant components, in particular the kSZ con- 
tribution. We use the Addison et al. (2012b) template spectrum 
in this paper and model the frequency dependence of the power 
spectrum as follows: 



tSZxCIB,. 



= -e 



tSZxCIB ^-.tSZxCIB template 



cm„ P 

3000 



3000 



CIB r . 
) ' 
3000 



"D 



tSZ„ iX 
3000 



,(32) 



where T> 



tSZxCIB template 



spectrum normalized to unity at 



is the Addison et al. (2012b) template 

CIB y . 



3000 and D 



tSZ,. 



and 



are given by Eqs. (27) and (31). The tSZxCIB con- 



tribution is therefore characterized by the dimensionless cross- 
correlation coefficient £ tSZxCIB . With the definition of Eq. (32), a 
positive value of £ tSZxCIB corresponds to an anti-correlation be- 
tween the CIB and the tSZ signals. 

Galactic cirrus For the masks used in the Planck CamSpec 
likelihood, Galactic dust makes a small contribution to 2D3000 
of around 5//K 2 to the 217 x 217 power spectrum, 1.5/iK 2 to 
the 143 x 217 spectrum, and around 0.5 fiK 2 to the 143 x 143 
spectrum. We subtract the Galactic dust contributions from 
these power spectra using a "universal" dust template spec- 
trum (at high multipoles this is accurately represented by a 
power law £)^ ust oc £~ 06 ). The template spectrum is based 
on an analysis of the 857 GHz Planck maps described in 
Planck Collaboration XV (2013), which uses mask-differenced 
power spectra to separate Galactic cirrus from an isotropic extra- 
galactic CIB contribution. This Galactic cirrus correction is kept 
fixed with an amplitude determined by template fitting the 217 
and 143 GHz Planck maps to the 857 GHz map, as described in 
Planck Collaboration XV (2013). Galactic dust contamination is 
ignored in the 100 x 100 spectrum. The Galactic dust template 
spectrum is actually a good fit to the dust contamination at low 
multipoles, ( <k 1000; however, we limit the effects of any in- 
accuracies in dust subtraction at low multipoles by truncating 
the 217 x 217 and 143 x 217 spectra at a minimum multipole 
of tmm = 500. (At multipoles ( % 1000, the Planck tempera- 
ture power spectra are signal dominated, so the 100 x 100 and 
143 x 143 spectra contain essentially all of the information on 
cosmology.) 

Compared to the contribution of Poisson point sources and 
the CIB, Galactic cirrus is a minor foreground component at 
217 GHz within our default mask, which retains 37% of the 
sky. However, the contribution of Galactic dust emission rises 
rapidly as more sky area is used. Extending the sky mask to 
65% of the sky (using the sequence of masks described in 
Planck Collaboration XV 2013), Galactic cirrus contributes to 
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Table 4. Astrophysical parameters used to model foregrounds in our analysis, plus instrumental calibration and beam parameters. 
We include the symbol for each parameter, the prior range adopted for the MCMC analysis and a summary definition (see text for 
details). Square brackets denote hard priors, parentheses indicate Gaussian priors. Note that the beam eigenmode amplitudes require 
a correlation matrix to describe fully their joint prior, and that all but /?[ are internally marginalized over rather than sampled over 
for the main MCMC runs. The bottom two blocks are only used in the analysis including the ACT and SPT high-f CMB data. 



Parameter 



Prior range 



Definition 



A PS 
^100 
4 PS 
A 143 
4 PS 
^217 
^S 

' 143x217 

a cm 

A 143 

4 cm 

^217 
„CIB 
'143x217 
yCIB 

A .sz 

^kSZ 
^tSZxCIB 

ClOO 
C217 



[-2 



[0, 360] 
[0, 270] 
[0,450] 

[0,1] 
[0, 20] 
[0, 80] 

[0,1] 
2] (0.7 ± 0.2) 
[0, 10] 
[0, 10] 

[0,1] 



Contribution of Poisson point-source power to D^jJ^ 100 for Planck (in /jK 2 ) 
As for A^,,, but at 143 GHz 
As forA^p, but at 217 GHz 

Point-source correlation coefficient for Planck between 143 and 217 GHz 

Contribution of CIB power to 2^ 43 * 143 at the Planck CMB frequency for 143 GHz (in fiK 2 ) 

As for A™ , but for 2 1 7 GHz 

CIB correlation coefficient between 143 and 217 GHz 

Spectral index of the CIB angular power (D c oc ^ yCIB ) 
Contribution of tSZ to £>£gg 143 at 143 GHz (in /M 2 ) 
Contribution of kSZ to D3000 (in A'K 2 ) 
Correlation coefficient between the CIB and tSZ (see text) 



[0.98, 1.02] (1.0006 ±0.0004) Relative power spectrum calibration for Planck between lOOGHzand 143GHz 
[0.95, 1.05] (0.9966 ± 0.0015) Relative power spectrum calibration for Planck between 217 GHz and 143 GHz 

(0 ± 1) Amplitude of the jth beam eigenmode (j = 1-5) for the rth cross-spectrum (1 = 1-4) 



"148 
A PS. ACT 
A 21B 

PS 
r l50x220 
4 ACTe 
'Must 
A ACTs 
^dusl 

ACTe 
^148 
..ACTe 
^217 

ACTs 
•M48 
..ACTs 
^217 



[0, 30] 
[0,200] 
[0,1] 
[0, 5] (0.8 ± 0.2) 
[0, 5] (0.4 ± 0.2) 
[0.8,1.3] 
[0.8, 1.3] 
[0.8,1.3] 
[0.8,1.3] 



Contribution of Poisson point-source power to D^m m f° r ACT (in /uK 2 ) 
As for A^j ACT , but at 218 GHz 

Point-source correlation coefficient between 150 and 220 GHz (for ACT and SPT) 
Contribution from Galactic cirrus to D3000 at 150 GHz for ACTe (in (iK 2 ) 
As A^ T <=, but for ACTs 

Map-level calibration of ACTe at 148 GHz relative to Planck 143 GHz 
Asyf£, Te ,butat217GHz 

Map-level calibration of ACTs at 148 GHz relative to Planck 143 GHz 
Asyff 8 T s jbut a t 217 GHz 



"95 

,PS. SPT 
A 150 
4 PS. SPT 
A 220 

PS 
'95x150 
r PS 

'95x220 

IS 

-'220 



[0, 30] 
[0, 30] 
[0,200] 
[0,1] 
[0,1] 
[0.8,1.3] 
[0.8, 1.3] 
[0.8,1.3] 



Contribution of Poisson point-source power to D^g 5 for SPT (in /jK 2 ) 
As for A P f SPT , but at 150 GHz 
As for A P f SPT , but at 220 GHz 

Point-source correlation coefficient between 95 and 150 GHz for SPT 
As r^| xl50 , but between 95 and 220 GHz 

Map-level calibration of SPT at 95 GHz relative to Planck 143 GHz 
As for y^, but at 150 GHz 
As for yS pT , but at 220 GHz 



D30Q0 around 50yuK 2 at 217 GHz (rising to around 200 ^K 2 on 
the scale of the first acoustic peak) and becomes a major fore- 
ground component, with an amplitude close to the net contri- 
bution of Poisson point sources and the clustered CIB. There is 
therefore a trade-off between limiting the signal-to-noise at 143 
and 217 GHz, by restricting the sky area, and potential system- 
atic errors associated with modelling Galactic cirrus over a large 
area of sky (i.e., sensitivity to the assumption of a "universal" 
dust template spectrum). We have chosen to be conservative in 
this first cosmological analysis of Planck by limiting the sky area 
at 143 and 217 GHz so that dust contamination is a minor fore- 
ground at high multipoles. As a further test of the importance 
of Galactic cirrus, we have analysed a Planck likelihood that re- 
tains only 24.7% of the sky (see Planck Collaboration XV 2013) 
at 217 GHz. Within this mask the CIB dominates over Galactic 
cirrus at multipoles I ~Z 500. There is a signal-to-noise penalty 
in using such a small area of sky at 217 GHz, but otherwise the 
results from this likelihood are in broad agreement with the re- 
sults presented here. With the conservative choices adopted in 



this paper, Galactic cirrus has no significant impact on our cos- 
mological results. 

We follow R12 and subtract a small-scale dust contribution 
of £>f st = 2.19jt/K 2 (^/3000r L2 from the R12 220 GHz spec- 
trum. This correction was determined by cross-correlating the 
SPT data with model 8 of Finkbeiner et al. (1999). For the ACT 
data we marginalize over a residual Galactic dust component 
Df st = A^J e/s (^/3000r - 7 , with different amplitudes for the 
southern and equatorial spectra, imposing Gaussian priors and 
frequency scaling as described in Dunkley et al. (2013). 

Notice that the spectral index of the SPT dust correction 
is significantly steeper than the dust correction applied to the 
Planck spectra. In future analyses it would be useful to derive 
more accurate dust corrections for the high-resolution CMB data 
by cross-correlating the SPT and ACT maps with the Planck 
545 and 857 GHz maps. Since the dust corrections are relatively 
small for the high-resolution data used here, we adopt the cor- 
rection described above in this paper. 

In application of the likelihood to Planck data alone, the 
model for unresolved foregrounds and relative calibrations con- 



17 



Planck Collaboration: Cosmological parameters 




Fig. 6. Comparison of the posterior distributions of 
the foreground parameters for Planck+WP (red) and 
P/a«c/t+ WP+highL (black). 



tains 13 parameters. In addition, we can solve for up to 20 beam 
eigenmode amplitudes (five amplitudes for each of the four spec- 
tra used in the Planck likelihood; see Planck Collaboration XV 
2013). In practice, we find that (usually) only the first beam 
eigenmode for the 100 x 100 spectrum, /?[, has a posterior dis- 
tribution that differs perceptibly from the prior, and we obtain 
nearly identical results on both foreground and cosmological pa- 
rameters if we treat only the amplitude of this eigenmode as 
a parameter and analytically marginalize over the rest. This is 
the default adopted in this paper. (The analytic marginalization 
improves stability of the minimisation for best-fit searches, and 
makes the Planck likelihood less cumbersome for the user.) 

The addition of ACT and SPT data introduces 17 extra pa- 
rameters. We provide a summary of the 50 foreground and nui- 
sance parameters in Table 4, including the prior ranges adopted 
in our MCMC analysis 14 . The choice of priors for many of these 
parameters is, to a large extent, subjective. They were chosen 
at an early stage in the Planck analysis to reflect "theoretically 
plausible" allowed ranges of the foreground parameters and to be 
broad compared to the results from high-resolution CMB exper- 
iments (which evolved over the course of this analysis as results 
from more ACT and SPT data were published). The foreground 
parameters from ACT and SPT depend on the assumptions of 
the underlying cosmology, and hence it is possible to introduce 
biases in the solutions for extensions to the base ACDM cosmol- 



14 Note that the foreground, calibration and beam parameters are all 
"fast" parameters as regards the MCMC sampling, and their inclusion 
has only a small impact on the computational speed. Marginalising over 
19 of the Planck beam parameters therefore leads only to 0(1) improve- 
ments in speed. 



ogy if overly restrictive foreground priors are imposed on the 
Planck data. Using the priors summarized in Table 4, the con- 
sistency between the Planck-alone results and the solutions for 
Planck combined with ACT and SPT provides a crude (but in- 
formative) measure of the sensitivity of cosmological results on 
the foreground model. Appendix B discusses the effects on ex- 
tended ACDM models of varying the priors on minor foreground 
components. 

4.3. The base ACDM model 

Cosmological and foreground parameters for the base six- 
parameter ACDM model are listed in Table 5, where here we 
give best-fit values and 68% confidence limits. The first two 
columns list the parameters derived from the Planck+WP anal- 
ysis discussed in Sect. 3, and are repeated here for easy refer- 
ence. The next two columns list the results of combining the 
Planck+WP likelihoods with the ACT and SPT likelihoods fol- 
lowing the model described above. We refer to this combination 
as "Planck+WP+highU' in this paper. The remaining columns 
list the parameter constraints combining the Planck+ WP+highL 
likelihood with the Planck lensing and BAO likelihoods (see 
Sect. 5). Table 5 lists the cosmological parameters for the base 
ACDM model and a selection of derived cosmological parame- 
ters. These parameters are remarkably stable for such data com- 
binations. We also list the values of the parameters describing the 
Planck foregrounds. A full list of all parameter values, includ- 
ing nuisance parameters, is given in the Explanatory Supplement 
(Planck Collaboration ES 2013). 

A comparison of the foreground parameter constraints from 
Planck+WP and P/ancfc+WP+highL is shown in Fig. 6; the 
corresponding cosmological parameter constraints are shown in 
Fig. 4. 

We can draw the following general conclusions. 

- The cosmological parameters for the base ACDM model are 
extremely insensitive to the foreground model described in 
the previous subsection. The addition of the ACT and SPT 
data causes the posterior distributions of cosmological pa- 
rameters to shift by much less than one standard deviation. 

- With Planck data alone, the CIB amplitude at 217 GHz is 
strongly degenerate with the 217 GHz Poisson point source 
amplitude. This degeneracy is broken by the addition of the 
high-resolution CMB data. This degeneracy must be borne 
in mind when interpreting Planck-only solutions for CIB pa- 
rameters; the sum of the Poisson point source and CIB con- 
tributions are well constrained by Planck at 217 GHz (and 
in good agreement with the map-based CIB Planck analy- 
sis reported in Planck Collaboration 2013), whereas the in- 
dividual contributions are not. Another feature of the CIB 
parameters is that we typically find smaller values of the CIB 
spectral index, y CIB , in Planck-alone solutions compared to 
Planck+highL solutions (which can be seen in Fig. 6). This 
provided additional motivation to treat y CIB as a parameter 
in the Planck likelihood rather than fixing it to a particu- 
lar value. There is evidence from the Planck spectra (most 
clearly seen by differencing the 217 x 217 and 143 X 143 
spectra) that the CIB spectrum at 217 GHz flattens in slope 
over the multipole range 500 < £ < 1000. This will be ex- 
plored in further detail in future papers (see also Appendix 
B). 

- The addition of the ACT and SPT data constrains the ther- 
mal SZ amplitude, which is poorly determined by Planck 
alone. In the Planck-alone analysis, the tSZ amplitude is 
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1 500 2000 



Fig. 7. Power spectrum residual plots, for assessing the quality of the foreground modelling. For each cross-spectrum, there are two 
sub-figures. The upper show the residuals with respect to the Planck+WP best-fit solution (from Table 5). The lower sub-figure 
show the residuals with respect to the Planck+ WP+highL solution The upper panel in each sub-figure shows the residual between 
the measured power spectrum and the best-fit (lensed) CMB power spectrum. The lower panels show the residuals after further 
removing the best-fit foreground model. The lines in the upper panels show the various foreground components. Major foreground 
components are shown by the solid lines, colour coded as follows: total foreground spectrum (red); Poisson point sources (orange); 
clustered CTB (blue); thermal SZ (green); and Galactic cirrus (purple). Minor foreground components are shown by the dotted 
lines colour coded as follows: kinetic SZ (green); tSZxCIB cross-correlation (purple). We also show residuals for the two spectra 
100 x 143 and 100 x 217 that are not used in the Planck likelihood. For these, we have assumed Poisson point-source correlation 
coefficients of unity. 
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I I 

Fig. 8. SPT power spectra at high I using the foreground model developed in this paper. The SPT R12 power spectra for each 
frequency combination are shown by the blue points, together with lcr error bars. The foreground components, determined from the 
Planck+ WP+highL analysis of ACDM models, are shown in the upper panels using the same colour coding as in Fig. 7. Here, the 
spectrum of the besft-fit CMB is shown in red and the total spectra are the upper green curves. The lower panel in each sub-figure 
shows the residuals with respect to the best-fit base ACDM cosmology +foreground model. The^- 2 values of the residuals, and the 
number of SPT bandpowers, are listed in the lower panels. 
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Fig. 9. As Fig. 8, but for the ACT south and ACT equatorial power spectra. 
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Planck+WP 


Hanot+WP+highL 


P/anrt+lensing+WP+highL 


P/anat+WP+highL+BAO 


Parameter 


Best fit 68% limits 


Best fit 68% limits 


Best fit 68% limits 


Best fit 68% limits 
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IOOAmc 



ln(10 lu A s ) . 
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/i 100 . . . . 
a PS 

A U3 . . . . 
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A t 

A 2 

A i 



1 1 i u 

M43 

a cm 



,.CIB 

' 143x217 



0.022032 0.02205 ± 0.00028 

0.12038 0.1199 ±0.0027 

1.04119 1.04131 ±0.00063 
0.0925 
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0.11847 0.1 186 ±0.0022 
1.04146 1.04144 ±0.00061 
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11.1 ± 1.1 
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11.1 ± 1.1 


11.42 


11.1 ± 1.1 


11.52 


11.3 ± 1.1 
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67.04 


67.3 ± 1.2 


67.15 


67.3 ± 1.2 


67.94 


67.9 ± 1.0 


67.77 


67.80 ± 0.77 
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13.798 ±0.037 


100fl» 


1.04136 
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1.04159 ±0.00060 
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1.04162 ±0.00056 
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Table 5. Best-fit values and 68% confidence limits for the base ACDM model. Beam and calibration parameters, and addi- 
tional nuisance parameters for "highL" data sets are not listed for brevity but may be found in the Explanatory Supplement 
(Planck Collaboration ES 2013). 



strongly degenerate with the Poisson point source ampli- 
tude at 100 GHz. This degeneracy is broken when the high- 
resolution CMB data are added to Planck. 



The last two points are demonstrated clearly in Fig. 7, which 
shows the residuals of the Planck spectra with respect to the 
best-fit cosmology for the Planck+WP analysis compared to the 
Planck+ WP+highL fits. The addition of high-resolution CMB 
data also strongly constrains the net contribution from the kSZ 
and tSZxCIB components (dotted lines), though these compo- 
nents are degenerate with each other (and tend to cancel). 

Although the foreground parameters for the Planck+WP fits 
can differ substantially from those for Planck+ WP+highL, the 
total foreground spectra are rather insensitive to the addition of 
the high-resolution CMB data. For example, for the 217x217 
spectrum, the differences in the total foreground solution are less 
than 10 /J.K 2 at t - 2500. The net residuals after subtracting both 
the foregrounds and CMB spectrum (shown in the lower panels 
of each sub-plot in Fig. 7) are similarly insensitive to the addi- 
tion of the high-resolution CMB data. The foreground model is 
sufficiently complex that it has a high "absorptive capacity" to 
any smoothly-varying frequency-dependent differences between 
spectra (including beam errors). 



Table 6. Goodness-of-fit tests for the Planck spectra. The Ay 2 = 
X 2 - N( is the difference from the mean assuming the model is 
correct, and the last column expresses Ay 2 in units of the disper- 
sion ^2N(. 



Spectrum 


p . 


r 

L max 


x 1 


X 2 INt 




100 x 100 


50 


1200 


1158 


1.01 


0.14 


143 x 143 


50 


2000 


1883 


0.97 


-1.09 


217x217 


500 


2500 


2079 


1.04 


1.23 


143x217 


500 


2500 


1930 


0.96 


-1.13 


All 


50 


2500 


2564 


1.05 


1.62 



To quantify the consistency of the model fits shown in Fig. 7 
for Planck we compute the \ 2 statistic 

x 2 = ^(Cf a - Cf MB - C?)M£(C$" - Cf MB - C% (33) 

tv 

for each of the spectra, where the sums extend over the mul- 
tipole ranges £ m i n and { lmx used in the likelihood, Aiec is 
the covariance matrix for the spectrum C^ ata (including cor- 
rections for beam eigenmodes and calibrations), Cf MB is the 
best-fit primordial CMB spectrum and C^ g is the best-fit fore- 
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ground model appropriate to the data spectrum. We expect 
X 1 to be approximately Gaussian distributed with a mean of 
Ne - { m ax - din + 1 and dispersion s/2N{. Results are sum- 
marized in Table 6 for the PZanc&:+ WP+highL best-fit parame- 
ters of Table 5. (The^ 2 values for the Planck+WP fit are almost 
identical.) Each of the spectra gives an acceptable global fit to 
the model, quantifying the high degree of consistency of these 
spectra described in Planck Collaboration XV (2013). (Note that 
Planck Collaboration XV 2013 presents an alternative way of in- 
vestigating consistency between these spectra via power spec- 
trum differences.) 

Figures 8 and 9 show the fits and residuals with respect to 
the best-fit P/ancfc+WP+highL model of Table 5, for each of the 
SPT and ACT spectra. The SPT and ACT spectra are reported 
as band-powers, with associated window functions [W^ PT ({)/{] 
and Wr CT (Q. The definitions of these window functions differ 
between the two experiments. 

For SPT, the contribution of the CMB and foreground spectra 
in each band is 



+ cf). 



(34) 



(Note that this differs from the equations give in R12 and S12.) 

For ACT, the window functions operate on the power spec- 
tra: 



c 6 = 2< CT w(Q CMB + cf). 



(35) 



In Fig. 9. we plot Db - tbi^b + l)Q,/(27r), where lb is the effec- 
tive multipole for band b. 

The upper panels of each of the sub-plots in Figs. 8 
and 9 show the spectra of the best-fit CMB, and the total 
CMB+foreground, as well as the individual contributions of the 
foreground components using the same colour codings as in 
Fig. 7. The lower panel in each sub-plot shows the residuals 
with respect to the best-fit cosmology+foreground model. For 
each spectrum, we list the value of x 2 , neglecting correlations 
between the (broad) ACT and SPT bands, together with the num- 
ber of data points. The quality of the fits is generally very good. 
For SPT, the residuals are very similar to those inferred from 
Fig. 3 of R 12. The SPT 150 x 220 spectrum has the largest x 1 
(approximately a 1.8<x excess). This spectrum shows systematic 
positive residuals of a few pK 2 over the entire multipole range. 
For ACT, the residuals and x 1 values are close to those plotted 
in Fig. 4 of Dunkley et al. (2013). All of the ACT spectra plot- 
ted in Fig. 9 are well fit by the model (except for some residuals 
at multipoles £ % 2000, which are also seen by Dunkley et al. 
2013). 

Having determined a solution for the best-fit foreground and 
other "nuisance" parameters, we can correct the four spectra 
used in the Planck likelihood and combine them to reconstruct a 
"best-fit" primary CMB spectrum and covariance matrix as de- 
scribed in Planck Collaboration XV (2013). This best-fit Planck 
CMB spectrum is plotted in the upper panels of Figs. 1 and 10 
for P/ancfc+ WP+highL foreground parameters. The spectrum in 
Fig. 10 has been band-averaged in bins of width AC ~ 31 using 
a window function Wb(l): 



&b = J] W b (€)i) { 



T>\-\ 



w h (0 



0, 



Anin — ^ ^ ^max' 

otherwise. 



(36a) 



(36b) 



Here, £ . and vL _ Y denote the minimum and maximum mul- 

' mm max 

tipole ranges of band b, and Al®, is the covariance ma- 
trix of the best-fit spectrum £)(, computed as described 
in Planck Collaboration XV (2013), and to which we have added 
corrections for beam and foreground errors (using the curvature 
matrix of the foreground model parameters from the MCMC 
chains). The solid lines in the upper panels of Figs. 1 and 10 
show the spectrum for the best-fit ACDM cosmology. The resid- 
uals with respect to this cosmology are plotted in the lower 
panel. To assess the goodness-of-fit, we compute^ 2 : 



X 



Cf MB ), 



(37) 



using the covariance matrix for the best-fit data spectrum (in- 
cluding foreground and beam errors' 3 ). The results are given in 
the last line of Table 6 labelled "All." The lower panel of Fig. 10 
shows the residuals with respect to the best-fit cosmology (on 
an expanded scale compared to Fig. 1). There are some visually 
striking residuals in this plot, particularly in the regions I ~ 800 
and I ~ 1300-1500 (where we see "oscillatory" behaviour). 
As discussed in detail in Planck Collaboration XV (2013), these 
residuals are reproducible to high accuracy across Planck de- 
tectors and across Planck frequencies; see also Fig. 7. There is 
therefore strong evidence that the residuals at these multipoles, 
which are in the largely signal dominated region of the spec- 
trum, are real features of the primordial CMB sky. These fea- 
tures are compatible with statistical fluctuations of a Gaussian 
ACDM model, and are described accurately by the covariance 
matrix used in the Planck likelihood. As judged by the x 2 statis- 
tic listed in Table 6, the best fit reconstructed Planck spectrum is 
compatible with the base ACDM cosmology to within 1.6cr 16 . 

To the extremely high accuracy afforded by the Planck data, 
the power spectrum at high multipoles is compatible with the 
predictions of the base six parameter ACDM cosmology. This is 
the main result of this paper. Fig. 1 does, however, suggest that 
the power spectrum of the best-fit base ACDM cosmology has 
a higher amplitude than the observed power spectrum at multi- 
poles ( < 30. We will return to this point in Sect. 7. 

Finally, Fig. 1 1 shows examples of Planck TE and EE spec- 
tra. These are computed by performing a straight average of the 
(scalar) beam-corrected 143 x 143, 143 x 217, and 217 x 217 
cross-spectra (ignoring auto-spectra). There are 32 TE and ET 
cross-spectra contributing to the mean TE spectrum plotted in 
Fig. 11, and six EE spectra contributing to the mean EE spec- 
trum. Planck polarization data, including LFI and 353 GHz data 
not shown here, will be analysed in detail, and incorporated 
into a Planck likelihood, following this data release. The pur- 
pose of presenting these figures here is twofold: first, to demon- 
strate the potential of Planck to deliver high quality polariza- 
tion maps and spectra, as described in the Planck "blue-book" 
(Planck Collaboration 2005); and, second, to show the consis- 
tency of these polarization spectra with the temperature spec- 
trum shown in Fig. 10. As discussed in Planck Collaboration VI 
(2013) and Planck Collaboration XV (2013), at present, the HFI 
polarization spectra at low multipoles (I < 200) are affected by 



15 Though the x~ value is similar if foreground and beam errors are 
not included in the covariance matrix. 

16 Planck Collaboration XXII (2013) describes a specific statistical 
test designed to find features in the primordial power spectrum. This test 
responds to the extended "dip" in the Planck power spectrum centred at 
about £ ~ 1800, tentatively suggesting 2.4-3. lcr evidence for a feature. 
A more extensive analysis is required to assess whether the statistical 
significance assigned to this result is accurate. 
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Fig. 10. Planck TT power spectrum. The points in the upper panel show the maximum-likelihood estimates of the primary CMB 
spectrum computed as described in the text for the best-fit foreground and nuisance parameters of the PZanc£+WP+highL fit listed 
in Table 5. The red line shows the best-fit base ACDM spectrum. The lower panel shows the residuals with respect to the theoretical 
model. The error bars are computed from the full covariance matrix, appropriately weighted across each band (see Eqs. 36a and 
36b), and include beam uncertainties and uncertainties in the foreground model parameters. 
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Fig. 11. Planck TE (left) and EE spectra (right) computed as described in the text. The red lines show the polarization spectra from 
the base ACDM PZanc£+ WP+highL model, which is fitted to the TT data only. 
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systematic errors that cause biases. For the HFI channels used in 
Fig. 1 1, there are two primary sources of systematic error arising 
from non-linear gain-like variations, and residual bandpass mis- 
matches between detectors. However, these systematics rapidly 
become unimportant at higher multipoles 17 . 

The errors on the mean TE and EE spectra shown in Fig. 
1 1 are computed from the analytic formulae given in Efstathiou 
(2006), using an effective beam-width adjusted to reproduce the 
observed scatter in the polarization spectra at high multipoles. 
The spectra are then band-averaged as in Eq. (37). The error bars 
shown in Fig. 1 1 are computed from the diagonal components of 
the band-averaged covariance matrices. 

The solid lines in the upper panels of Fig. 11 show 
the theoretical TE and EE spectra expected in the best-fit 
Planck+ WP+highL ACDM model (i.e., the model used to com- 
pute the theory TT spectrum plotted in Fig. 10). These theoret- 
ical spectra are determined entirely from the TT analysis and 
make no use of the Planck polarization data. As with the TT 
spectra, the ACDM model provides an extremely good match 
to the polarization spectra. Furthermore, polarized foreground 
emission is expected to be unimportant at high multipoles (e.g., 
Tucci & Toffolatti 2012) and so no foreground corrections have 
been made to the spectra in Fig. 11. The agreement between 
the polarization spectra and the theoretical spectra therefore pro- 
vides strong evidence that the best-fit cosmological parameters 
listed in Table 5 are not strongly affected by the modelling of 
unresolved foregrounds in the TT analysis. 



5. Comparison of the Planck base ACDM model 
with other astrophysical data sets 

Unlike CMB data, traditional astrophysical data sets - e.g., 
measurements of the Hubble parameter, type la supernovae 
(SNe la), and galaxy redshift surveys - involve complex phys- 
ical systems that are not understood at a fundamental level. 
Astronomers are therefore reliant on internal consistency tests 
and empirical calibrations to limit the possible impact of system- 
atic effects. Examples include calibrating the metallicity depen- 
dence of the Cepheid period luminosity relation, calibrating the 
colour-decline-rate-luminosity relation of Type la supernovae, 
or quantifying the relationship between the spatial distributions 
of galaxies and dark matter. In addition, there are more mundane 
potential sources of error, which can affect certain types of astro- 
physical observations (e.g., establishing consistent photometric 
calibration systems). We must be open to the possibility that un- 
known, or poorly quantified, systematic errors may be present 
in the astrophysical data, especially when used in combination 
with the high precision data from Planck. 

We have seen in the previous section that the base ACDM 
model provides an acceptable fit to the Planck TT power spec- 
tra (and the Planck TE and EE spectra) and also to the ACT 
and SPT temperature power spectra. The cosmological parame- 
ters of this model are determined to high precision. We therefore 
review whether these parameters provide acceptable fits to other 
astrophysical data. If they do not, then we need to assess whether 
the discrepancy is a pointer to new physics, or evidence of some 
type of poorly understood systematic effect. Unless stated other- 
wise, we use the Planck+ WP+highL parameters listed in Table 5 
as the default "Planck" parameters for the base ACDM model. 
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"The main focus of current work on Planck polarization is to reduce 
the effects of these systematics on the polarization maps at large angular 
scales. 
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I 

Fig. 12. Planck measurements of the lensing power spec- 
trum compared to the prediction for the best-fitting 
Planck+ WP+highL ACDM model parameters. In the top 
panel, the data points are the measured bandpowers and ±1 cr 
error ranges from the diagonal of the covariance matrix. The 
measured bandpowers are compared to the Cy in the best-fit 

model (black line). The grey region shows the 1 cr range in Cy 
due to ACDM parameter uncertainties. The lower panel shows 
the differences between the bandpower amplitudes A, and the 
predictions for their expectation values in the best-fit model, 

. theory 



5.1. CMB lensing measured by Planck 

Weak gravitational lensing by large-scale structure subtly al- 
ters the statistics of the CMB anisotropies, encoding informa- 
tion about the late-time Universe which is otherwise degen- 
erate in the primary anisotropies laid down at last-scattering 
(see Lewis & Challinor 2006 for a review). The lensing deflec- 
tions are given by the gradient of the lensing potential <p{h), 
which corresponds to an integrated measure of the matter distri- 
bution along the line of sight with peak sensitivity to structures 
around redshift 2. The rms deflection is expected to be around 
2.5 arcmin and to be coherent over several degrees. We include 
the effect of lensing on the temperature power spectrum in all our 
parameter analysis, but for some results we also include the lens- 
ing information encoded in the non-Gaussian trispectrum (con- 
nected 4-point function) of the CMB. Lensing generates a non- 
zero trispectrum, which, at leading order, is proportional to the 
power spectrum Cy* of the lensing potential (Hu 2001). 

In Planck Collaboration XVII (2013), we present a detailed 
analysis of CMB lensing with Planck data, including estima- 
tion of Cy from the trispectrum computed from Planck's maps. 
This paper also describes the construction of a lensing like- 
lihood. Briefly, we first reconstruct an estimate of the lens- 
ing potential using near-optimal quadratic estimators, follow- 
ing Okamoto & Hu (2003), with various Galactic and point- 
source masks. The empirical power spectrum of this reconstruc- 
tion, after subtraction of the Gaussian noise bias (i.e., the dis- 
connected part of the 4-point function), is then used to esti- 
mate Cy in bandpowers. The associated bandpower errors are 
estimated from simulations. The lensing power spectrum is es- 
timated from channel-coadded Planck maps at 100, 143 and 
217 GHz in the multipole range t = 10-1000, and also from a 
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minimum-variance combination of the 143 and 217 GHz maps. 
An empirical correction for the shot-noise trispectrum of un- 
resolved point sources is made to each spectrum, based on the 
measured amplitude of a generalized kurtosis of the appropriate 
maps. Additionally, the A^ (1) bias of Kesden et al. (2003), com- 
puted for a fiducial ACDM spectrum determined from a pre- 
publication analysis of the Planck data, is subtracted from each 
spectrum. This latter correction is proportional to Cy and ac- 
counts for sub-dominant couplings of the trispectrum, which mix 
lensing power over a range of scales into the power spectrum 
estimates. Excellent internal consistency of the various Cy esti- 
mates is found over the full multipole range. 

The Planck lensing likelihood is based on reconstruc- 
tions from the minimum-variance combination of the 143 and 
217 GHz maps with 30% of the sky masked. Conservatively, 
only multipoles in the range ( = 40-400 are included, with a 
bandpower width M = 45. The range £ = 40^-00 captures 90% 
of the signal-to-noise on a measurement of the amplitude of a 
fiducial Cy* , while minimizing the impact of imperfections in 
modelling the effect of survey anisotropies on the large-scale <p 
reconstruction (the "mean-field" of Planck Collaboration XVII 
2013), and the large Gaussian noise bias on small scales. Note, 
however, that by restricting the range of angular scales we do 
lose some ability to distinguish between scale-dependent mod- 
ifications of Cy, such as from massive neutrinos, and almost 
scale-independent modifications, such as from changes in the 
equation of state of unclustered dark energy or spatial curva- 
ture. Correlated uncertainties in the beam transfer functions, 
point-source corrections, and the cosmology dependence of the 
N m bias, are modelled in the covariance matrix. The resulting 
bandpower correlations are all small, less than 4%, but broad. 
Bandpower correlations induced by masking are estimated to be 
less than 5% for neighbouring bins and are neglected. The likeli- 
hood is modelled as a Gaussian in the bandpowers with a fiducial 
(i.e., parameter-independent) covariance. For verification of this 
approximation, see Schmittfull et al. (2013). 

The connected four-point function is related to the fully- 
reduced trispectrum Tff(L) by 



(T, 



LM 

£3 £\ L 
mi ni4 —M 
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mi ni2 M 



T/y(L) + perms, (38) 



(Hu 2001). In the context of lensing reconstruction, the CMB 
trispectrum due to lensing takes the form 



TT s-^TT , 



(39) 



where C T ( T is the lensed temperature power spectrum and F( { l( 2 
is a geometric mode-coupling function (Hu 2001; Hanson et al. 
2011). Our estimates of Cy derive from the measured trispec- 
trum. They are normalized using the fiducial lensed power spec- 
trum to account for the factors of C T ( T in Eq. (39). In the like- 
lihood, we renormalize the parameter-dependent Cy to account 
for the mismatch between the parameter-dependent Cj T and that 
in the fiducial model. Since the best-fit ACDM model we con- 
sider in this section has a lensed temperature power spectrum 
that is very close to that of the fiducial model, the renormalisa- 
tion factor differs from unity by less than 0.25%. 

The estimated lensing power spectrum Cy is not indepen- 
dent of the measured temperature power spectrum Cj T , but the 



dependence is very weak for Planck, and can be accurately 
ignored (Schmittfull et al. 2013; Planck Collaboration XVII 
2013). As discussed in detail in Schmittfull et al. (2013), there 
are several effects to consider. First, the reconstruction noise in 
the estimated <p derives from chance correlations in the unlensed 
CMB. If, due to cosmic variance, the unlensed CMB fluctuates 
high at some scale, the noise in the reconstruction will gen- 
erally increase over a broad range of scales. Over the scales 
relevant for Planck lensing reconstruction, the correlation be- 
tween the measured Cy and C r J from this effect is less than 
0.2% and, moreover, is removed by a data-dependent Gaussian 
noise bias removal that we adopt following Hanson et al. (201 1) 
and Namikawa et al. (2012). The second effect derives from cos- 
mic variance of the lenses. If a lens on a given scale fluctu- 
ates high, the estimated Cy will fluctuate high at that scale. 
In tandem, there will be more smoothing of the acoustic peaks 
in the measured C r J , giving broad-band correlations that are 
negative at acoustic peaks and positive at troughs. The max- 
imum correlation is around 0.05%. If we consider estimating 
the amplitude of a fiducial lensing power spectrum indepen- 
dently from the smoothing effect of Cj T and the measured Cy 
in the range £ = 40^-00, the correlation between these esti- 
mates due to the cosmic variance of the lenses is only 4%. This 
amounts to a mis-estimation of the error on a lensing ampli- 
tude in a joint analysis of Cy and C T t T , treated as indepen- 
dent, of only 2%. For physical parameters, the mis-estimation 
of the errors is even smaller: Schmittfull et al. (2013) estimate 
around 0.5% from a Fisher analysis. A third negligible effect 
is due to the T— <f> correlation sourced by the late integrated 
Sachs-Wolfe effect (see Planck Collaboration XIX 2013). This 
produces only local correlations between the measured Cy and 
Cj T which are less than 0.5% by £ = 40 and fall rapidly on 
smaller scales. They produce a negligible correlation between 
lensing amplitude estimates for the multipole ranges considered 
here. The T-<p correlation is potentially a powerful probe of 
dark energy dynamics (e.g., Verde & Spergel 2002) and modi- 
fied theories of gravity (e.g., Acquaviva et al. 2004). The power 
spectrum C ( * can be measured from the Planck data using the 
CMB 3-point function (Planck Collaboration XXIV 2013) or, 
equivalently, by cross-correlating the <p reconstruction with the 
large-angle temperature anisotropies Planck Collaboration XIX 
(2013) although the detection significance is only around 3cr. 
The power-spectrum based analysis in this paper discards the 
small amount of information in the T-<p correlation from Planck. 
In summary, we can safely treat the measured temperature and 
lensing power spectra as independent and simply multiply their 
respective likelihoods in a joint analysis. 

We note that ACT (Das et al. 2011, 2013) and 
SPT (van Engelen et al. 2012) have both measured the lensing 
power spectrum with significances of 4.6 cr and 6.3 cr, respec- 
tively, in the multipole ranges t = 75-2050 and i = 100-1500. 
The Planck measurements used here represent a 26 cr detection. 
We therefore do not expect the published lensing measurements 
from these other experiments to carry much statistical weight in 
a joint analysis with Planck, despite the complementary range 
of angular scales probed, and we choose not to include them in 
the analyses in this paper. 

In the lensing likelihood, we characterize the estimates of 
Cy* with a set of eight (dimensionless) amplitudes A,, where 



(40) 
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Table 7. Planck CMB lensing constraints. The A. eory are renor- 
malized power spectrum amplitudes in the best-fit ACDM model 
to Planck+WP+highL, and the errors cr(A,) are the square root 
of the diagonals of the A, covariance matrix. 
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0.08 
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0.77 


0.98 


0.10 


5 


220 


264 


0.88 


0.98 


0.13 


6 


265 


309 


0.88 


0.98 


0.16 


7 


310 


354 


1.10 


0.98 


0.18 


8 


355 


400 


0.75 


0.98 


0.21 



Here, S c t is a binning operation with 

Bf = — " —j . (41) 

for I within the band defined by a minimum multipole f' min and a 
maximum C l mm . The inverse of the weighting function, Vc, is an 
approximation to the variance of the measured Cy* and Cy*' fid 
is the lensing power spectrum of the fiducial model, which is 
used throughout the analysis. The A,- are therefore near-optimal 
estimates of the amplitude of the fiducial power spectrum within 
the appropriate multipole range, normalized to unity in the fidu- 
cial model. Given some parameter-dependent model Cy, the ex- 
pected values of the A, are 

<A,> = Af°* = J] Bf [l + A*(Cf )] 2 , (42) 

i 

where the term involving A^iCj 7 ), which depends on the 
parameter-dependent C T t T , accounts for the renormalisation 
step described above. The lensing amplitudes A, are com- 
pared to the ,4* eory f or the best-fitting ACDM model to the 
Planck+ WP+highL data combination (i.e., not including the 
lensing likelihood) in Table 7. The differences between A,- and 
^theoiy are pj otte( j m tne bottom panel of Fig. 12 while in the 

top panel the bandpower estimates are compared to Cy in the 

best-fitting model. The Planck measurements of Cy are con- 
sistent with the prediction from the best-fit ACDM model to 
Planck+ WP+highL. Using the full covariance matrix, we find 
X 1 = 10.9 with eight degrees of freedom, giving an accept- 
able probability to exceed of approximately 21%. It is worth 
recalling here that the parameters of the ACDM model are 
tightly constrained by the CMB 2-point function (as probed by 
our Planck+ WP+highL data combination) which derives from 
physics at z » 1 100 seen in angular projection. It is a significant 
further vindication of the ACDM model that its predictions for 
tlie evolution of structure and geometry at much lower redshifts 
(around z — 2) fit so well with Planck's CMB lensing measure- 
ments. 

The discussion above does not account for the small spread 
in the Cy predictions across the Planck+ WP+highL ACDM 
posterior distribution. To address this, we introduce a parameter 
A^ which, at any point in parameter space, scales the lensing 
trispectrum. Note that Ay* does not alter the lensed temperature 




Amplitude relative to physical 



Fig. 13. Marginalized posterior distributions for Ay 

(dashed) and Al (solid). For Ay we use the data combi- 
nation P/anc£+lensing+WP+WghL. For Al we consider 
P/a«cfe+lensing+ WP+highL (red), Planck+ WP+highL (green), 
Planck+WP (blue) and P/a«c£-lowL+highL+Tprior (cyan; see 
text). 

power spectrum, so it can be used to assess directly how well the 
ACDM predictions from Cj T agree with the lensing measure- 
ments; in ACDM we have Ay* = 1. The marginalized posterior 

distribution for A^ in a joint analysis of Planck+ WP+highL and 
the Planck lensing likelihood is given in Fig. 13. The agreement 
with Ay* — 1 is excellent, with 

A M> _ o 99 ± o 05 (68% . p/ a „ c/ t+lensing+WP+highL). (43) 

The significance of the detection of lensing using Ay* in ACDM 
is a little less than the 26 cr detection of lensing power reported 
in Planck Collaboration XVII (2013), due to the small spread in 
Cy from ACDM parameter uncertainties. 

Lensing also affects the temperature power spectrum, pri- 
marily by smoothing the acoustic peaks and troughs on the scales 
relevant for Planck. The most significant detection of the lensing 
effect in the power spectrum to date is from SPT Introducing a 
parameter Al (Calabrese et al. 2008) which takes C, — > A\jCy 
when computing the lensed temperature power spectrum (we 
shall shortly extend the action of this parameter to include the 
computation of the lensing trispectrum), Story et al. (2012) re- 
port A L = 0.86+° -\l (68%; SPT+WMAP-7). Results for A L from 
Planck in combination with WMAP low-£ polarization and the 
high-^ power spectra from ACT and SPT are also shown in 
Fig. 13. Where we include the Planck lensing measurements, 
we define Al to scale the explicit Cy' in Eq. (39), as well as 
modulating the lensing effect in the temperature power spec- 
trum. Figure 13 reveals a preference for Al > 1 from the Planck 
temperature power spectrum (plus WMAP polarization). This is 
most significant when combining with the high-^ experiments 
for which we find 

A L = 1.23 ± 0.11 (68%; Planck+ WP+highL), (44) 

i.e., a 2 cr preference for Al > 1. Including the lensing measure- 
ments, the posterior narrows but shifts to lower Al, becoming 
consistent with Al = 1 at the 1 cr level as expected from the Ay* 
results. 
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Fig. 14. Effect of allowing Al to vary on the degeneracies between Q^h 1 and n s (left) and Q. m h 2 and n s (right). In both panels the 
data combination is P/anc£+ WP+highL. The contours enclose the 68% and 95% confidence regions in the base ACDM model with 
Al = 1. The samples are from models with variable Al and are colour-coded by the value of Al. 



We do not yet have a full understanding of what is driving 
the preference for high Al in the temperature power spectrum. 
As discussed in Appendix B, the general preference is stable to 
assumptions about foreground modelling and cuts of the Planck 
data in the likelihood. To gain some insight, we consider the 
range of multipoles that drive the preference for Al > 1 . For our 
favoured data combination of PZancA:+ WP+highL, Ax 2 = -4.6 
going from the best-fit Al = 1 model to the best-fit model with 
variable Al- The improvement in fit comes only from the low-/' 
temperature power spectrum (Ay 2 = —1.3) and the ACT+SPT 
data (Ax 2 = -3.5); for this data combination, there is no prefer- 
ence for high Al from the Planck temperature data at intermedi- 
ate and high multipoles (A^ 2 = +0.15). The situation at low-{ is 
similar if we exclude the high-^ experiments, with Ax 1 = - 1.9 
there, but there is then a preference for the high Al best-fit from 
the Planck data on intermediate and small scales (Ay 2 = -3.6). 
However, as discussed in Sect. 4, there is more freedom in the 
foreground model when we exclude the high-^ data, and this can 
offset smooth differences in the CMB power spectra such as the 
transfer of power from large to small scales by lensing that is 
enhanced for Al > 1 . 

Since the low-f temperature data seem to be partly responsi- 
ble for pulling A L high, we consider the effect of removing the 
low-f likelihood from the analysis. In doing so, we also remove 
the WMAP large-angle polarization which we compensate by in- 
troducing a simple prior on the optical depth; we use a Gaussian 
with mean 0.09 and standard deviation 0.013, similar to the con- 
straint from WMAP polarization (Hinshaw et al. 2012). We de- 
note this data combination, including the high-^ experiments, by 
P/a«c£-lowL+highL+rprior and show the posterior for Al in 
Fig. 13. As anticipated, the peak of the posterior moves to lower 
A L giving A L = 1-17^'j' (68% CL) and the A^ 2 between the 
best-fit model (now at Al = 1.13) and the Al = 1 model remains 
negligible for the Planck data. 

Since varying Al alone does not alter the power spectrum 
on large scales, why should the low-/' data prefer higher A L ? 
The reason is due to a chain of parameter degeneracies that 
are illustrated in Fig. 14, and the deficit of power in the mea- 
sured C[ s on large scales compared to the best-fit ACDM model 
(see Fig. 1 and Sect. 7). In models with a power-law primor- 
dial spectrum, the temperature power spectrum on large scales 



can be reduced by increasing n s . The effect of an increase in 
rc s on the relative heights of the first few acoustic peaks can 
be compensated by increasing a>\, and reducing a> m , as shown 
by the contours in Fig. 14. However, on smaller scales, corre- 
sponding to modes that entered the sound horizon well before 
matter-radiation equality, the effects of baryons on the mid-point 
of the acoustic oscillations (which modulates the relative heights 
of even and odd peaks) is diminished since the gravitational po- 
tentials have pressure-damped away during the oscillations in 
the radiation-dominated phase (e.g., Hu & White 1996, 1997). 
Moreover, on such scales the radiation-driving at the onset of the 
oscillations that amplifies their amplitude happens early enough 
to be unaffected by small changes in the matter density. The net 
effect is that, in models with Al = 1, the extent of the degeneracy 
involving n s , a)\, and u> m is limited by the higher-order acoustic 
peaks, and there is little freedom to lower the large-scale tem- 
perature power spectrum by increasing n s while preserving the 
good fit at intermediate and small scales. Allowing Al to vary 
changes this picture, letting the degeneracy to extend to higher 
n s , as shown by the samples in Fig. 14. The additional smooth- 
ing of the acoustic peaks due to an increase Al can mitigate the 
effect of increasing n s around the fifth peak, where the signal- 
to-noise for Planck is still high. This allows one to decrease 
the spectrum at low £, while leaving it essentially unchanged 
on those smaller scales where Planck still has good sensitiv- 
ity. Above i ~ 2000, the best-fit Al model has a little more 
power than the base model (around 3 /J.K 2 at i — 2000), while 
the Planck, ACT, and SPT data have excess power over the best- 
fit Al = 1 ACDM+foreground model at the level of a few yuK 2 
(see Sect. 4). It is plausible that this may drive the preference 
for high A L in the x 1 of the high-/" experiments. We note that a 
similar 2 <x preference for Al > 1 is also found combining ACT 
and WMAP data (Sievers et al. 2013) and, as we find here, this 
tension is reduced when the lensing power spectrum is included 
in the fit. 

To summarize, there is no preference in the Planck lensing 
power spectrum for Al > 1. The general preference for high 
Al from the CMB power spectra in our favoured data combi- 
nation (P/anc£+ WP+highL) is mostly driven by two effects: the 
difficulty that ACDM models have in fitting the low-£ spectrum 
when calibrated from the smaller-scale spectrum; and, plausibly, 
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from excess residuals at the juK 2 level in the high-^ spectra rela- 
tive to the best-fit Al = 1 ACDM+foregrounds model on scales 
where extragalactic foreground modelling is critical. 

5.2. Baryon acoustic oscillations 

Baryon acoustic oscillations (BAO) in the matter power spec- 
trum were first detected in analyses of the 2dF Galaxy 
Redshift Survey (Cole et al. 2005) and the SDSS redshift sur- 
vey (Eisenstein et al. 2005). Since then, accurate BAO measure- 
ments have been made using a number of different galaxy red- 
shift surveys, providing constraints on the distance luminosity 
relation spanning the redshift range 0.1 ^ z ^ 0.7 18 . Here we use 
the results from four redshift surveys: the SDSS DR7 BAO mea- 
surements at effective redshifts z e ff = 0.2 and z e ff = 0.35, anal- 
ysed by Percival et al. (2010); the z = 0.35 SDSS DR7 measure- 
ment at z e ff = 0.35 reanalyzed by Padmanabhan et al. (2012); the 
WiggleZ measurements at z e s = 0.44, 0.60 and 0.73 analysed by 
Blake et al. (2011); the BOSS DR9 measurement at z eff = 0.57 
analyzed by Anderson et al. (2013); and the 6dF Galaxy Survey 
measurement at z = 0.1 discussed by Beutler et al. (201 1). 
BAO surveys measure the distance ratio 

r sfoirag) 

d * = n~TT ' (45) 
Dv(z) 

where r s (zd iag ) is the comoving sound horizon at the baryon drag 
epoch (when baryons became dynamically decoupled from the 
photons) and Dy(z) is a combination of the angular-diameter dis- 
tance, Da(z), and the Hubble parameter, H{z), appropriate for the 
analysis of spherically-averaged two-point statistics: 



(l+zfD 2 A (z) 



H(z) 



1/3 



(46) 



In the ACDM cosmology, the angular diameter distance to red- 
shift z is 

D A (z) = ^rD A . 
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c 



1 



H \n K \u 2 (i + z) 



sin^f|Q^| 1/2 x(z,Q m ,Q A )], (47) 



where 



x(z, O m , C2 A ) 
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dzf 



(48) 



[fi m (l+z') 3 +^(l+z') 2 + "A] 1/2 ' 
and sin# = sinh for O.^ > and sin^ — sin for Q.k < 0. Note 
that the luminosity distance, Z\, relevant for the analysis of Type 
la supernovae (see Sect. 5.4) is related to the angular diameter 
distance via D h - (c/H )D h - D A (l + z) 2 . 

Different groups fit and characterize BAO features in differ- 
ent ways. For example, the WiggleZ team encode some shape 
information on the power spectrum to measure the acoustic pa- 
rameter A(z), introduced by Eisenstein et al. (2005), 



A(z) = 



Dyfe)^// 2 



cz 



(49) 



18 Detections of a BAO feature have recently been reported in the 
three-dimensional correlation function of the Lya forest in large sam- 
ples of quasars at a mean redshift of z ~ 2.3 (Buscaetal. 2012; 
Slosar et al. 2013). These remarkable results, probing cosmology well 
into the matter-dominated regime, are based on new techniques that are 
less mature than galaxy BAO measurements. For this reason, we do not 
include Lya BAO measurements as supplementary data to Planck. For 
the models considered here and in Sect. 6, the galaxy BAO results give 
significantly tighter constraints than the Lyo- results. 
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Fig. 15. Acoustic-scale distance ratio r s /D v (z) divided by the 
distance ratio of the Planck base ACDM model. The points are 
colour-coded as follows: green star (6dF); purple squares (SDSS 
DR7 as analyzed by Percival et al. 2010); black star (SDSS DR7 
as analyzed by Padmanabhan et al. 2012); blue cross (BOSS 
DR9); and blue circles (WiggleZ). The grey band shows the ap- 
proximate ± 1 <x range allowed by Planck (computed from the 
CosmoMC chains). 



which is almost independent of u> m . To simplify the presenta- 
tion, Fig. 15 shows estimates of r s /Dy(z) and 1 <x errors, as 
quoted by each of the experimental groups, divided by the ex- 
pected relation for the Planck base ACDM parameters. Note 
that the experimental groups use the approximate formulae of 
Eisenstein & Hu (1998) to compute Zd rag an d r s (zdragX though 
they fit power spectra computed with Boltzmann codes, such 
as camb, generated for a set of fiducial-model parameters. The 
measurements have now become so precise that the small differ- 
ence between the Eisenstein & Hu (1998) approximations and 
the accurate values of Zdrag and rd iag = r s (zdrag) returned by camb 
need to be taken into account. In CosmoMC we multiply the ac- 
curate numerical value of r s (zdrag) by a constant factor of 1.0275 
to match the Eisenstein-Hu approximation in the fiducial model. 
This correction is sufficiently accurate over the range of a> m and 
Wb allowed by the CMB in the base ACDM cosmology (see e.g. 
Mehta et al. 2012) and also for the extended ACDM models dis- 
cussed in Sect. 6. 

The Padmanabhan et al. (2012) result plotted in Fig. 15 is 
a reanalysis of the z e ff = 0.35 SDSS DR7 sample discussed 
by Percival et al. (2010). Padmanabhan et al. (2012) achieve a 
higher precision than Percival et al. (2010) by employing a re- 
construction technique (Eisenstein et al. 2007) to correct (par- 
tially) the baryon oscillations for the smearing caused by galaxy 
peculiar velocities. The Padmanabhan et al. (2012) results are 
therefore strongly correlated with those of Percival et al. (2010). 
We refer to the Padmanabhan et al. (2012) "reconstruction- 
corrected" results as SDSS(R). A similar reconstruction tech- 
nique was applied to the BOSS survey by Anderson et al. (2013) 
to achieve 1.6% precision in Dy(z = 0.57)/r s , the most precise 
determination of the acoustic oscillation scale to date. 

All of the BAO measurements are compatible with the base 
ACDM parameters from Planck. The grey band in Fig. 15 
shows the +1 cr range in the acoustic-scale distance ratio com- 
puted from the Planck+ WP+highL CosmoMC chains for the base 
ACDM model. To get a qualitative feel for how the BAO mea- 
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Table 8. Approximate constraints with 68% errors on £2 m and 
Ho (in units of km s _1 Mpc~') from BAO, with a> m and Wb fixed 
to the best-fit Planck+WP+highL values for the base ACDM 
cosmology. 



Sample fi m H a 



6dF 


305 +0032 

U ' JUJ -0.026 


68.3 +1 ~ 


SDSS 


A 9QC+0.019 

U " S ' J -0.M7 


69.5!|? 


SDSS(R) 


293 +0015 

"""■'-0.013 


69-6!!;? 


WiggleZ 


309 +0041 

"■■' -0.035 


67.8!^ 


BOSS 


315+ 0015 


67.2+^ 


6dF+SDSS+BOSS+WiggleZ 


307+ 010 




6dF+SDSS(R)+BOSS 


305 + 009 


68.4+j;° 


6dF+SDSS(R)+BOSS+WiggleZ 


305 +o m 


68.4!!» 



surements constrain parameters in the base ACDM model, we 

form^- 2 , 

X 2 BA0 = (x - x^/C-^ix ~ * ACDM ), (50) 

where x is the data vector, jc acdm denotes the theoretical pre- 
diction for the ACDM model and C B l A0 is the inverse covari- 
ance matrix for the data vector x. The data vector is as fol- 
lows: D v (0.106) = (457 + 27)Mpc (6dF); r s /D v (0.20) = 
0.1905 ± 0.0061, r s /ZJ v (0.35) = 0.1097 ± 0.0036 (SDSS); 
A(0.44) = 0.474 ± 0.034, A(0.60) = 0.442 ± 0.020, A(0.73) = 
0.424 + 0.021 (WiggleZ); £> v (0.35)/r s = 8.88 + 0.17 (SDSS(R)); 
andD v (0.57)/r s = 13.67 + 0.22, (BOSS). The off-diagonal com- 
ponents of C B ^ for the SDSS and WiggleZ results are given 
in Percival et al. (2010) and Blake et al. (2011). We ignore any 
covariances between surveys. Since the SDSS and SDSS(R) re- 
sults are based on the same survey, we include either one set of 
results or the other in the analysis described below, but not both 
together. 

The Eisenstein-Hu values of r s for the Planck and WMAP-9 
base ACDM parameters differ by only 0.9%, significantly 
smaller than the errors in the BAO measurements. We can obtain 
an approximate idea of the complementary information provided 
by BAO measurements by minimizing Eq. (50) with respect to 
either Q m or Hq, fixing a> m and a>b to the CMB best-fit parame- 
ters. (We use the Planck+ WP+highL parameters from Table 5.) 
The results are listed in Table 8 19 . 

As can be seen, the results are very stable from survey to 
survey and are in excellent agreement with the base ACDM 
parameters listed in Tables 2 and 5. The values of Xbao ^ 
also reasonable. For example, for the six data points of the 
6dF+SDSS(R)+BOSS+WiggleZ combination, we find xl A0 = 
4.3, evaluated for the Planck+ WP+highL best-fit ACDM param- 
eters. 

The high value of Q m is consistent with the parameter anal- 
ysis described by Blake et al. (201 1) and with the "tension" dis- 
cussed by Anderson et al. (2013) between BAO distance mea- 
surements and direct determinations of Hq (Riess et al. 2011; 
Freedman et al. 2012). Furthermore, if the errors on the BAO 
measurements are accurate, the constraints on D. m and Ho (for 
fixed aj m and cj^,) are of comparable accuracy to those from 
Planck. 

19 As an indication of the accuracy of Table 8, the full likelihood 
results for the P/a«ofc+WP+6dF+SDSS(R)+BOSS BAO data sets give 
CK = 0.308 ± 0.010 and H = 67.8 ± 0.8 km s 1 Mpc~\ for the base 
ACDM model. 
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Fig. 16. Comparison of Ho measurements, with estimates of 
±1 <x errors, from a number of techniques (see text for details). 
These are compared with the spatially-flat ACDM model con- 
straints from Planck and WMAP-9. 

The results of this section show that BAO measurements are 
an extremely valuable complementary data set to Planck. The 
measurements are basically geometrical and free from complex 
systematic effects that plague many other types of astrophysical 
measurements. The results are consistent from survey to survey 
and are of comparable precision to Planck. In addition, BAO 
measurements can be used to break parameter degeneracies that 
limit analyses based purely on CMB data. For example, from 
the excellent agreement with the base ACDM model evident in 
Fig. 15, we can infer that the combination of Planck and BAO 
measurements will lead to tight constraints favouring £1% = 
(Sect. 6.2) and a dark energy equation-of-state parameter, w = 
-1 (Sect. 6.5). 

Finally, we note that we choose to use the 
6dF+SDSS(R)+BOSS data combination in the likelihood 
analysis of Sect. 6. This choice includes the two most accu- 
rate BAO measurements and, since the effective redshifts of 
these samples are widely separated, it should be a very good 
approximation to neglect correlations between the surveys. 

5.3. The Hubble constant 

A striking result from the fits of the base ACDM model to Planck 
power spectra is the low value of the Hubble constant, which is 
tightly constrained by CMB data alone in this model. From the 
Planck+ WP+highL analysis we find 

H = (67.3+1.2) km s" 1 Mpc" 1 (68%; Planck+ WP+highL).(51) 

A low value of Ho has been found in other CMB experi- 
ments, most notably from the recent WMAP-9 analysis. Fitting 
the base ACDM model, Hinshaw et al. (2012) find 

H = (70.0 + 2.2)kms" 1 Mpc^ 1 (68%; WMAP-9), (52) 

consistent with Eq. (51) to within 1 <x. We emphasize here that 
the CMB estimates are highly model dependent. It is important 
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Fig. 17. MCMC samples and contours in the r t -D. m h 2 plane (left) and the D A (z*)-Cl m h 2 plane (right) for ACDM models analysed 
with P/a«cfe+WP+highL. The lines in these plots show the expected degeneracy directions in the base ACDM cosmology. Samples 
are colour-coded by the values of Ob^ 2 (left) and Hq (right). 



therefore to compare with astrophysical measurements of Hq, 
since any discrepancies could be a pointer to new physics. 

There have been remarkable improvements in the preci- 
sion of the cosmic distance scale in the last decade or so. 
The final results of the Hubble Space Telescope (HST) Key 
Project (Freedman et al. 2001), which used Cepheid calibrations 
of secondary distance indicators, resulted in a Hubble constant 
of H = (72 ± 8)kms"' Mpc" 1 (where the error includes esti- 
mates of both 1 cr random and systematic errors). This estimate 
has been used widely in combination with CMB observations 
and other cosmological datasets to constrain cosmological pa- 
rameters (e.g., Spergel et al. 2003, 2007). It has also been recog- 
nized that an accurate measurement of Ho with around 1 % pre- 
cision, when combined with CMB and other cosmological data, 
has the potential to reveal exotic new physics, for example, a 
time-varying dark energy equation of state, additional relativistic 
particles, or neutrino masses (see e.g., Suyu et al. 2012b, and ref- 
erences therein). Establishing a more accurate cosmic distance 
scale is, of course, an important problem in its own right. The 
possibility of uncovering new fundamental physics provides an 
additional incentive. 

Two recent analyses have greatly improved the precision of 
the cosmic distance scale. Riess et al. (2011) use HST observa- 
tions of Cepheid variables in the host galaxies of eight SNe la to 
calibrate the supernova magnitude-redshift relation. Their "best 
estimate" of the Hubble constant, from fitting the calibrated SNe 
magnitude-redshift relation, is 

Ho = (73.8 ± 2.4) km s" 1 Mpc" 1 (Cepheids+SNe la), (53) 

where the error is 1 cr and includes known sources of systematic 
errors. At face value, this measurement is discrepant with the 
Planck estimate in Eq. (51) at about the 2.5 cr level. 

Freedman et al. (2012), as part of the Carnegie Hubble 
Program, use Spitzer Space Telescope mid-infrared observations 
to recalibrate secondary distance methods used in the HST Key 
Project. These authors find 

H = [74.3 ±1.5 (statistical) ± 2.1 (systematic)] km s -1 Mpc -1 

(Carnegie HP). (54) 

We have added the two sources of error in quadrature in the 
error range shown in Fig. 16. This estimate agrees well with 



Eq. (53) and is also discordant with the Planck value (Eq. 16) 
at about the 2.5 cr level. The error analysis in Eq. (54) does not 
include a number of known sources of systematic error and is 
very likely an underestimate. For this reason, and because of the 
relatively good agreement between Eqs. (53) and (54), we do not 
use the estimate in Eq. (54) in the likelihood analyses described 
in Sect. 6. 

The dominant source of error in the estimate in Eq. (53) 
comes from the first rung in the distance ladder. Using the 
megamaser-based distance to NGC4258, Riess et al. (201 1) find 
(74.8 + 3.1)kms _1 Mpc~'. Using parallax measurements for 10 
Milky Way Cepheids, they find (75.7 ± 2.6)kms~ 1 Mpc~', and 
using Cepheid observations and a revised distance to the Large 
Magellanic Cloud, they find (71.3 ± 3.8) km s^ 1 Mpc -1 . These 
estimates are consistent with each other, and the combined esti- 
mate in Eq. (53) uses all three calibrations. The fact that the er- 
ror budget of measurement (53) is dominated by the "first-rung" 
calibrators is a point of concern. A mild underestimate of the 
distance errors to these calibrators could eliminate the tension 
with Planck. 

Figure 16 includes three estimates of Hq based on "geomet- 
rical" methods 20 . The estimate labelled "UGC 3789" shows the 
result Hq = (71.6 + 5.7)kms _1 Mpc~' from VLBI observations 
of water masers orbiting within the accretion disc of UCG 3789 
(Reidetal. 2012). The point labelled "RXJ1 131-1231" shows 
the estimate Hq = 78.7^ ? kms -1 Mpc~' derived from gravita- 
tional lensing time delay measurements of the system RXJ1 131- 
1231, observed as part of the "COSmological MOnitoring of 
GRAvitational Lenses" (COSMOGRAIL) project (Suyu et al. 
2012a). Finally, the point labelled SZ clusters shows the value 
H = 76.9^° 7 7 kms- 1 Mpc -1 (Bonamente et al. 2006), derived 
by combining tSZ and X-ray measurements of rich clusters 
of galaxies (see Carlstrom et al. 2002, and references therein). 
These geometrical methods bypass the need for local distance 
calibrators, but each has its own sources of systematic error 
that need to be controlled. The geometrical methods are con- 
sistent with the Cepheid-based methods, but at present, the er- 
rors on these methods are quite large. The COSMOGRAIL mea- 
surement (which involved a "blind" analysis to prevent experi- 



20 Note that each of these estimates is weakly dependent on the as- 
sumed background cosmology. 
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menter bias) is discrepant at about 2.5 <x with the Planck value 
in Eq. (51). We note here a number of other direct measurements 
of H (Sandage et al. 2006; Oguri 2007; Jones et al. 2005) that 
give lower values than the measurements summarized in Fig. 16. 

The tension between the CMB-based estimates and the astro- 
physical measurements of H () is intriguing and merits further dis- 
cussion. In the base ACDM model, the sound horizon depends 
primarily on Cl m h 2 (with a weaker dependence on Q.hh 2 ). This is 
illustrated by the left-hand panel of Fig. 17, which shows sam- 
ples from the P/ancfc+WP+highL MCMC chains in the r„-Q m /i 2 
plane colour coded according to Qt,/i 2 . The acoustic scale param- 
eter 6* is tightly constrained by the CMB power spectrum, and so 
a change in r, must be matched by a corresponding shift in the 
angular diameter distance to the last scattering surface Da(z*)- 
In the base ACDM model, Da depends on Ho and Q m /i 2 , as 
shown in the right-hand panel of Fig. 17. The 2.7 km s" 1 Mpc~' 
shift in H) between Planck and WMAP-9 is primarily a con- 
sequence of the slightly higher matter density determined by 
Planck (Q m /z 2 = 0.143 + 0.003) compared to WMAP-9 (Q. m h 2 = 
0.136 + 0.004). A shift of around 7kms~' Mpc~', necessary to 
match the astrophysical measurements of Ho would require an 
even larger change in Cl m h 2 , which is disfavoured by the Planck 
data. The tension between Planck and the direct measurements 
of Ho cannot be easily resolved by varying the parameters of 
the base ACDM model. Section 6 will explore whether there are 
any extensions to the base ACDM model that can relieve this 
tension. In that section, results labelled "Ho" include a Gaussian 
prior on Hq based on the Riess et al. (201 1) measurement given 
in Eq. (53). 



5.4. Type la supernovae 

In this subsection, we analyse two SNe la samples: the sample of 
473 SNe as reprocessed by Conley et al. (2011), which we will 
refer to as the "SNLS" compilation; and the updated Union2.1 
compilation of 580 SNe described by Suzuki et al. (2012). 



5.4.1. The SNLS compilation 

The SNLS "combined" compilation consists of 123 SNe la 
at low redshifts, 242 SNe la from the three-year Supernova 
Legacy Survey (SNLS; see Regnault et al. 2009; Guy et al. 
2010; Conley et al. 201 1), 93 intermediate redshift SNe la from 
the Sloan Digital Sky Survey (SDSS; Holtzman et al. 2008; 
Kessler et al. 2009) and 14 objects at high redshift observed with 
the Hubble Space Telescope (HST; Riess et al. 2007). 

The "combined" sample of Conley et al. (2011) combines 
the results of two light-curve fitting codes, SiFTO (Conley et al. 
2008) and SALT2 (Guy et al. 2007), to produce a peak appar- 
ent B-band magnitude, nts, stretch parameter s and colour C for 
each supernova. To explore the impact of light-curve fitting, we 
also analyse separately the SiFTO and SALT2 parameters. The 
SiFTO and SALT2 samples differ by a few SNe from the com- 
bined sample because of colour and stretch constraints imposed 
on the samples. We also use ancillary data, such as estimates of 
the stellar masses of the host galaxies and associated covariance 
matrices, as reported by Conley et al. (201 1) 21 . 



21 https : //tspace . library . utoronto . ca/handle/1807/ 
25398. We use the CosmoMC module downloadable from this site. 



In this section, we focus exclusively on the base ACDM 
model (i.e., w = — 1 and Q.k = 0). For a flat Universe, the ex- 
pected apparent magnitudes are then given by 

mf DM = 51og 10 D L (z hel ,zcMB,O m ) - a(s - 1) +(3C + M B , (55) 

where is the dimensionless luminosity distance 22 and Mb 
absorbs the Hubble constant. As in Sullivan et al. (2011), we ex- 
press values of the parameter(s) Mb in terms of an effective ab- 
solute magnitude 

M B = M B -51og 10 ^J-25, (56) 

for a value of Ho = 70 km s _1 Mpc~' . 

The likelihood for this sample is then constructed as in 
Conley et al. (2011) and Sullivan et al. (2011): 

* 2 Ne = (M B - Mf™) T C s l(M B - Mf DM ), (57) 

where Mb is the vector of effective absolute magnitudes and 
CsNe is the sum of the non-sparse covariance matrices of 
Conley et al. (2011) quantifying statistical and systematic er- 
rors. As in Sullivan et al. (2011), we divide the sample accord- 
ing to the estimated stellar mass of the host galaxy and solve for 
two parameters, M l B for M host < 10 10 M G and M 2 for M host > 
10 ll) M o . We adopt the estimates of the "intrinsic" scatter in m B 
for each SNe sample given in Table 4 of Conley et al. (201 1). 

Fits to the SNLS combined sample are shown in the left- 
hand panel of Fig. 18. The best-fit parameters for the combined, 
SiFTO and SALT2 samples are given in Table 9. In the base 
ACDM model, the SNe data provide a constraint on O m , inde- 
pendent of the CMB. As can be seen from Table 9 (and also 
in the analyses of Conley et al. 2011 and Sullivan et al. 2011), 
the SNLS combined compilation favours a lower value of Q. m 
than we find from the CMB. The key question, of course, is 
whether the SNe data are statistically compatible with the Planck 
data. The last three rows of Table 9 give the best-fit SNe param- 
eters constraining Sl m to the Planck+ WP+highL best-fit value 
£2 m = 0.317. The grey bands in Fig. 18 show the magnitude 
residuals expected for a ±2 cr variation in the value of Q m al- 
lowed by the CMB data. The CMB band lies systematically low 
by about 0.1 magnitude over most of the redshift range shown in 
Fig. 18a. 

Table 9 also lists the^ 2 values for the Q m = 0.317 fits 23 . The 
likelihood ratio for the SiFTO fits is 

= exp( ifck -*SNe + CMB oj) * 0.074. (58) 

This is almost a 2 cr discrepancy. (The discrepancy would appear 
to be much more significant if only the diagonal statistical errors 
were included in the covariance matrix in Eq. 57). The likeli- 
hood ratio for the combined sample is slightly larger (0.095) and 
is larger still for the SALT2 sample (0.33). In summary, there 
is some tension between the SNLS compilations and the base 
ACDM value of Q. m derived from Planck. The degree of tension 
depends on the light-curve fitter and is stronger for the SiFTO 
and combined SNLS compilations. 

22 Note that the luminosity distance depends on both the heliocentric, 
Zheb and CMB frame, zcmb, redshifts of the SNe. This distinction is 
important for low-redshift objects. 

23 We caution the reader that, generally, the ;^ Nc obtained from 
Eq. (57) will differ from that quoted in the online parameter tables in 
cases where the SNLS data is importance sampled. For importance sam- 
pling, we modified the SNLS likelihood to marginalize analytically over 
the a and fi parameters. 
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Fig. 18. Magnitude residuals relative to the base ACDM model that best fits the SNLS combined sample (left) and the Union2.1 
sample (right). The error bars show the 1 <x (diagonal) errors on m B . The filled grey regions show the residuals between the expected 
magnitudes and the best-fit to the SNe sample as Q m varies across the ±2 <x range allowed by Planck+WP +highL in the base 
ACDM cosmology. The colour coding of the SNLS samples are as follows: low redshift (blue points); SDSS (green points); SNLS 
three-year sample (orange points); and HST high redshift (red points). 



Table 9. Best-fit parameters for the SNLS compilations. 



Data set 


NsNe 


K 


K 




P 




x 1 


SNLS combined 


472 


-19.16 


-19.21 


1.425 


3.256 


0.227 


407.8 


SNLS SiFTO 


468 


-19.15 


-19.20 


1.352 


3.375 


0.223 


414.9 


SNLS SALT2 


473 


-19.15 


-19.20 


1.698 


3.289 


0.247 


376.7 


SNLS combined (CMB Q m ) 


472 


-19.12 


-19.18 


1.417 


3.244 


0.317 


412.5 


SNLS SiFTO (CMB fi m ) 


468 


-19.12 


-19.18 


1.339 


3.351 


0.317 


420.1 


SNLS SALT2 (CMB £2 m ) 


473 


-19.12 


-19.18 


1.691 


3.302 


0.317 


378.9 



5.4.2. The Union2.1 compilation 

The Union2.1 compilation (Suzuki et al. 2012) is the latest ap- 
plication of a scheme for combining multiple SNe data sets 
described by Kowalski et al. (2008). The Union2.1 compilation 
contains 19 data sets and includes early high-redshift SNe data 
(e.g., Riess et al. 1998; Perlmutter et al. 1999) as well as recent 
data from the HST Cluster Supernova Survey (Amanullah et al. 
2010; Suzuki et al. 2012). The SNLS and Union2.1 compilations 
contain 256 SNe in common and are therefore not independent. 

The SALT2 model (Guy et al. 2007) is used to fit the light 
curves returning a B-band magnitude at maximum light, a light- 
curve shape parameter and a colour correction. As in Eq. (55), 
the theoretically-predicted magnitudes include nuisance param- 
eters a and (3 multiplying the shape and colour corrections, and 
an additional nuisance parameter 5 describing the variation of 
SNe luminosity with host galaxy mass (see Eq. 3 of Suzuki et al. 
2012). The CosmoMC module associated with the Union2.1 sam- 
ple 24 holds the nuisance parameters fixed (a = 0.1218, /3 = 
2.4657, and 6 = -0.03634) and computes a x 2 via Eq. (57) us- 
ing a fixed covariance matrix that includes a model for system- 
atic errors. An analysis of the base ACDM model then requires 
minimization with respect to only two parameters, O m and Mb 
(or equivalently, Mb). It is a slight concern that the Union2.1 
CosmoMC module does not allow the nuisance parameters to vary 
along with the cosmological parameters. Conley et al. (2011) 
show that fixing the SNLS a and /3 parameters introduces small 
biases in cosmological parameters (see Fig. 1 1 of Conley et al. 
201 1). However, these biases are small enough that they are un- 
likely to affect any of our conclusions in any significant way. 

24 http : // supernova . lbl . gov/Union. 



Maximizing the Union2.1 likelihood, we find best-fit param- 
eters of Q m = 0.296 and M B = -19.272 (defined as in Eq. 56 for 
a value of H = 70kms- 1 Mpc" 1 ) and* 2 nion21 = 545.11 (580 
SNe). The magnitude residuals with respect to this fit are shown 
in the right-hand panel of Fig. 18. Notice that the scatter in this 
plot is significantly larger than the scatter of the SNLS compi- 
lation (left-hand panel) reflecting the more diverse range of data 
and the lower precision of some of the earlier SNe data used in 
the Union2.1 compilation. Nevertheless, the Union2.1 best-fit is 
close to (and clearly compatible with) the Planck base ACDM 
value of Q. m . 

5.4.3. SNe: Summary 

The results of this subsection are summarized in Fig. 19. This 
shows the posterior distributions for Q m in the base ACDM cos- 
mology, marginalized over nuisance parameters, for each of the 
SNe samples. These distributions are broad (with the Union2.1 
distribution somewhat broader than the SNLS distributions) and 
show substantial overlap. There is no obvious inconsistency be- 
tween the SNe samples. The posterior distribution for Q m in the 
base ACDM model fit to Planck+WP +highL is shown by the 
narrow green curve. This is consistent with the Union2.1 and 
SNLS SALT2 results, but is in some tension with the distribu- 
tions from the SNLS combined and SNLS SiFTO samples. As 
we will see in Sect. 6, Planck combined with Planck lensing and 
BAO measurements overwhelm SNe data for most of the exten- 
sions of the ACDM model considered in this paper. However, the 
results presented here suggest that there could be residual sys- 
tematic errors in the SNe data that are not properly accounted 
for in the covariance matrices. Hints of new physics based on 
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Fig. 19. Posterior distributions for Q m (assuming a flat cosmol- 
ogy) for the SNe compilations described in the text. The poste- 
rior distribution for Q. m from the PZa«cA:+WP+highL fits to the 
base ACDM model is shown by the solid green line. 

combining CMB and SNe data should therefore be treated with 
caution. 

5.5. Additional data 

In this subsection we review a number of other astrophysical data 
sets that have sometimes been combined with CMB data. These 
data sets are not used with Planck in this paper, either because 
they are statistically less powerful than the data reviewed in pre- 
vious subsections and/or they involve complex physics (such as 
the intra-cluster gas in rich clusters of galaxies) which is not yet 
well understood. 

5.5.1. Shape information on the galaxy/matter power 
spectrum 

Reid et al. (2010) present an estimate of the dark matter 
halo power spectrum, Phaio(^), derived from 110,756 lumi- 
nous red galaxies (LRGs) from the SDSS 7th data release 
(Abazajian et al. 2009). The sample extends to redshifts z ~ 0.5, 
and is processed to identify LRGs occupying the same dark 
matter halo, reducing the impact of redshift-space distortions 
and recovering an approximation to the halo density field. The 
power spectrum Phaio(^) is reported in 45 bands, covering the 
wavenumber range 0.02 /iMpc" 1 < k < 0.2/iMpc-'. The win- 
dow functions, covariance matrix and CosmoMC likelihood mod- 
ule are available on the NASA LAMBDA web site 25 . 

The halo power spectrum is plotted in Fig. 20. The blue line 
shows the predicted halo power spectrum from our best-fit base 
ACDM parameters convolved with the Reid et al. (2010) win- 
dow functions. Here we show the predicted halo power spec- 
trum for the best-fit values of the "nuisance" parameters bo 
(halo bias), a\, and «2 (defined in equation 15 of Reid et al. 
2010) which relate the halo power spectrum to the dark mat- 
ter power spectrum (computed using camb). The Planck model 
gives X\jrg ~ 46.9 for 42 degrees of freedom and is an ac- 
ceptable, though marginally worse, fit than the best-fit model 



http : //lambda . gsf c . nasa . gov/ toolbox/lrgdr. 
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Fig. 20. Band-power estimates of the halo power spectrum, 
Phaio(k), from Reid et al. (2010) together with 1 <x errors. (Note 
that these data points are strongly correlated.) The line shows 
the predicted spectrum for the best-fit Planck+WP+highL base 
ACDM parameters. 

of Reid et al. (2010), which has xt RG = 40.0. Interestingly, the 
main differences between the two models are at wavenumbers 
k ^ 0.1 /zMpc~', where the nonlinear corrections to the matter 
power spectrum become important. 

Figure 20 shows that the Planck parameters provide a good 
match to the shape of the halo power spectrum. However, we do 
not use these data (in this form) in conjunction with Planck. The 
BAO scale derived from these and other data is used with Planck, 
as summarized in Sect. 5.2. As discussed by Reid et al. (2010, 
see their figure 5) there is very little additional information on 
cosmology once the BAO features are filtered from the spec- 
trum, and hence little to be gained by adding this information to 
Planck. The corrections for nonlinear evolution, though small in 
the wavenumber range 0.1-0.2 ZiMpc -1 , add to the complexity 
of using shape information from the halo power spectrum. 

5.5.2. Cosmic shear 

Another key cosmological observable is the distortion of distant 
galaxy images by the gravitational lensing of large-scale struc- 
ture, often called cosmic shear. The shear probes the (nonlinear) 
matter density projected along the line of sight with a broad ker- 
nel. It is thus sensitive to the geometry of the Universe and the 
growth of large-scale structure, with a strong sensitivity to the 
amplitude of the matter power spectrum. 

The most recent, and largest, cosmic shear data sets are 
provided by the CFHTLenS survey (Hey mans et al. 2012; 
Erben et al. 2012), which covers 26 154deg 2 in five optical 
bands with accurate shear measurements and photometric 
redshifts. The CFHTLenS team has released several cosmic 
shear results which are relevant to this paper. Benjamin et al. 
(2012) present results from a two-bin tomographic analysis, 
Heymans et al. (2013) from a finely binned tomographic anal- 
ysis, and Kitching et al. (2013) from a 3D analysis. 

Heymans et al. (2013) estimate shear correlation func- 
tions associated with six redshift bins. Assuming a flat, 
ACDM model, from the weak lensing data alone they find 
cr 8 (fi m /0.27)°- 46±0 - 02 = 0.774 + 0.04 (68% errors) which is con- 



' Approximately 61% of the survey is fit for cosmic shear science. 
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sistent with the constraint found by Benjamin et al. 2012. For 
comparison, we find 

o- g (Q m /0.27) a46 = 0.89 ± 0.03 (68%; Planck+ WP+highL), 

(59) 

which is discrepant at the 2<x level. Combining the tomo- 
graphic lensing data with CMB constraints from WMAP-1, 
Heymans et al. (2013) are able to constrain the individual pa- 
rameters of the flat, ACDM model to be Q m = 0.255 ±0.014 
and h = 0.717 ± 0.016. The best-fit Planck value of Q m is 4 cr 
away from this value, while h is nearly 3 cr discrepant. As might 
be expected, given the good agreement between the Planck and 
BAO distance scales, the best-fit CFHTLenS ACDM cosmology 
is also substantially discrepant with the BOSS data, predicting 
a distance ratio to z — 0.57 which is 5% lower than measured 
by BOSS (Anderson et al. 2013). This is discrepant at approx- 
imately the 3 cr level, comparable to the discrepancy with the 
Planck values. 

Kitching et al. (2013) have presented constraints on cosmo- 
logical models from a 3D analysis of cosmic shear, in com- 
bination with WMAP data. The 3D analysis limits the radial 
and angular modes used to k < 1.5/iMpc -1 and 400 < 
I < 4000, respectively, which is more conservative than the 
cuts used by Benjamin et al. (2012) and Heymans et al. (2013). 
Kitching et al. (2013) do not present results for our fiducial six- 
parameter ACDM model and the lensing data alone do not pro- 
vide strong constraints on the cosmologies they consider. When 
combined with CMB constraints from WMAP-1, for a flat CDM 
model with a non-evolving dark energy equation of state they 
findQ m = 0.293 ± 0.033, h = 0.715 ±0.038, cr 8 = 0.872 ±0.057, 
and w = -0.968 ± 0.12. For the same model Heymans et al. 
(2013) find Q m = 0.256+°;^', h = 0.74 ±0.13, cr 8 = 0.81 ±0.10, 
and w = -1.05 ± 0.33, which are consistent given the large er- 
rors in the latter case. By comparison the best-fit parameters for 
such a model from P/anc£+ WP+highL+BAO are Q m = 0.29, 
h = 0.70, cr 8 = 0.87, and w - —1.11 in excellent agreement with 
the 3D shear+WMAP-7 values. 

The source of the discrepancies between Planck and the 
CFHTLenS tomographic analyses is at present unclear, and fur- 
ther work will be needed to resolve them. 

5.5.3. Counts of rich clusters 

For the base ACDM model we find cr 8 = 0.828 ± 0.012 from 
Planck+ WP+highL. This value is in excellent agreement with 
the WMAP-9 value of cr 8 = 0.821 ± 0.023 (Hinshaw et al. 
2012). There are other ways to probe the power spectrum nor- 
malization, in addition to the cosmic shear measurements dis- 
cussed above. For example, the abundances of rich clusters of 
galaxies are particularly sensitive to the normalization (see e.g., 
Komatsu & Seljak 2002). Recently, a number of studies have 
used tSZ-cluster mass scaling relations to constrain combina- 
tions of cr 8 and Q m (e.g., Benson et al. 2013; Reichardt et al. 
2013; Hasselfield et al. 2013) including an analysis of a sam- 
ple of Planck tSZ clusters reported in this series of papers 
(Planck Collaboration XX 2013) 27 . 

The Planck analysis uses a relation between cluster mass 
and tSZ signal based on comparisons with X-ray mass measure- 
ments computed assuming hydrostatic equilibrium. Because the 



27 There is additionally a study of the statistical properties of the 
P/anc^-derived Compton-y map (Planck Collaboration XXI 2013) from 
which other parameter estimates can be obtained. 



intra-cluster gas is not expected to be in perfect hydrostatic equi- 
librium, Planck Collaboration XX (2013) assume a bias of 20% 
between the X-ray derived masses and the true cluster masses. 
For this value they find cr 8 (Q m /0.27) 03 = 0.79 ± 0.01 for the 
base ACDM model. In comparison, for the same model we find 

cr 8 (Q m /0.27) 3 = 0.87 ± 0.02 (68%; Planck+ WP+highL), 

(60) 

which is a significant (3 cr) discrepancy that remains unex- 
plained. Qualitatively similar results, though at lower precision, 
are found from analyses of SPT clusters [cr 8 (Q m /0.27)°- 3 
0.77 ± 0.04]. Key difficulties with this type of measurement, 
as discussed in Planck Collaboration XX (2013), include ad- 
equately modelling selection biases and calibrating cluster 
masses. These effects are discussed in the analysis of ACT clus- 
ters by Hasselfield et al. (2013), who adopt a number of ap- 
proaches, including folding in dynamical mass measurements, 
to calibrate biases in clusters mass estimates. Some of these ap- 
proaches give joint cr 8 -Q m constraints, consistent with the base 
ACDM parameters reported here. 

At this stage of our understanding of the biases and scat- 
ter in the cluster mass calibrations, we believe that for the pur- 
poses of this paper it is premature to use cluster counts to- 
gether with CMB measurements to search for new physics. 
Planck Collaboration XX (2013) explore a number of possi- 
bilities for reducing the tension between Planck CMB mea- 
surements and tSZ cluster counts, including non-zero neutrino 
masses. 



6. Extensions to the base ACDM model 
6. 1 . Grid of models 

To explore possible deviations from ACDM we have analysed 
an extensive grid of models that covers many well-motivated ex- 
tensions of ACDM. As in the exploration of the base ACDM 
cosmology, we have also considered a variety of data combina- 
tions for each model. For models involving more than one ad- 
ditional parameter we restrict ourselves to Planck+WP combi- 
nations in order to obtain tighter constraints by leveraging the 
relative amplitude of the power spectrum at very low I and 
high t. Most models are run with Planck, Planck+WP, and 
Planck+ WP+highL; additionally all are importance sampled 
with Planck lensing (Sect. 5.1), BAO (Sect. 5.2), SNe (Sect. 5.4), 
and the Riess et al. (2011) direct Hq measurement (Sect. 5.3). 
For models where the non-CMB data gives a large reduction in 
parameter volume (e.g. Q# models), we run separate chains in- 
stead of importance sampling. 

These runs provide no compelling evidence for deviations 
from the base ACDM model, and indeed, as shown in Table 10 
and Fig. 21, the posteriors for individual extra parameters gen- 
erally overlap the fiducial model within one standard deviation. 
The inclusion of BAO data shrinks further the allowed scope for 
deviation. The parameters of the base ACDM model are rela- 
tively robust to inclusion of additional parameters, but the errors 
on some do broaden significantly when additional degeneracies 
open up, as can be seen in Fig. 21 

The full grid results are available online 28 . Here we summa- 
rize some of the key results, and also consider a few additional 
extensions. 



28 http : / /www . sciops . esa . int/index . php?pro j ect= 
planck&page=Planck_Legacy_Archive 
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Fig. 21. 68% and 95% confidence regions on one-parameter extensions of the base ACDM model for Planck+WP (red) and 
Planck+WP +"& AO (blue). Horizontal dashed lines correspond to the fixed base model parameter value, and vertical dashed lines 
show the mean posterior value in the base model for Planck+WP. 
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Planck+WP Planck+WP+BAO Planck+WP+bighL PZanc£+WP+highL+BAO 



Parameter Best fit 95% limits Best fit 95% limits Best fit 95% limits Best fit 95% limits 



"a 


-0.0105 


-0 037 +0 043 

"■"■"-0.049 


0.0000 


0000+ 0066 

"■""""-0.0067 


-0.0111 


-0 042 +0 043 

"■"^-0.048 


0.0009 


-0 0005 +0 0065 

"•""" -0.0066 


2m v [eV] 


0.022 


< 0.933 


0.002 


< 0.247 


0.023 


< 0.663 


0.000 


< 0.230 


K a 


3.08 


q c i +0.80 
-'-' -0.74 


3.08 


3.40!-? 


3.23 


3 36 +0 - 68 

-'-'"-0.64 


3.22 


3 30 +0 - 54 

-'••'"-0.51 


*p 


0.2583 


0.283!°;°« 


0.2736 


0.283!°;°« 


0.2612 


266 +0 04 ° 


0.2615 


267 +0 03S 

"• z,u/ -0.040 


dn^/dlnk 


-0.0090 


-0.013^1 


-0.0102 


-0 013 +0018 


-0.0106 


-0 015+ 0017 


-0.0103 


-0 014+ 016 

"•""-0.017 


'"0.002 


0.000 


< 0.120 


0.000 


< 0.122 


0.000 


< 0.108 


0.000 


< 0.111 


w 


-1.20 


_i 4.Q+0.65 
lw -0.57 


-1.076 


i i q+0.24 
1 - 1 "'-0.25 


-1.20 


-1 51 +0 - 62 

J -0.53 


-1.109 


-1 13 +0 - 23 

1 - 1J -0.25 



Table 10. Constraints on one-parameter extensions to the base ACDM model. Data combinations all include Planck combined with 
WMAP polarization, and results are shown for combinations with high-/" CMB data and BAO. Note that we quote 95% limits here. 



6.2. Early-Universe physics 

Inflationary cosmology offers elegant explanations of key fea- 
tures of our Universe, such as its large size and near spa- 
tially flat geometry. Within this scenario, the Universe un- 
derwent a brief period of accelerated expansion (Starobinsky 
1979, 1982; Kazanas 1980; Guth 1981; Sato 1981; Linde 1982; 
Albrecht & Steinhardt 1982) during which quantum fluctuations 
were inflated in scale to become the classical fluctuations that 
we see today. In the simplest inflationary models, the primor- 
dial fluctuations are predicted to be adiabatic, nearly scale- 
invariant and Gaussian (Mukhanov & Chibisov 1981; Hawking 
1982; Starobinsky 1982; Guth & Pi 1982; Bardeen et al. 1983), 
in good agreement with CMB observations and other probes of 
large-scale structure. 

Despite this success, the fundamental physics behind in- 
flation is not yet understood and there is no convincing evi- 
dence that rules out alternative scenarios for the early Universe. 
A large number of phenomenological models of inflation, 
some inspired by string theory, have been discussed in the 
literature (see Liddle & Lyth 2000; Bassett et al. 2006; Linde 
2008, for reviews), as well as alternatives to inflation includ- 
ing pre-big bang scenarios (e.g., Gasperini & Veneziano 1993; 
Khouryetal. 2001; Boyle et al. 2004; Creminelli & Senatore 
2007; Brandenberger 2012). Many of these models lead to dis- 
tinctive signatures, such as departures from Gaussianity, isocur- 
vature perturbations, or oscillatory features in the power spec- 
trum, that are potentially observable. The detection of such sig- 
natures would offer valuable information on the physics of the 
early Universe and is one of the main science goals of Planck. 

In this section we discuss basic aspects of the primor- 
dial power spectrum, such as the spectral index, departures 
from a pure power law, limits on tensor modes etc., and 
discuss the implications for inflationary cosmology. Tests 
of more complex models, such as multi-field inflation, are 
discussed in a separate paper (Planck Collaboration XXII 
2013). In Planck Collaboration XXIV (2013), the Planck maps 
are used to constrain possible deviations from Gaussianity 
via measurements of the bispectrum and trispectrum. 
Planck Collaboration XXIII (2013) considers departures 
from statistical isotropy and additional tests of non-Gaussianity. 



6.2.1. Scale dependence of primordial fluctuations 

The primordial fluctuations in the base ACDM model are pa- 
rameterized as a pure power law with a spectral index n s (Eq. 2). 
Prior to Planck, CMB observations have favoured a power 
law index with slope n s < 1, which is expected in simple 
single-field slow-roll inflationary models (see e.g., Mukhanov 
2007 and Eq. 66a below). The final WMAP nine-year data 
give n s = 0.972 + 0.013 at 68% confidence (Hinshaw et al. 
2012). Combining this with damping -tail measurements from 
ACT and SPT data gives n s = 0.968 + 0.009, indicating a de- 
parture from scale invariance at the 3 cr level. The addition of 
BAO data has resulted in a stronger preference for n s < 1 
(Anderson et al. 2013; Hinshaw et al. 2012; Story et al. 2012; 
Sievers et al. 2013). These constraints assume the basic six- 
parameter ACDM cosmological model. Any new physics that 
affects the damping tail of the CMB spectrum, such as additional 
relativistic particles, can alter these constraints substantially and 
still allow a precisely scale-invariant spectrum. 

With Planck, a robust detection of the deviation from scale 
invariance can now be made from a single set of CMB observa- 
tions spanning three decades in scale from i — 2 to I — 2500. 
We find 

n s = 0.959 + 0.007 (68%; Planck+ WP+highL), (61) 

for the base ACDM model, a roughly 6 cr departure from scale 
invariance. This is consistent with the results from previous 
CMB experiments cited above. The statistical significance of this 
result is high enough that the difference between a purely scale 
invariant spectrum can be seen easily in a plot of the power spec- 
trum. Figure 22 shows the Planck spectrum of Fig. 10 plotted as 
e 2 D[ compared to the base ACDM fit with n s = 0.96 (red dashed 
line) and to the best-fit base ACDM cosmology with n s — 1 . The 
n s — 1 model has more power at small scales and is strongly 
excluded by the Planck data. 

The unique contribution of Planck, compared to previous ex- 
periments, is that we are able to show that the departure from 
scale invariance is robust to changes in the underlying theoreti- 
cal model. For example, Figs. 21 and 23 show that the departure 
from scale invariance is not sensitive to the parameterization of 
the primordial fluctuations. Even if we allow a possible running 
of the spectral index (the parameter dn s /dlnk defined in equa- 
tion 2) and/or a component of tensor fluctuations, the Planck 
data favour a tilted spectrum at a high significance level. 
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Fig. 22. The Planck power spectrum of Fig. 10 plotted as £ 2 D( 
against multipole, compared to the best-fit base ACDM model 
with n s - 0.96 (red dashed line). The best-fit base ACDM model 
with « s constrained to unity is shown by the blue line. 



Our extensive grid of models allows us to investigate cor- 
relations of the spectral index with a number of cosmological 
parameters beyond those of the base ACDM model (see Figs. 
21 and 24). As expected, n s is uncorrected with parameters de- 
scribing late-time physics, including the neutrino mass, geom- 
etry, and the equation of state of dark energy. The remaining 
correlations are with parameters that affect the evolution of the 
early Universe, including the number of relativistic species, or 
the helium fraction. This is illustrated in Fig. 24: modifying the 
standard model by increasing the number of neutrinos species, 
or the helium fraction, has the effect of damping the small-scale 
power spectrum. This can be partially compensated by an in- 
crease in the spectral index. However, an increase in the neu- 
trino species must be accompanied by an increased matter den- 
sity to maintain the peak positions. A measurement of the matter 
density from the BAO measurements helps to break this degen- 
eracy. This is clearly seen in the upper panel of Fig. 24, which 
shows the improvement in the constraints when BAO measure- 
ments are added to the Planck+WP +highL likelihood. With the 
addition of BAO measurements we find more than a 3 cr devi- 
ation from n s = 1 even in this extended model, with a best-fit 
value of n s = 0.969 ± 0.010 for varying relativistic species. As 
discussed in Sect. 6.3, we see no evidence from the Planck data 
for non-standard neutrino physics. 

The simplest single-field inflationary models predict that the 
running of the spectral index should be of second order in infla- 
tionary slow -roll parameters and therefore small [dn^/dlnk ~ 
(n s - l) 2 ], typically about an order of magnitude below the 
sensitivity limit of Planck (see e.g., Kosowsky & Turner 1995; 
Baumann et al. 2009). Nevertheless, it is easy to construct in- 
flationary models that have a larger scale dependence (e.g., by 
adjusting the third derivative of the inflaton potential) and so it 
is instructive to use the Planck data to constrain dn s /dlnk. A 
test for dn s /d In k is of particularly interest given the results from 
previous CMB experiments. 

Early results from WMAP suggested a preference for a nega- 
tive running at the 1-2 cr level. In the final 9-year WMAP analy- 
sis no significant running was seen using WMAP data alone, with 
dnjdlnk = -0.019 ± 0.025 (68% confidence; Hinshaw et al. 
2012. Combining WMAP data with the first data releases from 
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Fig. 23. Upper. Posterior distribution for n s for the base ACDM 
model (black) compared to the posterior when a tensor compo- 
nent and running scalar spectral index are added to the model 
(red) Middle: Constraints (68% and 95%) in the n s -dn s /dlnk 
plane for ACDM models with running (blue) and additionally 
with tensors (red). Lower: Constraints (68% and 95%) on n s and 
the tensor-to-scalar ratio rn.002 for ACDM models with tensors 
(blue) and additionally with running of the spectral index (red). 
The dotted line show the expected relation between r and n s for 
a V(<p) cc (p 2 inflationary potential (Eqs. 66a and 66b); here N is 
the number of inflationary e-foldings as defined in the text. The 
dotted line should be compared to the blue contours, since this 
model predicts negligible running. All of these results use the 
Planck+WP +highL data combination. 



38 



Planck Collaboration: Cosmological parameters 



4.2 



3.6 



3.0 



2.4 



1 1 

Planck+WP 
P/anc/c+WP+highL 
P/andc+WP+highL+BAO 




_i_ 



0.92 0.94 0.96 0.98 1.00 1.02 



0.35 



0.30 



0.25 



0.20 



1 1 

Planck+\NP 
P/anc/t+WP+highL 
P/ancfc+WP+highL+BAO 




0.94 



0.96 



0.98 



1.00 



1.02 



Fig. 24. Constraints on n s for ACDM models with non-standard 
relativistic species, A^ff, (upper) and helium fraction, Y P , (lower). 
We show 68% and 95% contours for various data combinations. 
Note the tightening of the constraints with the addition of BAO 
data. 



ACT and SPT, Hinshaw et al. (2012) found a negative running 
at nearly the 2<x level with dnjdlnk = -0.022 ± 0.012 (see 
also Dunkley et al. 2011 and Keisleretal. 2011 for analysis 
of ACT and SPT with earlier data from WMAP). The ACT 
3-year release, which incorporated a new region of sky, gave 
dnjdlnk = -0.003 ± 0.013 (Sievers et al. 2013) when com- 
bined with WMAP 1 year data. With the wide field SPT data at 
150 GHz, a negative running was seen at just over the 2 cr level, 
dnjdlnk = -0.024 ± 0.011 (Hou et al. 2012). 

The picture from previous CMB experiments is therefore 
mixed. The latest WMAP data show a 1 cr trend for a running, 
but when combined with the S12 SPT data, this trend is ampli- 
fied to give a potentially interesting result. The latest ACT data 
go in the other direction, giving no support for a running spectral 
index when combined with WMAP 29 . 

The results from Planck data are as follows (see Figs. 21 and 

23): 



dnjdlnk = -0.013 ± 0.009 (68%; Planck+WP); 



(62a) 



29 The differences between the Planck results reported here and the 
WM4P-7+SPT results (Hou et al. 2012) are discussed in Appendix A. 



dnjdlnk = -0.015 ± 0.009 (68%; Planck+ WP+highL); (62b) 
dnjdlnk = -0.011 ± 0.008 (68%; P/ancAr+lensing 

+WP+highL). (62c) 

The consistency between (62a) and (62b) shows that these re- 
sults are insensitive to modelling of unresolved foregrounds. 
The preferred solutions have a small negative running, but not 
at a high level of statistical significance. Closer inspection of 
the best-fits shows that the change In x 1 when dnjdlnk is in- 
cluded as a parameter comes almost entirely from the low multi- 
pole temperature likelihood. In fact, the fits to the high multipole 
Planck likelihood have a slightly worse x 1 when dnjdlnkis in- 
cluded. The slight preference for a negative running is therefore 
driven by the spectrum at low multipoles I S5 50. The tendency 
for negative running is partly mitigated by including the Planck 
lensing likelihood (Eq. 62c). 

The constraints on dnjdlnk are broadly similar if tensor 
fluctuations are allowed in addition to a running of the spectrum 
(Fig. 23) . Adding tensor fluctuations, the marginalized posterior 
distributions for dnjdlnk give 

dnjdlnk = -0.021 + 0.011 (68%; Planck+WP), (63a) 
dnjdlnk = -0.022 + 0.010 (68%; Planck+WP+highL), (63b) 
dnjdlnk = -0.019 ± 0.010 (68%; Planck+lensing 

+WP+highL). (63c) 

As with Eqs. (62a)-(62c) the tendency to favour negative run- 
ning is driven by the low multipole component of the tempera- 
ture likelihood not by the Planck spectrum at high multipoles. 

This is one of several examples discussed in this section 
where marginal evidence for extensions to the base ACDM 
model are favoured by the TT spectrum at low multipoles. (The 
low multipole spectrum is also largely responsible for the pull of 
the lensing amplitude, Al, to values greater than unity discussed 
in Sect. 5.1). The mismatch between the best-fit base ACDM 
model and the TT spectrum at multipoles t ^ 30 is clearly vis- 
ible in Fig. 1 . The implications of this mismatch are discussed 
further in Sect. 7. 

Beyond a simple running, various extended parameter- 
izations have been developed by e.g., Bridle et al. (2003), 
Shafieloo & Souradeep (2008), Verde & Peiris (2008), and 
Hlozeketal. (2012), to test for deviations from a power-law 
spectrum of fluctuations. Similar techniques are applied to the 
Planck data in Planck Collaboration XXII (2013). 

6.2.2. Tensor fluctuations 

In the base ACDM model, the fluctuations are assumed to 
be purely scalar modes. Primordial tensor fluctuations could 
also contribute to the temperature and polarization power spec- 
tra (e.g., Grishchuk 1975; Starobinsky 1979; Basko & Polnarev 
1980; Crittenden et al. 1993, 1995). The most direct way of test- 
ing for a tensor contribution is to search for a magnetic -type par- 
ity signature via a large-scale B-mode pattern in CMB polar- 
ization (Zaldarriaga & Seljak 1997; Kamionkowski et al. 1997). 
Direct B-mode measurements are challenging as the expected 
signal is small; upper limits measured by BICEP and QUIET 
give 95% upper limits of ro,oo2 < 0.73 and room < 2.8 respec- 
tively (Chiang et al. 2010; QUIET Collaboration et al. 2012) 30 . 



30 As discussed in Planck Collaboration II (2013) and 
Planck Collaboration VI (2013), residual low-level polarization 
systematics in both the LFI and HFI data preclude a Planck S-mode 
polarization analysis at this stage. 
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Measurements of the temperature power spectrum can also 
be used to constrain the amplitude of tensor modes. Although 
such limits can appear to be much tighter than the limits from 
B-mode measurements, it should be borne in mind that such lim- 
its are indirect because they are derived within the context of a 
particular theoretical model. In the rest of this subsection, we 
will review temperature based limits on tensor modes and then 
present the results from Planck. 

Adding a tensor component to the base ACDM model, the 
WMAP 9-year results constrain ro.002 < 0.38 at 95% confidence 
(Hinshaw et al. 2012). Including small-scale ACT and SPT data 
this improves to ro.002 < 0.17, and to ro.002 < 0.12 with the 
addition of BAO data. These limits are degraded substantially, 
however, in models which allow running of the scalar spectral 
index in addition to tensors. For such models, the WMAP data 
give ro.002 < 0.50, and this limit is not significantly improved by 
adding high resolution CMB and BAO data. 

The precise determination of the fourth, fifth and sixth 
acoustic peaks by Planck now largely breaks the degener- 
acy between the primordial fluctuation parameters. For the 
P/anc£+ WP+highL likelihood we find 

ro.002 < 0.11 (95%; no running), (64a) 
ro.002 < 0.26 (95%; including running). (64b) 

As shown in Figs. 21 and 23, the tensor amplitude is weakly cor- 
related with the scalar spectral index; an increase in n s that could 
match the first three peaks cannot fit the fourth and higher acous- 
tic peak in the Planck spectrum. Likewise, the shape constraints 
from the fourth and higher acoustic peaks give a reduction in 
the correlations between a tensor mode and a running in the 
spectral index, leading to significantly tighter limits than from 
previous CMB experiments. These numbers in Eqs. (64a) and 
(64b) are driven by the temperature spectrum and change very 
little if we add non-CMB data such as BAO measurements. The 
Planck limits are largely decoupled from assumptions about the 
late-time evolution of the Universe and are close to the tightest 
possible limits achievable from the temperature power spectrum 
alone (Knox & Turner 1994; Knox 1995). 

These limits on a tensor mode have profound implications 
for inflationary cosmology. The limits translate directly to an up- 
per limit on the energy scale of inflation, 

V. = (1.94 x 10 16 GeV) 4 (r .oo2/0.12) (65) 

(Linde 1983; Lyth 1984), and to the parameters of "large-field" 
inflation models. Slow-roll inflation driven by a power law po- 
tential V(<p) oc (p a offers a simple example of large-field inflation. 
The field values in such a model must necessarily exceed the 
Planck scale mpi , and lead to a scalar spectral index and tensor 
amplitude of 

1 - n s * (a + 2)/2N, (66a) 
r * 4a /N, (66b) 

where N is the number of e-foldings between the end of inflation 
and the time that our present day Hubble scale crossed the infla- 
tionary horizon (see e.g., Lyth & Riotto 1999). The 95% confi- 
dence limits from the Planck data are now close to the predic- 
tions of a — 2 models for N ~ 50-60 e-folds (see Fig. 23). 
Large-field models with quartic potentials (e.g., Linde 1982) are 
now firmly excluded by CMB data. Planck constraints on power- 
law and on broader classes of inflationary models are discussed 
in detail in Planck Collaboration XXIV (2013). Improved lim- 
its on B-modes will be required to further constrain high field 
models of inflation. 



6.2.3. Curvature 

An explanation of the near flatness of our observed Universe 
was one of the primary motivations for inflationary cosmology. 
Inflationary models that allow a large number of e-foldings pre- 
dict that our Universe should be very accurately spatially flat 31 . 
Nevertheless, by introducing fine tunings it is possible to con- 
struct inflation models with observationally interesting open ge- 
ometries (e.g., Linde 1995; Bucheretal. 1995; Linde 1999) or 
closed geometries (Linde 2003). Even more speculatively, there 
has been interest in models with open geometries from consid- 
erations of tunnelling events between metastable vacua within 
a "string landscape" (Freivogel et al. 2006). Observational lim- 
its on spatial curvature therefore offer important additional con- 
straints on inflationary models and fundamental physics. 

CMB temperature power spectrum measurements suffer 
from a well-known "geometrical degeneracy" (Bond et al. 1997; 
Zaldarriaga et al. 1997). Models with identical primordial spec- 
tra, physical matter densities and angular diameter distance to 
the last scattering surface, will have almost identical CMB tem- 
perature power spectra. This is a near perfect degeneracy (see 
Fig. 25) and is broken only via the integrated Sachs-Wolfe (ISW) 
effect on large angular scales and gravitational lensing of the 
CMB spectrum (Stompor & Efstathiou 1999). The geometrical 
degeneracy can also be broken with the addition of probes of 
late time physics, including BAO, Type la supernova, and mea- 
surement of the Hubble constant (e.g., Spergel et al. 2007). 

Recently, the detection of the gravitational lensing of the 
CMB by ACT and SPT has been used to break the geomet- 
rical degeneracy, by measuring the integrated matter potential 
distribution. ACT constrained fi A = 0.61 + 0.29 (68% CL) 
in Sherwin et al. (2011), with the updated analysis in Das et al. 
(2013) giving Q. K = -0.031 + 0.026 (68% CL) (Sievers et al. 
2013). The SPT lensing measurements combined with seven 
year WMAP temperature spectrum improved this limit to Q.k - 
-0.0014 + 0.017 (68 % CL) (van Engelen et al. 2012). 

With Planck we detect gravitational lensing at 
about 26 cr through the 4-point function (Sect. 5.1 and 
Planck Collaboration XVII 2013). This strong detection of 
gravitational lensing allows us to constrain the curvature to 
percent level precision using observations of the CMB alone: 

IQ0D. K = -4.2^J (95%; Planck+ WP+highL); (67a) 
lOOO/f = -LOtJJ (95%; PZanc/fc+lensing 

+ WP+highL). (67b) 

These constraints are improved substantially by the addition 
of BAO data. We then find 

lOOQ/f = -0.05![S;g (95%; P/anofc+WP+highL+BAO), (68a) 
lOOQx = -O.lO^g (95%; PZanc/fc+lensing+WP 

+highL+BAO). (68b) 

These limits are consistent with (and slightly tighter than) the 
results reported by Hinshaw et al. (2012) from combining the 
nine-year WMAP data with high resolution CMB measurements 
and BAO data. We find broadly similar results to Eqs. (68a) and 
(68b) if the Riess et al. (201 1) Hq measurement, or either of the 
SNe compilations discussed in Sect. 5.4, are used in place of the 
BAO measurements. 

3 'The effective curvature within our Hubble radius should then be of 
the order of the amplitude of the curvature fluctuations generated during 
inflation, £l K ~ 0(1(T 5 ). 
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Fig. 25. The PZancA:+ WP+highL data combination (samples; colour-coded by the value of Ho) partially breaks the geometric degen- 
eracy between Q m and Q.\ due to the effect of lensing in the temperature power spectrum. These limits are significantly improved 
by the inclusion of the Planck lensing reconstruction (black contours). Combining also with BAO (right; solid blue contours) tightly 
constrains the geometry to be nearly flat. 



In summary, there is no evidence from Planck for any depar- 
ture from a spatially flat geometry. The results of Eqs. (68a) and 
(68b) suggest that our Universe is spatially flat to an accuracy of 
better than a percent. 

6.3. Neutrino physics and constraints on relativistic 
components 

A striking illustration of the interplay between cosmology and 
particle physics is the potential of CMB observations to con- 
strain the properties of relic neutrinos, and possibly of additional 
light relic particles in the Universe (see e.g., Dodelson et al. 
1996; Hu et al. 1995; Bashinsky & Seljak 2004; Ichikawa et al. 
2005; Lesgourgues & Pastor 2006; Hannestad 2010). In the fol- 
lowing subsections, we present Planck constraints on the mass of 
ordinary (active) neutrinos assuming no extra relics, on the den- 
sity of light relics assuming they all have negligible masses, and 
finally on models with both light massive and massless relics. 

6.3.1. Constraints on the total mass of active neutrinos 

The detection of solar and atmospheric neutrino oscillations 
proves that neutrinos are massive, with at least two species being 
non-relativistic today. The measurement of the absolute neutrino 
mass scale is a challenge for both experimental particle physics 
and observational cosmology. The combination of CMB, large- 
scale structure and distance measurements already excludes a 
large range of masses compared to beta-decay experiments. 
Current limits on the total neutrino mass £ m v (summed over the 
three neutrino families) from cosmology are rather model depen- 
dent and vary strongly with the data combination adopted. The 
tightest constraints for flat models with three families of neutri- 
nos are typically around 0.3 eV (95% CL; e.g., de Putter et al. 
2012). Since 2' w v must be greater than approximately 0.06 eV 
in the normal hierarchy scenario and 0.1 eV in the degener- 
ate hierarchy (Gonzalez -Garcia et al. 2012), the allowed neu- 
trino mass window is already quite tight and could be closed 
further by current or forthcoming observations (Jimenez et al. 
2010; Lesgourgues et al. 2013). 



Cosmological models, with and without neutrino mass, have 
different primary CMB power spectra. For observationally- 
relevant masses, neutrinos are still relativistic at recombina- 
tion and the unique effects of masses in the primary power 
spectra are small. The main effect is around the first acoustic 
peak and is due to the early integrated Sachs-Wolfe (ISW) ef- 
fect; neutrino masses have an impact here even for a fixed red- 
shift of matter-radiation equality (Lesgourgues & Pastor 2012; 
Hall & Challinor 2012; Hou et al. 2012; Lesgourgues et al. 
2013). To date, this effect has been the dominant one in con- 
straining the neutrino mass from CMB data, as demonstrated in 
Hou et al. (2012). As we shall see here, the Planck data move 
us into a new regime where the dominant effect is from gravi- 
tational lensing. Increasing neutrino mass, while adjusting other 
parameters to remain in a high-probability region of parameter 
space, increases the expansion rate at z ^ 1 and so suppresses 
clustering on scales smaller than the horizon size at the non- 
relativistic transition (Kaplinghat et al. 2003; Lesgourgues et al. 
2006). The net effect for lensing is a suppression of the CMB 
lensing potential and, for orientation, by i — 1000 the suppres- 
sion is around 10% in power for 2 m v = 0.66 eV. 

Here we report constraints assuming three species of degen- 
erate massive neutrinos. At the level of sensitivity of Planck, the 
effect of mass splittings is negligible, and the degenerate model 
can be assumed without loss of generality. 

Combining the Planck+WP +highL data, we obtain an upper 
limit on the summed neutrino mass of 



Y m v < 0.66 eV (95%; Planck+ WP+highL). (69) 



The posterior distribution is shown by the solid black curve in 
Fig. 26. To demonstrate that the dominant effect leading to the 
constraint is gravitational lensing, we remove the lensing infor- 
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Fig. 26. Marginalized posterior distributions for 2 m v 
in flat models from CMB data. We show results for 
P/anc£+ WP+highL without (solid black) and with (red) 
marginalization over Al, showing how the posterior is signifi- 
cantly broadened by removing the lensing information from the 
temperature anisotropy power spectrum. The effect of replacing 
the \cw-t temperature and (WMAP) polarization data with a 
t prior is shown in solid blue (P/anc£-lowL+highL+Tprior) 
and of further removing the high-/" data in dot-dashed blue 
(P/ancfc-lowL+rprior). We also show the result of including 
the lensing likelihood with P/a/icA:+ WP+highL (dashed black) 
and P/ancfc-lowL+highL+rprior (dashed blue). 

mation by marginalizing over Al 32 . We see that the posterior 
broadens considerably (see the red curve in Fig. 26) to give 

Y m v < 1 .08 eV [95%; Planck+ WP+highL (Al)], (70) 

taking us back close to the value of 1.3 eV (for Al = 1) from 
the nine-year WMAP data (Hinshaw et al. 2012), corresponding 
to the limit above which neutrinos become non-relativistic be- 
fore recombination. (The resolution of WMAP gives very little 
sensitivity to lensing effects.) 

As discussed in Sect. 5.1, the Planck+ WP+highL data com- 
bination has a preference for high Al- Since massive neutrinos 
suppress the lensing power (like a low A L ) there is a concern 
that the same tensions which drive Al high may give artificially 
tight constraints on £ m v . We can investigate this issue by re- 
placing the low-/ 1 data with a prior on the optical depth (as in 
Sect. 5.1) and removing the high-/' data. Posterior distributions 
with the r prior, and additionally without the high-/' data, are 
shown in Fig. 26 by the solid blue and dot-dashed blue curves, 
respectively. The constraint on £ m Y does not degrade much by 
replacing the low-/' data with the t prior only, but the degra- 
dation is more severe when the high-/" data are also removed: 
2> v <l-31eV(95%CL). 

Including the lensing likelihood (see Sect. 5.1) has a signif- 
icant, but surprising, effect on our results. Adding the lensing 

32 The power spectrum of the temperature anisotropies is predomi- 
nantly sensitive to changes in only one mode of the lensing potential 
power spectrum (Smith et al. 2006). It follows that marginalizing over 
the single parameter Al is nearly equivalent to marginalizing over the 
full amplitude and shape information in the lensing power spectrum as 
regards constraints from the temperature power spectrum. 



likelihood to the Planck+ WP+highL data combination weakens 
the limit on £ ni v , 

Yj m v< °- 85 eV ( 95% ; P^ncfe+lensing+WP+highL), (71) 

as shown by the dashed black curve in Fig. 26. This is representa- 
tive of a general trend that the Planck lensing likelihood favours 
larger £ m v than the temperature power spectrum. Indeed, if we 
use the data combination P/ancfc-lowL+highL+rprior, which 
gives a weaker constraint from the temperature power spectrum, 
adding lensing gives a best-fit away from zero m v = 0.66 eV; 
dashed blue curve in Fig. 26). However, the total x 2 at the 
best-fit is only 0.3 lower than in the best-fitting base model 
(E m v = 0.06 eV). The fit to the lensing data is rather better 
(Ax 2 = -3.2) while the fit to the Planck temperature spec- 
trum (excluding low-f) and high-/' is worse (Ax 2 = 0.4 and 2.6, 
respectively). There are rather large shifts in other cosmolog- 
ical parameters between these best-fit solutions corresponding 
to shifts along the acoustic-scale degeneracy direction for the 
temperature power spectrum. Note that, as well as the change 
in Hq (which falls to compensate the increase in £ m v at fixed 
acoustic scale), n s , u>\, and a> c change significantly keeping the 
lensed temperature spectrum almost constant. These latter shifts 
are similar to those discussed for Al in Sect. 5.1, with non-zero 
2 m v acting like Al < 1. The lensing power spectrum CT is 
lower by 9.1% for the higher-mass best fit at I — 400 and larger 
by 2.1% at I — 40, which is a similar trend to the residuals from 
the best-fit minimal-mass model shown in the bottom panel of 
Fig. 12. Planck Collaboration XVII (2013) explores the robust- 
ness of the CT estimates to various data cuts and foreground- 
cleaning methods. The first (I - 40-85) bandpower is the least 
stable to these choices, although the variations are not statis- 
tically significant. We have checked that excluding this band- 
power does not change the posterior for 2 m Y significantly, as 
expected since most of the constraining power on 2 m v comes 
from the bandpowers on smaller scales. At this stage, it is un- 
clear what to make of this mild preference for high masses from 
the 4-point function compared to the 2-point function. As noted 
in Planck Collaboration XVII (2013), the lensing measurements 
from ACT (Das et al. 2013) and SPT (van Engelen et al. 2012) 
show similar trends to those from Planck where they overlap 
in scale. With further Planck data (including polarization), and 
forthcoming measurements from the full 2500 deg 2 SPT temper- 
ature survey, we can expect more definitive results on this issue 
in the near future. 

Apart from its impact on the early-ISW effect and lensing 
potential, the total neutrino mass affects the angular-diameter 
distance to last scattering, and can be constrained through the 
angular scale of the first acoustic peak. However, this effect is 
degenerate with Q.^ (and so the derived Hq) in flat models and 
with other late-time parameters such as Q. K and w in more gen- 
eral models (Howlett et al. 2012). Late-time geometric measure- 
ments help in reducing this "geometric" degeneracy. Increasing 
the neutrino masses at fixed 6 m increases the angular-diameter 
distance for < z < z* and reduces the expansion rate at low red- 
shift (z $ 1) but increases it at higher redshift. The spherically- 
averaged BAO distance Dy(z) therefore increases with increas- 
ing neutrino mass at fixed G t , and the Hubble constant falls; see 
Fig. 8 of Hou et al. (2012). With the BAO data of Sect. 5.2, we 
find a significantly lower bound on the neutrino mass: 

Y m y < 0.23 eV (95%; PZancfc+WP+highL+BAO). (72) 

The ACDM model with minimal neutrino masses was shown in 
Sect. 5.3 to be in tension with recent direct measurements of Ho 
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which favour higher values. Increasing the neutrino mass will 
only make this tension worse and drive us to artificially tight 
constraints on £ rn v . If we relax spatial flatness, the CMB ge- 
ometric degeneracy becomes three-dimensional in models with 
massive neutrinos and the constraints on £ tn v weaken consider- 
ably to 

( 0.98 eV (95%; PZanc/fc+ WP+highL) 
) m v < { " (73) 

^ ( 0.32 eV (95%; Planck+ WP+highL+BAO). 

6.3.2. Constraints on A^fF 

As discussed in Sect. 2, the density of radiation in the Universe 
(besides photons) is usually parameterized by the effective neu- 
trino number N e ff. This parameter specifies the energy density 
when the species are relativistic in terms of the neutrino tem- 
perature assuming exactly three flavours and instantaneous de- 
coupling. In the Standard Model, A^fF = 3.046, due to non- 
instantaneous decoupling corrections (Mangano et al. 2005). 

However, there has been some mild preference for 
N e g > 3.046 from recent CMB anisotropy measurements 
(Komatsu et al. 2011; Dunkleyetal. 2011; Keisleretal. 2011; 
Archidiacono et al. 201 1; Hinshaw et al. 2012; Hou et al. 2012). 
This is potentially interesting, since an excess could be caused 
by a neutrino/anti-neutrino asymmetry, sterile neutrinos, and/or 
any other light relics in the Universe. In this subsection we dis- 
cuss the constraints on A^ff from Planck in scenarios where the 
extra relativistic degrees of freedom are effectively massless. 

The physics of how N s g is constrained by CMB anisotropies 
is explained in Bashinsky & Seljak (2004), Hou et al. (2011) 
and Lesgourgues et al. (2013). The main effect is that increasing 
the radiation density at fixed 8* (to preserve the angular scales of 
the acoustic peaks) and fixed z eq (to preserve the early-ISW ef- 
fect and so first-peak height) increases the expansion rate before 
recombination and reduces the age of the Universe at recombi- 
nation. Since the diffusion length scales approximately as the 
square root of the age, while the sound horizon varies propor- 
tionately with the age, the angular scale of the photon diffusion 
length, On, increases, thereby reducing power in the damping tail 
at a given multipole. Combining Planck, WMAP polarization and 
the high-f experiments gives 

N eB = 3.36!o;g4 (95%; Planck +WP+highL). (74) 

The marginalized posterior distribution is given in Fig. 27 (black 
curve). 

Increasing A^ff at fixed ft, and z e q necessarily raises the ex- 
pansion rate at low redshifts too. Combining CMB with distance 
measurements can therefore improve constraints (see Fig. 27) al- 
though for the BAO observable fd rag /.Dv(z) the reduction in both 
rdrag and Dy(z) with increasing N e s partly cancel. With the BAO 
data of Sect. 5.2, the A^ff constraint is tightened to 

AW = 3.30^ (95%; P/anc*+ WP+highL+BAO). (75) 

Our constraints from CMB alone and CMB+BAO are compati- 
ble with the standard value A^ff = 3.046 at the 1 cr level, giving 
no evidence for extra relativistic degrees of freedom. 

Since A^ff is positively correlated with Ho, the tension be- 
tween the Planck data and direct measurements of Hq in the base 
ACDM model (Sect. 5.3) can be reduced at the expense of high 
AW . The marginalized constraint is 

N eB = 3.621om ( 95% ; Planck+WP+highL+Ho). (76) 
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Fig. 27. Marginalized posterior distribution of N e ff for 
P/ancfc+WP+highL (black) and additionally BAO (blue), 
the Ho measurement (red), and both BAO and Ho (green). 

For this data combination, the x 2 f° r the best-fitting model al- 
lowing A'efF to vary is lower by 5.0 than for the base A^ff = 3.046 
model. The Hq fit is much better, with Ax 2 = -4.0, but there 
is no strong preference either way from the CMB. The low-/' 
temperature power spectrum does mildly favour the high N e g 
model (Ax 2 = -1.6) since A^ff is positively correlated with « s 
(see Fig. 24) and increasing n s reduces power on large scales. 
The rest of the Planck power spectrum is agnostic [Ax 2 = -0.5), 
while the high-^ experiments mildly disfavour high A^ft in our 
fits (Ax 2 = 1.3). Further including the BAO data pulls the cen- 
tral value downwards by around 0.5 cr (see Fig. 27): 

N e « = 3.52^ (95%; PZa«c/fc+WP+highL+// +BAO). (77) 

The x 2 at the best-fit for this data combination (N e g = 3.37) 
is lower by 3.6 than the best-fitting A^ff = 3.046 model. While 
the high N ef[ best-fit is preferred by Planck+WP (Ax 2 = -3.3) 
and the Hq data (Ax 2 = -2.8 giving an acceptable x 2 = 2.4 
for this data point), it is disfavoured by the high-f CMB data 
(Ax 2 = 2.0) and slightly by BAO (Ax 2 = 0.4). We conclude 
that the tension between direct Hq measurements and the CMB 
and BAO data in the base ACDM can be relieved at the cost of 
additional neutrino-like physics, but there is no strong preference 
for this extension from the CMB damping tail. 

Throughout this subsection, we have assumed that all the 
relativistic components parameterized by A^ff consist of ordi- 
nary free-streaming relativistic particles. Extra radiation com- 
ponents with a different sound speed or viscosity parame- 
ter (Hu 1998) can provide a good fit to ^re-Planck CMB 
data (Archidiacono et al. 2013), but are not investigated in this 
paper. 

6.3.3. Simultaneous constraints on N eff and either £ m v or 

v, sterile 

It is interesting to investigate simultaneous contraints on N e g and 
2 my, since extra relics could coexist with neutrinos of size- 
able mass, or could themselves have a mass in the eV range. 
Joint constraints on N s g and £ m v have been explored sev- 
eral times in the literature. These two parameters are known 
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to be partially degenerate when large-scale structure data are 
used (Hannestad & Raffelt 2004; Crotty et al. 2004), but their 
impact in the CMB is different and does not lead to significant 
correlations. 

Joint constraints on A^ff and £ m v are always model- 
dependent: they vary strongly with assumptions about how the 
total mass is split between different species (and they would also 
be different for models in which massive species have chem- 
ical potentials or a non-thermal phase-space distribution). We 
present here Planck constraints for two different models and de- 
scribe the scenarios that motivate them. 

First, as in the previous subsection we assume that the three 
active neutrinos share a mass m v — £ m v /3, and may coexist 
with extra massless species contributing to A^ff. In this model, 
when N e ff is greater than 3.046, AN eS = A^ eff - 3.046 gives the 
density of extra massless relics with arbitrary phase-space dis- 
tribution. When N e ff < 3.046, the temperature of the three active 
neutrinos is reduced accordingly, and no additional relativistic 
species are assumed. In this case, the CMB constraint is 



as m 



N a = 3 29 +0 - 67 
2>v < 0.60 eV 



(95%; P/anc/fc+ WP+highL). 



(78) 



These bounds tighten somewhat with the inclusion of BAO data, 
as illustrated in Fig. 28; we find 



N, 



— Q OO+0.54 
e ff J.JZ._ Q52 

£m„ < 0.28 eV 



(95%; Planck+ WP+highL+BAO). (79) 



We see that the joint constraints do not differ very much from 
the bounds obtained when introducing these parameters sepa- 
rately. The physical effects of neutrino masses and extra rela- 
tivistic relics are sufficiently different to be resolved separately 
at the level of accuracy of Planck. 

In the second model, we assume the existence of one mas- 
sive sterile neutrino, in addition to the two massless and one 
massive active neutrino of the base model. The active neutrino 
mass is kept fixed at 0.06 eV. In particle physics, this assump- 
tion can be motivated in several ways. For example, there has 
recently been renewed interest in models with one light sterile 
neutrino in order to explain the MiniBoone anomaly reported 
in Aguilar-Arevalo et al. (2012), as well as reactor and Gallium 
anomalies (Giunti et al. 2013). The statistical significance of 
these results is marginal and they should not be over-interpreted. 
However, they do motivate investigating a model with three ac- 
tive neutrinos and one heavier sterile neutrino with mass m ster ji e . 
If the sterile neutrino were to thermalize with the same tempera- 
ture as active neutrinos, this model would have N e s ~ 4. 

Since we wish to be more general, we assume that the ex- 
tra eigenstate is either: (i) thermally distributed with an arbi- 
trary temperature T s ; or (ii) distributed proportionally to ac- 
tive neutrinos with an arbitrary scaling factor Xs in which the 
scaling factor is a function of the active-sterile neutrino mix- 
ing angle. This second case corresponds the Dodelson-Widrow 
scenario (Dodelson & Widrow 1994). The two cases are in fact 
equivalent for cosmological observables and do not require sep- 
arate analyses (Colombi et al. 1996; Lesgourgues et al. 2013). 
Sampling the posterior with flat priors on N e g and m ster iie would 
not be efficient, since in the limit of small temperature T s , or 
small scaling factor x%, the mass would be unbounded. Hence we 
adopt a flat prior on the "effective sterile neutrino mass" defined 



v ff s teriie = (94. 1 w v , steri i e ) eV 33 . In the case of a thermally- 
distributed sterile neutrino, this parameter is related to the true 

mass via 



m eff = (T IT } 3 m tham ' dl = (AN B-") 3/4 m* ermaI 

m - sterile y 1 si 1 v) w sterile K^&fiJ '"sterile • 



(80) 



Here, recall that T v = (4/1 1)^ 3 !^ is the active neutrino temper- 
ature in the instantaneous-decoupling limit and that the effective 
number is defined via the energy density, A/V e ff = (T s /T v )*. In 
the Dodelson-Widrow case the relation is given by 



eff _ DW 

m v, sterile _ As'« s terile ■ 



(81) 



with AN e s = AV For a thermalized sterile neutrino with temper- 
ature T v (i.e., the temperature the active neutrinos would have if 
there were no heating at electron-positron annihilation), corre- 
sponding to AA'erF = 1, the three masses are equal to each other. 
Assuming flat priors on AW and m eff , ., with f7i eff , ., < 

° r cl1 v, sterile v, sterile 

3eV, we find the results shown in Fig. 28. The physical mass, 
'"sterile 31 m me thermal scenario is constant along the dashed lines 
in the figure and takes the indicated value in eV. The physical 
mass, mPy., , in the Dodelson-Widrow scenario is constant on 

sterile 

the dotted lines. For low AT e ff the physical mass of the neutrinos 
becomes very large, so that they become non-relativistic well be- 
fore recombination. In the limit in which the neutrinos become 
non-relativistic well before any relevant scales enter the horizon, 
they will behave exactly like cold dark matter, and hence are 
completely unconstrained within the overall total constraint on 
the dark matter density. For intermediate cases where the neutri- 
nos become non-relativistic well before recombination they be- 
have like warm dark matter. The approach to the massive limit 
gives the tail of allowed models with large »Jy ff sterile and low A^tr 
shown in Fig. 28, with increasing m v s stelile being compensated 
by decreased Q c /z 2 to maintain the total level required to give the 
correct shape to the CMB power spectrum. 

For low »iy tt sterile and AA^s away from zero the physical neu- 
trino mass is very light, and the constraint becomes similar to 
the massless case. The different limits are continuously con- 
nected, and given the complicated shape seen in Fig. 28 it is 
clearly not appropriate to quote fully marginalized parameter 
constraints that would depend strongly on the assumed upper 
limit on m t ., . Instead we restrict attention to the case where 

v, sterile 

the physical mass is 'w*^™ 1 < 10 eV, which roughly defines the 
region where (for the CMB) the particles are distinct from cold 
or warm dark matter. Using the Planck+ WP+highL (abbreviated 
to CMB below) data combination, this gives the marginalized 
one-parameter constraints 



A^eff < 3.91 
m eff t ., <0.59eV 

v, stenle 



(95%; CMB for m* e ™ al < 10 eV) . (82) 



Combining further with BAO these tighten to 



JVeff < 3.80 
>i efl \ , <0.42eV 

v, stenle 



(95%; CMB+BAO for m* e ™ al < 10 eV) . 



(83) 



These bounds are only marginally compatible with a fully ther- 
malized sterile neutrino (N e s ~ 4) with sub-eV mass m^^f ~ 
< 0.5 eV that could explain the oscillation anomalies. 



//?' 



elf 



, sterile 



33 The factor of 94. 1 eV here is the usual one in the relation between 
physical mass and energy density for non-relativistic neutrinos with 
physical temperature T v . 
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Fig. 28. Left: 2D joint posterior distribution between N e g and £ m v (the summed mass of the three active neutrinos) in models with 
extra massless neutrino-like species. Right: Samples in the ^Veff-'«y ft sterile plane, colour-coded by Q. c h 2 , in models with one massive 
sterile neutrino family, with effective mass fw^ ff sterile , and the three active neutrinos as in the base ACDM model. The physical mass 
of the sterile neutrino in the thermal scenario, 'w*"™ 1 , is constant along the grey dashed lines, with the indicated mass in eV. The 
physical mass in the Dodelson-Widrow scenario, wHSL, is constant along the dotted lines (with the value indicated on the adjacent 
dashed lines). 



The above contraints are also appropriate for the Dodelson- 
Widrow scenario, but for a physical mass cut of tnP?^ < 20 eV. 

The thermal and Dodelson-Widrow scenarios considered 
here are representative of a large number of possible models that 
have recently been investigated in the literature (Hamann et al. 
2011; Diamanti et al. 2012; Archidiacono et al. 2012; 
Hannestad et al. 2012). 

6.4. Big bang nucleosynthesis 

Observations of light elements abundances created during big 
bang nucleosynthesis (BBN) provided one of the earliest preci- 
sion tests of cosmology and were critical in establishing the ex- 
istence of a hot big bang. Up-to-date accounts of nucleosynthe- 
sis are given by Iocco et al. (2009) and Steigman (2012). In the 
standard BBN model, the abundance of light elements (parame- 
terized by T^ BN = 4« He /"b for helium-4 and y™ N = 10 5 n D /n H 
for deuterium, where n, is the number density of species i) can 
be predicted as a function of the baryon density «b, the number 
of relativistic degrees of freedom parameterized by N e s, and of 
the lepton asymmetry in the electron neutrino sector. Throughout 
this subsection, we assume for simplicity that lepton asymmetry 
is too small to play a role at BBN. This is a reasonable assump- 
tion, since Planck data cannot improve existing constraints on 
the asymmetry 34 . We also assume that there is no significant en- 



34 A primordial lepton asymmetry could modify the outcome of BBN 
only if it were very large (of the order of 10~ 3 or bigger). Such a large 
asymmetry is not motivated by particle physics, and is strongly con- 
strained by BBN. Indeed, by taking into account neutrino oscillations 
in the early Universe, which tend to equalize the distribution function 
of three neutrino species, Mangano et al. (2012) derived strong bounds 
on the lepton asymmetry. CMB data cannot improve these bounds, as 
shown by Castorina et al. (2012); an exquisite sensitivity to N e g would 
be required. Note that the results of Mangano et al. (2012) assume that 
N e g departs from the standard value only due to the lepton asymmetry. 
A model with both a large lepton asymmetry and extra relativistic relics 
could be constrained by CMB data. However, we will not consider such 
a contrived scenario in this paper. 



tropy increase between BBN and the present day, so that our 
CMB constraints on the baryon-to-photon ratio can be used to 
compute primordial abundances. 

To calculate the dependence of Tp BN and y™ N on the 
parameters co^, and N e g, we use the accurate public code 
PArthENoPE (Pisanti et al. 2008), which incorporates values 
of nuclear reaction rates, particle masses and fundamental 
constants, and an updated estimate of the neutron lifetime 
(r n = 880.1s; Beringer et al. 2012). Experimental uncertain- 
ties on each of these quantities lead to a theoretical error for 
Yp BN (o) h , N e ff) andy B p N (wb,A r e ff). For helium, the error is dom- 
inated by the uncertainty in the neutron lifetime, leading to 35 
cr(Tp BN ) = 0.0003. For deuterium, the error is dominated by 
uncertainties in several nuclear rates, and is estimated to be 
o-(Jdp N ) = °- 04 (Serpico et al. 2004). 

These predictions for the light elements can be confronted 
with measurements of their abundances, and also with CMB data 
(which is sensitive to co^, N e s, and Yp). We shall see below that 
for the base cosmological model with N e g = 3.046 (or even for 
an extended scenario with free N e ff) the CMB data predict the 
primordial abundances, under the assumption of standard BBN, 
with smaller uncertainties than those estimated for the measured 
abundances. Furthermore, the CMB predictions are consistent 
with direct abundance measurements. 

6.4.1. Observational data on primordial abundances 

The observational constraint on the primordial helium-4 frac- 
tion used in this paper is Y* BN = 0.2534 ± 0.0083 (68% CU) 
from the recent data compilation of Averetal. (2012), based 
on spectroscopic observations of the chemical abundances in 
metal-poor H n regions. The error on this measurement is domi- 
nated by systematic effects that will be difficult to resolve in the 
near future. It is reassuring that the independent and conserva- 



35 Serpico et al. (2004) quotes cr(r° BN ) = 0.0002, but since that 
work, the uncertainty on the neutron lifetime has been re-evaluated, 
from cr(T n ) = 0.8 s to ct(t„) = 1.1 s Beringer et al. (2012). 
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tive method presented in Mangano & Serpico (201 1) leads to an 
upper bound for Tp BN that is consistent with the above estimate. 
The recent measurement of the proto-Solar helium abundance 
by Serenelli & Basu (2010) provides an even more conservative 
upper bound, Yp < 0.294 at the 2 cr level. 

For the primordial abundance of deuterium, data points show 
excess scatter above the statistical errors, indicative of system- 
atic errors. The compilation presented in Iocco et al. (2009), 
based on data accumulated over several years, gives y^p N = 
2.87 + 0.22 (68% CL). Pettini & Cooke (2012) report an accu- 
rate deuterium abundance measurement in the z = 3.04984 low- 
metallicity damped Ly a system in the spectrum of QSO SDSS 
J1419+0829, which they argue is particularly well suited to deu- 
terium abundance measurements. These authors find y™ N = 
2.535 ± 0.05 (68% CL), a significantly tighter constraint than 
that from the Iocco et al. (2009) compilation. The Pettini-Cooke 
measurement is, however, a single data point, and it is im- 
portant to acquire more observations of similar systems to as- 
sess whether their error estimate is consistent with possible 
sources of systematic error. We adopt a conservative position 
in this paper and compare both the Iocco et al. (2009) and the 
Pettini & Cooke (2012) measurements to the CMB predictions 

We consider only the 4 He and D abundances in this paper. 
We do not discuss measurements of 3 He abundances since these 
provide only an upper bound on the true primordial 3 He frac- 
tion. Likewise, we do not discuss lithium. There has been a long 
standing discrepancy between the low lithium abundances mea- 
sured in metal-poor stars in our Galaxy and the predictions of 
BBN. At present it is not clear whether this discrepancy is caused 
by systematic errors in the abundance measurements, or has an 
"astrophysical" solution (e.g., destruction of primordial lithium) 
or is caused by new physics (see Fields 2011, for a recent re- 
view). 

6.4.2. Planck predictions of primordial abundances in 
standard BBN 

We first restrict ourselves to the base cosmological model, with 
no extra relativistic degrees of freedom beyond ordinary neutri- 
nos (and a negligible lepton asymmetry), leading to iV e ff = 3.046 
(Mangano et al. 2005). Assuming that standard BBN holds, and 
that there is no entropy release after BBN, we can compute 
the spectrum of CMB anisotropies using the relation Tp(wb) 
given by PArthENoPE. This relation is used as the default 
in the grid of models discussed in this paper; we use the 
CosmoMC implementation developed by Hamann et al. (2008). 
The Planck+ WP+highL fits to the base ACDM model gives the 
following estimate of the baryon density, 

w b = 0.02207 ± 0.00027 (68%; P/anc/fc+WP+highL), (84) 

as listed in Table 5. In Fig. 29, we show this bound together 
with theoretical BBN predictions for Tp BN (dL» b ) and y B p N (Wb)- 
The bound of Eq. (84) leads to the predictions 

Yp(a> h ) = 0.24725 ± 0.00032, (85a) 
y DP (w b ) = 2.656 ± 0.067, (85b) 

where the errors here are 68% and include theoretical errors that 
are added in quadrature to those arising from uncertainties in 
Mb- (The theoretical error dominates the total error in the case 
of lp.) For helium, this prediction is in very good agreement 

36 Note that, throughout this paper, our quoted CMB constraints on 
all parameters do not include the theoretical uncertainty in the BBN 
relation (where used). 
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Fig. 29. Predictions of standard BBN for the primordial abun- 
dance of 4 He (top) and deuterium (bottom), as a function of 
the baryon density. The width of the green stripes corresponds 
to 68% uncertainties on nuclear reaction rates. The horizontal 
bands show observational bounds on primordial element abun- 
dances compiled by various authors, and the red vertical band 
shows the Planck+ WP+highL bounds on cj\, (all with 68% er- 
rors). BBN predictions and CMB results assume N e $ = 3.046 
and no significant lepton asymmetry. 

with the data compilation of Averetal. (2012), with an error 
that is 26 times smaller. For deuterium, the CMB+BBN pre- 
diction lies midway between the best-fit values of Iocco et al. 
(2009) and Pettini & Cooke (2012), but agrees with both at ap- 
proximately the 1 <x level. These results strongly support stan- 
dard BBN and show that within the framework of the base 
ACDMmodel, Planck observations lead to extremely precise 
predictions of primordial abundances. 

6.4.3. Estimating the helium abundance directly from Planck 
data 

In the CMB analysis, instead of fixing Yp to the BBN predic- 
tion, Tp BN (wb), we can relax any BBN prior and let this pa- 
rameter vary freely. The primordial helium fraction has an influ- 
ence on the recombination history and affects CMB anisotropies 
mainly through the redshift of last scattering and the dif- 
fusion damping scale (Hu et al. 1995; Trotta & Hansen 2004; 
Ichikawa & Takahashi 2006; Hamann et al. 2008). Extending 
the base ACDM model by adding Y P as a free parameter with 
a flat prior in the range [0.1, 0.5], we find 

Yp = 0.266 ± 0.021 (68%; Planck+ WP+highL). (86) 

Constraints in the ip-Wb plane are shown in Fig. 30. This figure 
shows that the CMB data have some sensitivity to the helium 
abundance. In fact, the error on the CMB estimate of Yp is only 
2.7 times larger than the direct measurements of the primordial 
helium abundance by Averet al. (2012). The CMB estimate of 
Yp is consistent with the observational measurements adding fur- 
ther support in favour of standard BBN. 

6.4.4. Extension to the case with extra relativistic relics 

We now consider the effects of additional relativistic degrees of 
freedom on photons and ordinary neutrinos (obeying the stan- 
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Fig. 30. Constraints in the (o^—Yp plane from CMB and 
abundance measurements. The CMB constraints are for 
PZa«c/fc+ WP+highL (red 68% and 95% contours) in ACDM 
models with Tpallowed to vary freely. The horizontal band 
shows observational bounds on 4 He compiled by Aver et al. 
(2012) with 68% errors, while the grey region at the top of 
the figure delineates the conservative 95% upper bound inferred 
from Solar helium abundance by Serenelli & Basu (2010). The 
green stripe shows the predictions of standard BBN for the pri- 
mordial abundance of 4 He as a function of the baryon density 
(with 68% errors on nuclear reaction rates). Both BBN predic- 
tions and CMB results assume N e g = 3.046 and no significant 
lepton asymmetry. 

dard model of neutrino decoupling) by adding N e g as a free pa- 
rameter. In the absence of lepton asymmetry, we can predict the 
BBN primordial abundances as a function of the two parameters 
a>b and N e g. 

Figure 31 shows the regions in the u>b-N e g plane preferred 
by primordial abundance measurements, and by the CMB data if 
the standard BBN picture is correct. The regions allowed by the 
abundance measurements are defined by the y 1 statistic 



X ' 2 (u b ,N eB ) = 



\y((o b ,N eS ) - Jobs] 2 



cr 



obs 



(T. 



(87) 



theory 



where y(u\,, N e g) is the BBN prediction for either Yp or yp,p, the 
quantity y b s is the observed abundance, and the two errors in the 
denominator are the observational and theoretical uncertainties. 
Figure 31 shows the edges of the 68% preferred regions in the 
Wb-AW plane, given by x 2 = X 2 mia + 2.3. 

For the CMB data, we fit a cosmological model with seven 
free parameters (the six parameters of the base ACDM model, 
plus Neg) to the PZancA:+ WP+highL data, assuming that the pri- 
mordial helium fraction is fixed by the standard BBN prediction 
Yp{a)b, Neg). Figure 31 shows the joint 68% and 95% confidence 
contours in the co\,-N e g plane. The preferred regions in this plane 
from abundance measurements and the CMB agree remarkably 
well. The CMB gives approximately three times smaller error 
bars than primordial abundance data on both parameters. 

We can derive constraints on N e g from primordial element 
abundances and CMB data together by combining their likeli- 
hoods. The CMB-only confidence interval for N e g is 
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Fig. 31. Constraints in the a)\,-N e g plane from the CMB and 
abundance measurements. The blue stripes shows the 68% 
confidence regions from measurements of primordial element 
abundances assuming standard BBN: 4 He bounds compiled 
by Averetal. (2012); and deuterium bounds complied by 
Iocco et al. (2009) or measured by Pettini & Cooke (2012). We 
show for comparison the 68% and 95% contours inferred from 
PZancA:+ WP+highL, when N e g is left as a free parameter in the 
CMB analysis (and Yp is fixed as a function of a»b and N e g ac- 
cording to BBN predictions). These constraints assume no sig- 
nificant lepton asymmetry. 

When combined with the data reported respectively by 
Averetal. (2012), Iocco et al. (2009), and Pettini & Cooke 
(2012), the 68% confidence interval becomes 



A eff = 



3.41 ±0.30, Yp (Aver et al.), 

3.43 ± 0.34, y DP (Iocco et al.), 

3.02 ± 0.27, y DP (Pettini and Cooke). 



(89) 



Aeff = 3.36 ± 0.34 (68%; Planck+ WP+highL). 



(88) 



Since there is no significant tension between CMB and primor- 
dial element results, all these bounds are in agreement with the 
CMB-only analysis. The small error bar derived from combining 
the CMB with the Pettini & Cooke (2012) data point shows that 
further deuterium observations combined with Planck data have 
the potential to pin down the value of A^ff to high precision. 

6.4.5. Simultaneous constraints on both N e g and Y P 

In this subsection, we discuss simultaneous constraints on both 
A^eff and Yp by adding them to the six parameters of the base 
ACDM model. Both N e g and Y P have an impact on the damp- 
ing tail of the CMB power spectrum by altering the ratio A:" 1 /r„, 
where k D l is the photon diffusion length at last scattering and 
r» is the sound horizon there. There is thus an approximate de- 
generacy between these two parameters along which the ratio is 
nearly constant. The extent of the degeneracy is limited by the 
characteristic phase shift of the acoustic oscillations that arises 
due to the free streaming of the neutrinos (Bashinsky & Seljak 
2004). As discussed by Hou et al. (2011), the early ISW effect 
also partly breaks the degeneracy, but this is less important than 
the effect of the phase shifts. 

The joint posterior distribution for N e g and Yp from 
the Planck+ WP+highL likelihood is shown in Fig. 32, with 
each MCMC sample colour-coded by the value of the 
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Fig. 32. 2D joint posterior distribution for N e g and Yp with both 
parameters varying freely, determined from Planck+ WP+highL 
data. Samples are colour-coded by the value of the angular ra- 
tio 6d/0*, which is constant along the degeneracy direction. The 
Neg—Yp relation from BBN theory is shown by the dashed curve. 
The vertical line shows the standard value N e g = 3.046. The 
region with Yp > 0.294 is highlighted in grey, delineating the re- 
gion that exceeds the 2 <x upper limit of the recent measurement 
of initial Solar helium abundance (Serenelli & Basu 2010), and 
the blue horizontal region is the 68% confidence region from 
the Aver et al. (2012) compilation of 4 He measurements. 

observationally-relevant angular ratio 6 D /9 t oc (k D r t )~ l . The 
main constraint on N e g and Yp comes from the precise measure- 
ment of this ratio by the CMB, leaving the degeneracy along the 
constant 9 D /0 t direction. The relation between N e g and Y P from 
BBN theory is shown in the figure by the dashed curve 37 . The 
standard BBN prediction with A^ff = 3.046 is contained within 
the 68% confidence region. The grey region is for Yp > 0.294 
and is the 2 <x conservative upper bound on the primordial he- 
lium abundance from Serenelli & Basu (2010). Most of the sam- 
ples are consistent with this bound. The inferred estimates of N e g 
and Yp from the Planck+ WP+highL data are 



Yp 



3.33 



+0.59 
■0.83 



0.254 



+0.041 
-0.033 



(68%; Planck+ WP+highL), 
(68%; Planck+ WP+highL). 



(90a) 
(90b) 



With Yp allowed to vary, N e g is no longer tightly constrained 
by the value of 8u/0*. Instead, it is constrained, at least in part, 
by the impact that varying A^sf has on the phase shifts of the 
acoustic oscillations. As discussed in Hou et al. (2012), this ef- 
fect explains the observed correlation between N e g and 6 m , which 
is shown in Fig. 33. The correlation in the ACDM+A^ff model 
is also plotted in the figure showing that the N e g—Yp degeneracy 
combines with the phase shifts to generate a larger dispersion in 
6* in such models. 

6.5. Dark energy 

A major challenge for cosmology is to elucidate the nature of the 
dark energy driving the accelerated expansion of the Universe. 
Perhaps the most straightforward explanation is that dark en- 
ergy is a cosmological constant. An alternative is dynamical dark 



37 For constant N e g, the variation due to the uncertainty in the baryon 
density is too small to be visible, given the thickness of the curve. 
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Fig. 33. 2D joint posterior distribution between N e g and 6* for 
ACDM models with variable N e g (blue) and variable N e g and Yp 
(red). Both cases are for Planck+ WP+highL data. 

energy (Wetterich 1988; Ratra & Peebles 1988; Caldwell et al. 
1998), usually based on a scalar field. In the simplest models, 
the field is very light, has a canonical kinetic energy term and 
is minimally coupled to gravity. In such models the dark energy 
sound speed equals the speed of light and it has zero anisotropic 
stress. It thus contributes very little to clustering. We shall only 
consider such models in this subsection. 

A cosmological constant has an equation of state w = pip - 
-1, while scalar field models typically have time varying w with 
W > — 1 . The analysis performed here is based on the "parameter- 
ized post-Friedmann" (PPF) framework of Hu & Sawicki (2007) 
and Hu (2008) as implemented in camb (Fang et al. 2008b, a) and 
discussed earlier in Sect. 2. This allows us to investigate both re- 
gions of parameter space in which w < — 1 (sometimes referred 
to as the "phantom" domain) and models in which w changes 
with time. 

Figure 34 shows the marginalized posterior distributions for 
w for an extension of the base ACDM cosmology to models with 
constant w. We present results for Planck+WP and in combi- 
nation with SNe or BAO data. (Note that adding in the high-/' 
data from ACT and SPT results in little change to the posteriors 
shown in Fig. 34.) As expected, the CMB alone does not strongly 
constrain w, due to the two-dimensional geometric degeneracy 
in these models. We can break this degeneracy by combining 
the CMB data with lower redshift distance measures. Adding in 
BAO data tightens the constraints substantially, giving 



w = -1.13 



+0.24 
-0.25 



(95%; Planck+WP+BAO), 



in good agreement with a cosmological constant (w 
Using supernovae data leads to the constraints 



-1.09 ± 0.17 (95%; Planck+ WP+Union2.1), 
(95%; Planck+WP+SNLS), 



-1.13 



+0.13 
-0.14 



(91) 

(92a) 
(92b) 



The combination with SNLS data favours the phantom domain 
(w < -1) at 2 cr, while the Union2.1 compilation is more consis- 
tent with a cosmological constant. 

If instead we combine Planck+WP with the Riess et al. 
(201 1) measurement of Hq, we find 



w = -1.24 



+0.18 
-0.19 



(95%; Planck+WP+H ), 



(93) 
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Fig. 34. Marginalized posterior distributions for the dark en- 
ergy equation of state parameter w (assumed constant), for 
Planck+WP alone (green) and in combination with SNe data 
(SNSL in blue and the Union2.1 compilation in red) or BAO 
data (black). A flat prior on w from -3 to -0.3 was as- 
sumed and, importantly for the CMB-only constraints, the prior 
[20, 100]kms _1 Mpc~' on Hq. The dashed grey line indicates 
the cosmological constant solution, w = - 1 . 



which is in tension with w — - 1 at more than the 2 a level. 

The results in Eqs. (91-93) reflect the tensions between the 
supplementary data sets and the Planck base ACDM cosmology 
discussed in Sect. 5. The BAO data are in excellent agreement 
with the Planck base ACDM model, so there is no significant 
preference for w + — 1 when combining BAO with Planck. In 
contrast, the addition of the Ho measurement, or SNLS SNe data, 
to the CMB data favours models with exotic physics in the dark 
energy sector. These trends form a consistent theme throughout 
this section. The SNLS data favours a lower Q. in the ACDM 
model than Planck, and hence larger dark energy density today. 
The tension can be relieved by making the dark energy fall away 
faster in the past than for a cosmological constant, i.e., w < — 1. 

The constant w models are of limited physical interest. If 
w + -1 then it is likely to change with time. To investigate 
this we consider the simple linear relation in Eq. (4), w(a) = 
wo + w a (l - a), which has often been used in the literature 
(Chevallier & Polarski 2001; Linder 2003). This parameteriza- 
tion approximately captures the low-redshift behaviour of light, 
slowly-rolling minimally-coupled scalar fields (as long as they 
do not contribute significantly to the total energy density at early 
times) and avoids the complexity of scanning a large number of 
possible potential shapes and initial conditions. The dynamical 
evolution of w(a) can lead to distinctive imprints in the CMB 
(Caldwell et al. 1998) which would show up in the Planck data. 

Figure 35 shows contours of the joint posterior distribution in 
the wo-w fl plane using P/ancA:+WP+BAO data (colour-coded ac- 
cording to the value of Hq). The points are coloured by the value 
of Ho, which shows a clear variation with wq and w a reveal- 
ing the three-dimensional nature of the geometric degeneracy in 




Fig. 35. 2D marginalized posterior distribution for wq and w a 
for Planck+WP+BAO data. The contours are 68% and 95%, 
and the samples are colour-coded according to the value of // () . 
Independent flat priors of -3 < wq < -0.3 and -2 < w a < 2 
are assumed. Dashed grey lines show the cosmological constant 
solution Wq = — 1 and w a = 0. 



such models. The cosmological constant point (wq, w a ) - (—1,0) 
lies within the 68% contour and the marginalized posteriors for 
wo and w a are 



= -1.04t 7 J 

w a < 1.32 



(95%; P/anc/fc+WP+BAO), 
(95%; Planck+WP+BAO). 



(94a) 
(94b) 



Including the Hq measurement in place of the BAO data moves 
(wo, w a ) away from the cosmological constant solution towards 
negative w a at just under the 2 cr level. 

Figure 36 shows likelihood contours for (wo,w a ), now 
adding SNe data to Planck. As discussed in detail in Sect. 5, 
there is a dependence of the base ACDM parameters on the 
choice of SNe data set, and this is reflected in Fig. 36. The re- 
sults from the WancA:+ WP+Union2. 1 data combination are in 
better agreement with a cosmological constant than those from 
the Planck+WP+SNLS combination. For the latter data combi- 
nation, the cosmological constant solution lies on the 2 cr bound- 
ary of the (wo, w a ) distribution. 

Dynamical dark energy models might also give a non- 
negligible contribution to the energy density of the Universe 
at early times. Such early dark energy (EDE; Wetterich 2004) 
models may be very close to ACDM recently, but have a non- 
zero dark energy density fraction, £2 e , at early times. Such mod- 
els complement the (wq, w a ) analysis by investigating how much 
dark energy can be present at high redshifts. EDE has two main 
effects: it reduces structure growth in the period after last scat- 
tering; and it changes the position and height of the peaks in the 
CMB spectrum. 

The model we adopt here is that of Doran & Robbers (2006): 



Q de (fl) 



Q e (l 



,-3w 



Q0 e + Q;) fl 3,,„ 



+ fie(l -a- iw ") 



(95) 



It requires two additional parameters to those of the base ACDM 
model: £2 e , the dark energy density relative to the critical den- 
sity at early times (assumed constant in this treatment); and the 
present-day dark energy equation of state parameter wq. Here flJJ, 
is the present matter density and £2^ e = 1 -QjJ, is the present dark 
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Fig. 36. 2D marginalized posterior distributions for wq and 
w a , for the data combinations Plcmck+WP+BAO (grey), 
Planck+ WP+Union2.1 (red) and Planck+ WP+SNLS (blue). 
The contours are 68% and 95%, and dashed grey lines show the 
cosmological constant solution. 

energy abundance (for a fiat Universe). Note that the model of 
Eq. (95) has dark energy present over a large range of redshifts; 
the bounds on D. e can be substantially weaker if dark energy is 
only present over a limited range of redshifts (Pettorino et al. 
2013). The presence or absence of dark energy at the epoch of 
last scattering is the dominant effect on the CMB anisotropies 
and hence the constraints are insensitive to the addition of low 
redshift supplementary data such as BAO. 

The most precise bounds on EDE arise from the analysis 
of CMB anisotropies (Doran et al. 2001; Caldwell et al. 2003; 
Calabrese et al. 2011; Reichardt et al. 2012; Sievers et al. 
2013; Houetal. 2012; Pettorino et al. 2013). Using 
Planck+ WP+highL, we find 



O e < 0.009 (95%; Planck+ WP+highL). 



(96) 



(The limit for Planck+WP is very similar: Q e < 0.010.) These 
bounds are consistent with and improve the recent ones of 
Houetal. (2012), who give O f < 0.013 at 95% CL, and 
Sievers et al. (2013), who find Q e < 0.025 at 95% CL. 

In summary, the results on dynamical dark energy (except for 
those on early dark energy discussed above) are dependent on 
exactly what supplementary data are used in conjunction with 
the CMB data. (Planck lensing does not significantly improve 
the constraints on the models discussed here.) Using the direct 
measurement of Ho, or the SNLS SNe sample, together with 
Planck we see preferences for dynamical dark energy at about 
the 2 cr level reflecting the tensions between these data sets and 
Planck in the ACDM model. In contrast, the BAO measurements 
together with Planck give tight constraints which are consistent 
with a cosmological constant. Our inclination is to give greater 
weight to the BAO measurements and to conclude that there is 
no strong evidence that the dark energy is anything other than a 
cosmological constant. 

6.6. Dark matter annihilation 

Energy injection from dark matter (DM) annihilation can 
change the recombination history and affect the shape of 
the angular CMB spectra (Chen & Kamionkowski 2004; 



Padmanabhan & Finkbeiner 2005; Zhang et al. 2006; 
Mapelli et al. 2006). As recently shown in several papers 
(see e.g., Galli et al. 2009, 2011; Giesen et al. 2012; Hutsi et al. 
2011; Natarajan 2012) CMB anisotropies offer an opportunity 
to constrain DM annihilation models. 

High-energy particles injected in the high-redshift thermal 
gas by DM annihilation are typically cooled down to the keV 
scale by high energy processes; once the shower has reached 
this energy scale, the secondary particles produced can ion- 
ize, excite or heat the thermal gas (Shull & van Steenberg 1985; 
Valdes et al. 2010); the first two processes modify the evolution 
of the free electron fraction jt e , while the third affects the tem- 
perature of the baryons. 

The rate of energy release, dE/dt, per unit volume by a relic 
annihilating DM particle is given by 



(z) = 2gp c 2 C 2 Q2(l+ Z )Va„n(z), 



dE 
dt 

where p 



(97) 



Pann(z) = f(Z) 



is, in principle, a function of redshift z, defined as 

(crv) 



(98) 



where (crv) is the thermally averaged annihilation cross-section, 
m x is the mass of the DM particle, p c is the critical density of 
the Universe today, g is a degeneracy factor equal to 1/2 for 
Majorana particles and 1 /4 for Dirac particles (in the following, 
constraints will refer to Majorana particles), and the parameter 
f(z) indicates the fraction of energy which is absorbed overall 
by the gas at redshift z. We note that the presence of the brackets 
in (crv) denote a thermal average over the velocity distribution 
of particles. 

In Eq. (98), the factor f(z) depends on the details of the 
annihilation process, such as the mass of the DM particle and 
the annihilation channel (see e.g., Slatyer et al. 2009). The func- 
tional shape of f(z) can be taken into account using gen- 
eralized parameterizations (Finkbeiner et al. 2012; Hutsi et al. 
2011). However, as shown in Galli et al. (2011), Giesen et al. 
(2012), and Finkbeiner et al. (2012) it is possible to neglect the 
redshift dependence of f(z) to first approximation, since current 
data shows very little sensitivity to variations of this function. 
The effects of DM annihilation can therefore be well parameter- 
ized by a single constant parameter, p ann , that encodes the de- 
pendence on the properties of the DM particles. 

We compute here the theoretical angular power in the pres- 
ence of DM annihilations, by modifying the RECFAST routine 
in the camb code as in Galli et al. (2011) and by making use 
of the package CosmoMC for Monte Carlo parameter estimation. 
We checked that we obtain the same results by using the CLASS 
Boltzmann code (Lesgourgues 2011a) and the Monte Python 
package (Audren et al. 2012), with DM annihilation effects cal- 
culated either by RECFAST or HyRec (Ali-Haimoud & Hirata 
201 1), as detailed in Giesen et al. (2012). Besides p Rnn , we sam- 
ple the parameters of the base ACDM model and the fore- 
ground/nuisance parameters described in Sect. 4. 

From Planck+WP we find 



< 5.4 x 10 



b m 3 s" 



kg" 1 (95; Planck+WP). 



(99) 



This constraint is weaker than that found from the full 
WMAP9 temperature and polarization likelihood, p. dan < 1 .2 x 
10" 6 m 3 s _1 kg _1 because the Planck likelihood does not yet in- 
clude polarization information at intermediate and high multi- 
poles. In fact, the damping effect of DM annihilation on the 
CMB temperature power spectrum is highly degenerate with 
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other cosmological parameters, in particular with the scalar 
spectral index and the scalar amplitude, as first shown by 
Padmanabhan & Finkbeiner (2005). As a consequence, the con- 
straint on the scalar spectral index is significantly weakened 
when is allowed to vary, n s = 0.984*9qB, to be compared to 
the constraint listed in Table 2 for the base ACDM cosmology, 
n s = 0.9603 ± 0.0073. 

These degeneracies can be broken by polarization data. The 
effect of DM annihilation on polarization is in fact an overall 
enhancement of the amplitude at large and intermediate scales, 
and a damping at small scales (see e.g., Fig. 1 in Galli et al. 2009 
or Fig. 3 in Giesen et al. 2012). We thus expect the constraint to 
improve significantly with the forthcoming Planck polarization 
data release. We verified that adding BAO, HST or highL data to 
Planck+WP improves the constraints only marginally, as these 
datasets are not able to break the degeneracy between p. dnn and 
n s . 

On the other hand, we observe a substantial improvement 
in the constraints when we combine the Planck+WP data with 
the Planck lensing likelihood data. For this data combination we 
find an upper limit of 

p mm < 3.1xl0~ 6 m 3 s" 1 kg" 1 (95%; f7anc/t+lensing+WP).(100) 

The improvement over Eq. (99) comes from the constraining 
power of the lensing likelihood on A s and n s , that partially breaks 
the degeneracy with p^n- 

Our results are consistent with previous work and show no 
evidence for DM annihilation. Future release of Planck polar- 
ization data will help to break the degeneracies which currently 
limit the accuracy of the constraints presented here. 

6. 7. Constraints on a stochastic background of primordial 
magnetic fields 

Large-scale magnetic fields of the order of a few pG observed 
in galaxies and galaxy clusters may be the product of the am- 
plification during structure formation, of primordial magnetic 
seeds (Ryu et al. 2012). Several models of the early Universe 
predict the generation of primordial magnetic fields (hereafter 
PMF), either during inflation or during later phase transitions 
(see Widrow 2002 and Widrow et al. 2012 for reviews). 

PMF have an impact on cosmological perturbations and 
in particular on CMB anisotropy angular power spectra 
(Subramanian 2006), that can be used to constrain the PMF 
amplitude. In this section we will derive the constraints from 
Planck data on a stochastic background of PMF. We are mainly 
interested in constraints from CMB temperature anisotropies. 
Therefore, we will not consider the effect of Faraday ro- 
tation on CMB polarization anisotropies (Kosowsky & Loeb 
1996; Kosowsky et al. 2005) nor non-Gaussianities associ- 
ated with PMF (Brown & Crittenden 2005; Caprini et al. 2009; 
Seshadri & Subramanian 2009; Trivedi et al. 2010). We will re- 
strict the analysis reported here to the non-helical case. 

A stochastic background of PMF is modelled as a fully in- 
homogeneous component whose energy-momentum tensor is 
quadratic in the fields. We assume the usual magnetohydrody- 
namics limit, in which PMF are frozen and the time evolu- 
tion is simply given by the dilution with cosmological expan- 
sion, B(k,Tj) = B{k)/a{j]) 2 . We model the PMF with a simple 
power-law power spectrum: P B (k) = Ak" B , with a sharp cut off 
at the damping scale k D , as computed in Jedamzik et al. (1998) 
and Subramanian & Barrow (1998), to model the suppression of 
PMF on small scales. 



It is customary to specify the amplitude of the PMF power 
spectrum with B^, the root mean square of the field smoothed 
over length scale A, defined such that 

B\ = J ^-e- k2A2 P B (k). (101) 

Given our assumed model and conventions, PMF are fully de- 
scribed by two parameters: the smoothed amplitude By, and the 
spectral index n B . For specificity, we set A — 1 Mpc and hence 
use B\ as the parameter. 

The components of the energy momentum tensor of PMF 
source all types of linear cosmological perturbations, i.e., scalar, 
vector and tensor. In particular, the source terms are given by 
the magnetic energy density and anisotropic stress for scalar 
magnetized perturbations, whereas vector and tensor modes are 
sourced only by the magnetic anisotropic stress. In addition, both 
scalar and vector perturbations are affected by the Lorentz force; 
PMF induce a Lorentz force on baryons modifying their evolu- 
tion and in particular their velocity, but during the tight coupling 
regime between matter and radiation the Lorentz force also has 
an indirect effect on photons. 

For the computation of magnetized angular power spec- 
tra, we use the analytic approximations for the PMF energy- 
momentum tensor components given in Paoletti & Finelli 
(201 1). We consider here the regular mode for magnetized scalar 
perturbations, with the initial conditions of Paoletti et al. (2009) 
and Shaw & Lewis (2010) (which describes the singular passive 
mode, depending on the generation time of PMF). 

Previous analyses show that the main impact of PMF on 
the CMB anisotropy angular power spectrum is at small angular 
scales, well into the Silk damping regime. The dominant mode 
is the magnetized vector mode which peaks at ( ~ 2000-3000 
(Mack et al. 2002; Lewis 2004). The scalar magnetized mode is 
the dominant PMF contribution on large and intermediate angu- 
lar scales (Giovannini & Kunze 2008; Finelli et al. 2008). The 
tensor contribution is always subdominant with respect to the 
other two and it is negligible for the purposes of this analysis. 

We include the scalar and vector magnetized contributions 
to the angular power spectrum within the MCMC analysis in or- 
der to derive the constraints on the PMF amplitude and spectral 
index using Planck TT data. We vary the magnetic parameters 
B[/nG and ns, in addition to the other cosmological parameters 
of the base ACDM cosmology (this analysis assumed massless 
neutrinos, rather than the default value of a single eigenstate of 
mass 0.06 eV used in the rest of this paper). We adopt as prior 
ranges for the parameters [0 , 10] for Z?i/nG and [-2.99 , 3] for 
the spectral index n B . The lower bound n B > -3 is necessary to 
avoid infrared divergences in the PMF energy momentum tensor 
correlators. 

We perform analyses with Planck+WP and 
Planck+WP +highL likelihood combinations. Results are 
shown in Fig. 37. We find that the cosmological parameters 
are in agreement with those estimated assuming no PMF, 
confirming that the magnetic parameters are not degenerate 
with the cosmological parameters of the base ACDM model. 
The constraints on PMF with Planck+WP likelihood are 
B\ < 4.1 nG, with a preference for negative spectral indices 
at the 95% confidence level. These limits are improved using 
the Planck+ WP+highL likelihood to B x < 3.4 nG with n B < 
preferred at the 95% confidence level. The new constraints 
are consistent with, and slightly tighter, than previous limits 
based on combining WMAP-1 data with high resolution CMB 
data (see e.g. Paoletti & Finelli 2011; Shaw & Lewis 2012; 
Paoletti & Finelli 2012). 
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Table 11. Constraints on the cosmological parameters of the base ACDM model with the addition of a varying fine-structure 
constant. We quote ±1 <x errors. Note that for WMAP there is a strong degeneracy between Hq and a, which is why the error on 
ajaiQ is much larger than for Planck. 



Planck+WP Planck+WP+BAO WMAP-9 



Q. b h 2 0.02206 ± 0.00028 0.02220 ± 0.00025 0.02309 ± 0.00130 

Slji 2 0.1174±0.0030 0.1161 +0.0028 0.1148 + 0.0048 

t 0.095 + 0.014 0.097 + 0.014 0.089 + 0.014 

H 65.2+ 1.8 66.7+ 1.1 73.9+ 10.9 

n s 0.975 + 0.012 0.969 + 0.012 0.973 + 0.014 

log(10 10 A s ) 3.106 + 0.029 3.100 + 0.029 3.090+ 0.039 

a/a 0.9936 + 0.0043 0.9989 + 0.0037 1.008 + 0.020 




50 0.96 0.98 1.00 1.02 1.04 020 0.96 0.98 1.00 1.02 1.04 0.950 0.975 1.000 1.025 1.050 

a/aQ ol/olq ol/olq 

Fig. 38. Left: Likelihood contours (68% and 95%) in the a/a -H plane for the WMAP-9 (red), Planck+WP (blue), Planck+WP+H 
(purple), and Planck+WP+BAO (green) data combinations. Middle: As left, but in the a/ao—Clbh 2 plane. Right: Marginalized 
posterior distributions of a/ao for these data combinations. 




Fig. 37. Constraints on the PMF amplitude obtained with 
Planck+WP (red line) and P/a/ic£+ WP+highL (black line). 



6.8. Constraints on variation of the fine-structure constant 

The ACDM model assumes the validity of General Relativity 
on cosmological scales, as well as the physics of the standard 
model of particle physics. One possible extension, which may 
have motivations in fundamental physics, is to consider varia- 
tions of dimensionless constants. Such variations can be con- 
strained through tests on astrophysical scales (Uzan 2003, 201 1). 

A number of physical systems have been used, spanning dif- 
ferent time scales, to set constraints on variations of the funda- 
mental constants. These range from atomic clocks in the lab- 
oratory at a redshift z = to BBN at z ~ 10 8 . However, 
apart from the claims of varying a based on high resolution 
quasar absorption-line spectra (Webb et al. 2001; Murphy et al. 
2003) 38 , there is no other evidence for time-variable fundamen- 
tal constants. 

CMB temperature anisotropies have been used extensively 
to constrain the variation of fundamental constants over cos- 
mic timescales. The temperature power spectrum is sensitive 
to the variation of the fine-structure constant a, the electron-to- 
proton mass ratio ju, and the gravitational constant a„ = Gm^lfic. 
A variation of G can affect the Friedmann equation, and also 
raises the issue of consistency in the overall theory of grav- 
ity. However, a variation of the non-gravitational constants (a 
and m e ) is more straightforward, mostly inducing a modifica- 
tion of the interaction between light and atoms (shifts in the en- 
ergy levels and binding energy of hydrogen and helium). This 
induces a modification of the ionization history of the Universe. 
In particular, a variation of a modifies the redshift of recombi- 
nation through the shift in the energy levels and the Thomson 



See however Srianand et al. (2004, 2007). 
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scattering cross-section. An increase in a induces a shift of the 
position of the first acoustic peak, which is inversely propor- 
tional to the sound horizon at last scattering. The larger redshift 
of last scattering also produces a larger early ISW effect, and 
hence a higher amplitude of the first acoustic peak. Finally, an 
increase in a decreases diffusive damping at high multipoles. 
For earlier studies of varying constants using the CMB (see e.g., 
Kaplinghat et al. 1999; Avelino et al. 2000; Nakashima et al. 
2008, 2010; Menegoni et al. 2009; Landau & Scoccola 2010). 

The analysis presented here focuses solely on the time vari- 
ation of the fine-structure constant a, in addition to the param- 
eters of the base ACDM model, using a modified form of the 
RECFAST recombination code (Hannestad 1999; Martins et al. 
2004; Rocha et al. 2004). Selected results are given in Table 1 1, 
which compares parameter constraints from Planck and from our 
own analysis of the full WMAP-9 TT, TE and EE likelihood. 
From CMB data alone, Planck improves the constraints from a 
2% variation in a to about 0.4%. Planck thus improves the limit 
by a factor of around five, while the constraints on the param- 
eters of the base ACDM model change very little with the ad- 
dition of a time-varying a. These results are in good agreement 
with earlier forecasts (Rocha et al. 2004). 

Given the apparent tension between the base ACDM param- 
eters from Planck and direct measurements of Ho discussed in 
Sect. 5.3), we include further information from the Ho prior and 
BAO data (see Sect. 5.2). Fig. 38 compares the constraints in the 
(a/ao, Ho) and (a/ao, £\/z 2 ) planes and also shows the marginal- 
ized posterior distribution of a/ao for the various data combina- 
tions. 

The constraint on a is slightly improved by including the 
BAO data (via a tightening of the parameters of the base ACDM 
model). Note that the central value of the prior on Ho is outside 
the 95% confidence region, even for the Planck+WP+Ho combi- 
nation. Adding a varying a does not resolve the tension between 
direct measurements of Ho and the value determined from the 
CMB. 

In summary, Planck data improve the constraints on a/ao, 
with respect to those from WMAP-9 by a factor of about five. Our 
analysis of Planck data limits any variation in the fine-structure 
constant from z ~ 10 3 to the present day to be less than approxi- 
mately 0.4%. 



7. Discussion and conclusions 39 

The most important conclusion from this paper is the excel- 
lent agreement between the Planck temperature power spec- 
trum at high multipoles with the predictions of the base ACDM 
model. The base ACDM model also provides a good match to 
the Planck power spectrum of the lensing potential, Ct , and to 
the TE and EE power spectra at high multipoles. 

The high statistical significance of the Planck detection of 
gravitational lensing of the CMB leads to some interesting sci- 
ence conclusions using Planck data alone. For example, grav- 
itational lensing breaks the "geometrical degeneracy" and we 
find that the geometry of the Universe is consistent with spatial 
flatness to percent-level precision using CMB data alone. The 
Planck lensing power spectrum also leads to an interesting con- 
straint on the reionization optical depth of t = 0.089 + 0.032, 
independent of CMB polarization measurements at low multi- 
poles. 



39 Unless otherwise stated, we quote 68% confidence limits in this 
section for the Planck+WP+highh data combination. 



The parameters of the base ACDM model are determined to 
extremely high precision by the Planck data. For example, the 
scalar spectral index is determined as n s = 0.9585 ± 0.0070, 
a 6cr deviation from exact scale invariance. Even in the base 
ACDM model, we find quite large changes in some param- 
eters compared to previous CMB experiments 40 . In particu- 
lar, from Planck we find a low value of the Hubble constant, 
H) - (67.3 ± 1.2)kms~' Mpc~', and a high matter density, 
£2 m = 0.315 ± 0.016. If we accept that the base ACDM model 
is the correct cosmology, then as discussed in Sect. 5 Planck is 
in tension with direct measurements of the Hubble constant (at 
about the 2.5 cr level) and in mild tension with the SNLS Type 
la supernova compilation (at about the 2 cr level). For the base 
ACDM model, we also find a high amplitude for the present-day 
matter fluctuations, <x 8 = 0.828 + 0.012, in agreement with pre- 
vious CMB experiments. This value is higher than that inferred 
from counts of rich clusters of galaxies, including our own anal- 
ysis of Planck cluster counts (Planck Collaboration XX 2013). 

One possible interpretation of these tensions is that some 
sources of systematic error are not completely understood in 
some astrophysical measurements. The fact that the Planck re- 
sults for the base ACDM model are in such good agreement with 
BAO data, which are based on a simple geometrical measure- 
ment, lends support to this view. An alternative explanation is 
that the base ACDM model is incorrect. In summary, at high 
multipoles, the base ACDM cosmology provides an excellent 
fit to the spectra from Planck, ACT and SPT (for all frequency 
combinations), as illustrated in Figs. 7, 8 and 9, but the parame- 
ters derived from the CMB apparently conflict with some types 
of astrophysical measurement. 

Before summarizing our results on extensions to the base 
ACDM model, it is worth making some remarks on foreground 
modelling and the impact of this modelling on our error esti- 
mates. The addition of CMB data at high multipoles helps to 
constrain the model of unresolved foregrounds, in particular, the 
contribution from "minor" components, such as the kinetic SZ, 
which are poorly constrained from Planck alone. For the base 
ACDM model, the cosmological parameters are not limited by 
foreground modelling, as illustrated 41 in Fig. 4. As discussed in 
Appendix B, foreground modelling becomes more important in 
analysing extended CDM models, particularly those that have 
strong parameter degeneracies that are broken only via precision 
measurements of the damping tail in the CMB spectrum. As a 
crude measure of the importance of foreground modelling, we 
can compare parameter values with and without inclusion of the 
ACT and SPT data at high multipoles. A large shift in parameter 
values indicates a possible sensitivity to foreground modelling, 
and so any such result should be treated with caution. We have 
thus normally adopted the Planck+WP+highL likelihood com- 
bination as offering the most reliable results for extensions to the 
base ACDM cosmology. 

From an analysis of an extensive grid of models, we find 
no strong evidence to favour any extension to the base ACDM 
cosmology, either from the CMB temperature power spectrum 
alone, or in combination with the Planck lensing power spec- 
trum and other astrophysical data sets. 



40 The particular tension in the case of the SPT S12 results is dis- 
cussed in detail in Appendix A. 

4 'Even in the restricted case of the base ACDM model, parameters 
can shift as a result of small changes to the theoretical assumptions. An 
example is given in Sect. 3.2, where we show that changing from our 
default assumption of 2 My = 0.06 eV to 2 ni v = 0, causes an upward 
shift of 0.4 cr in the value of Hq. 
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Maximum multipole moment, £, 



Fig. 39. Left: Planck TT spectrum at low multipoles with 68% ranges on the posteriors. The "rainbow" band show the best fits to 
the entire Planck+WP likelihood for the base ACDM cosmology, colour-coded according to the value of the scalar spectral index 
n s . Right: Limits (68% and 95%) on the relative amplitude of the base ACDM fits to the Planck+WP likelihood fitted only to the 
Planck TT likelihood over the multipole range 2 < I < ^ max . 



We find the following notable results using CMB data alone: 

- The deviation of the scalar spectral index from unity is ro- 
bust to the addition of tensor modes and to changes in the 
matter content of the Universe. For example, adding a tensor 
component we find n s = 0.9600 ± 0.0072, a 5.5 cr departure 
from n s = 1 . 

- A 95% upper limit on the tensor-to-scalar ratio of ro.002 < 
0.11. The combined contraints on n s and ro.002 are on the 
borderline of compatibility with single-field inflation with a 
quadratic potential (Fig. 23). 

- A 95% upper limit on the summed neutrino mass of £ m v < 
0.66 eV. 

- A determination of the effective number of neutrino-like rel- 
ativistic degrees of freedom of N s « = 3.36 ±0.34, compatible 
with the standard value of 3.046. 

- The results from Planck are consistent with the results of 
standard big bang nucleosynthesis. In fact, combining the 
CMB data with the most recent results on the deuterium 
abundance, leads to the constraint N e « = 3.02 + 0.27, again 
compatible with the standard value of 3.046. 

- New limits on a possible variation of the fine-structure 
constant, dark matter annihilation and primordial magnetic 
fields. 

We also find a number of marginal (around 2(f) results, 
perhaps indicative of internal tension within the Planck data. 
Examples include the preference of the (phenomenological) 
lensing parameter for values greater than unity (Al = 1.23+0.11; 
Eq. 44) and for negative running (dn s /dlnk = -0.015±0.09; Eq. 
62b). In Planck Collaboration XXII (2013), the Planck data indi- 
cate a preference for anti-correlated isocurvature modes and for 
models with a truncated power spectrum on large scales. None 
of these results have a decisive level of statistical significance, 
but they can all be traced to an unusual aspect of the tempera- 
ture power spectrum at low multipoles. As can be seen in Fig. 
1, and on an expanded scale in the left-hand panel of Fig. 39, 
the measured power spectrum shows a dip relative to the best-fit 
base ACDM cosmology in the multipole range 20 ^ t 30 and 
an excess at I — 40. The existence of "glitches" in the power 
spectrum at low multipoles was noted by the WMAP team in the 



first-year papers (Hinshaw et al. 2003; Spergel et al. 2003) and 
acted as motivation to fit an inflation model with a step-like fea- 
ture (Peiris et al. 2003). Similar investigations have been carried 
out by a number of authors, (see e.g., Mortonson et al. 2009, and 
references therein). At these low multipoles, the Planck spec- 
trum is in excellent agreement with the WMAP nine-year spec- 
trum (Planck Collaboration XV 2013), so it is unlikely that any 
of the features such as the low quadrupole or "dip" in the multi- 
pole range 20-30 are caused by instrumental effects or Galactic 
foregrounds. These are real features of the CMB anisotropics. 

The Planck data, however, constrain the parameters of the 
base ACDM model to such high precision that there is little re- 
maining flexibility to fit the low-multipole part of the spectrum. 
To illustrate this point, the right-hand panel of Fig. 39 shows the 
68% and 95% limits on the relative amplitude of the base ACDM 
model (sampling the chains constrained by the full likelihood) 
fitted only to the Planck TT likelihood over the multipole range 
2 < I < ( mdx . From multipoles ^ max « 25 to multipoles { max « 
35, we see more than a 2 cr departure from values of unity. (The 
maximum deviation from unity is 2.7 cr at I — 30.) It is difficult 
to know what to make of this result, and we present it here as a 
"curiosity" that needs further investigation. The Planck temper- 
ature data are remarkably consistent with the predictions of the 
base ACDM model at high multipoles, but it is also conceivable 
that the ACDM cosmology fails at low multipoles. There are 
other indications, from both WMAP and Planck data for "anoma- 
lies" at low multipoles (Planck Collaboration XXIII 2013), that 
may be indicative of new physics operating on the largest scales 
in our Universe. Interpretation of large-scale anomalies (includ- 
ing the results shown in Fig. 39) is difficult in the absence of a 
theoretical framework. The problem here is assessing the role of 
a posteriori choices, i.e., that inconsistencies attract our atten- 
tion and influence our choice of statistical test. Nevertheless, we 
know so little about the physics of the early Universe that we 
should be open to the possibility that there is new physics be- 
yond that assumed in the base ACDM model. Irrespective of the 
interpretation, the unusual shape of the low multipole spectrum 
is at least partly responsible for some of the 2 cr effects seen in 
the analysis of extensions to the ACDM model discussed in Sect. 
6. 
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Supplementary information from astrophysical data sets has 
played an important role in the analysis of all previous CMB 
experiments. For Planck the interpretation of results combined 
with non-CMB data sets is not straightforward (as a consequence 
of the tensions discussed in Sect. 5). For the base ACDM model, 
the statistical power of the Planck data is so high that we find 
very similar cosmological parameters if we add the Riess et al. 
(2011) constraint on Hq, or either of the two SNe samples, to 
those derived from the CMB data alone. In these cases, the solu- 
tions simply reflect the tensions discussed in Sect. 5, for exam- 
ple, including the Hq measurement with the Planck+WP likeli- 
hood we find Ho = (68.6 ± 1.2)kms~' Mpc~', discrepant with 
the direct measurement at the 2.2 cr level. 

The interpretation becomes more complex for extended 
models where astrophysical data is required to constrain param- 
eters that cannot be determined accurately from CMB measure- 
ments alone. As an example, it is well known that CMB data 
alone provide weak constraints on the dark energy equation of 
state parameter w (see Fig. 34). The addition of BAO data to 
the CMB data gives a tight constraint of w — —1.13 ± 0.12 42 . 
However, adding the SNLS SNe data gives w = -1.135 ± 0.069 
and adding the Hq measurement gives w = -1.244 ± 0.095. 
Adding either of the two data sets which show tension with the 
CMB measurements for the base ACDM model, draws the solu- 
tions into the phantom domain (w < - 1) at about the 2 <x level. In 
contrast, if we use the BAO data in addition to the CMB, we find 
no evidence for dynamical dark energy; these data are compati- 
ble with a cosmological constant, as assumed in the base ACDM 
model. 

The impact of additional astrophysical data is particularly 
complex in our investigation of neutrino physics (Sect. 6.3). 
We will use the effective number of relativistic degrees of free- 
dom, Af e ff as an illustration. From the CMB data alone, we find 
N eS = 3.36 + 0.34. Adding BAO data gives N eS = 3.30 ± 0.27. 
Both of these values are consistent with the standard value of 
3.046. Adding the Hq measurement to the CMB data gives 
N e g = 3.62±0.25 and relieves the tension between the CMB data 
and Hq at the expense of new neutrino-like physics (at around 
the 2.3 cr level). It is possible to alleviate the tensions between 
the CMB, BAO, H and SNLS data by invoking new physics 
such as an increase in N e ff. However, none of these cases are 
favoured significantly over the base ACDM model by the Planck 
data (and they are often disfavoured). Any preference for new 
physics comes almost entirely from the astrophysical data sets. 
It is up to the reader to decide how to interpret such results, but it 
is simplistic to assume that all astrophysical data sets have accu- 
rately quantified estimates of systematic errors. We have there- 
fore tended to place greater weight on the CMB and BAO mea- 
surements in this paper rather than on more complex astrophys- 
ical data. 

Our overall conclusion is that the Planck data are remark- 
ably consistent with the predictions of the base ACDM cosmol- 
ogy. However, the mismatch with the temperature spectrum at 
low multipoles, evident in Figs. 1 and 39, and the existence of 
other "anomalies" at low multipoles, is possibly indicative that 
the model is incomplete. The results presented here are based on 
a first, and relatively conservative, analysis of the Planck data. 
The 2014 data release will use data obtained over the full mis- 
sion lifetime of Planck, including polarization data. It remains 
to be seen whether these data, together with new astrophysical 



42 The addition of the Planck lensing measurements tightens this fur- 
ther to w = - 1 .08!" 1L- 



data sets and CMB polarization measurements, will offer any 
convincing evidence for new physics. 
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Fig. A.l. Best-fit temperature power spectra for the base ACDM 
model for Planck+ WP+highL and for SPT+WMAP-7 (up- 
per panel). The parameters defining these models are listed 
in Table A.l. The lower panel shows the residuals of the 
SPT+WMAP-7 best fit with respect to the Planck best fit. 
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Fig.A.2. The acoustic-scale distance ratio r s /Dy(z) divided by 
the distance ratio of the best-fit WMAP-7 + SPT base ACDM 
cosmology of Table A.l. The points are colour-coded as fol- 
lows: green star (6dF); pink squares (SDSS DR7 as analysed 
by Percival etal. 2010); black star (SDSS DR7 as analysed by 
Padmanabhan et al. 2012); blue cross (BOSS DR9); and blue 
circles (WiggleZ). Error bars show 1 cr errors on the data points. 
The grey band shows the ±1 cr range allowed by the WMAP- 
7+SPT data. 



Appendix A: Comparison of the Planck and 
WMAP7+SPT S12 base ACDM cosmologies 

The main result of this paper is the precise agreement between 
the Planck power spectra and the predictions of the base ACDM 
model. The exploration of the grid of extended CDM models 
described in Sect. 6 shows that the Planck data provide no evi- 
dence for any new physics beyond that incorporated in the base 
ACDM model. This conclusion differs from the recent results of 
S12 and (Hou et al. 2012, , hereafter H12) based on an analysis 
of the temperature power spectrum computed from SPT obser- 
vations at 150 GHz combined with WMAP-7. The S12 and H12 
results extend the earlier work of Keisleret al. (2011) from an 
analysis of 790 deg 2 of sky to a total field area of 2540 deg 2 . 
S12 discuss constraints on the base ACDM model and sim- 
ple extensions such as ACDM models with non-zero curvature. 
More exotic models are discussed by HI 2, who point out that 
the SPT power spectrum declines slightly faster with increasing 
multipole than expected in the base ACDM model. When com- 
bined with additional data, H12 find tentative evidence for new 
physics. Although they cannot find a unique explanation for their 
observations, certain models appear to be favoured over the base 
ACDM model. For example, by combining WMAP-7 and SPT 
they find that a running spectral index is favoured over a con- 
stant spectral index at the 2.2 cr level, a significance that grows 
to 2.7 cr when BAO and Hq data are added. 

The analysis of Sect. 6 shows that neither this (nor other) ex- 
tensions to the base ACDM model are favoured by the Planck 
data. It is therefore important to understand why the Planck and 
WMAP-7+SPT S12 data give different results. This is the pur- 
pose of this Appendix. 

We begin by comparing the parameters of the base ACDM 
model of this paper with those of S12. These are listed in Table 
A. 1 . One can see immediately that the best-fit value of D. c h 2 from 
Planck is higher than the WMAP-7+SPT best-fit value by about 
2.6 cr (where cr is the WMAP-7+SPT standard deviation listed 



in Table A.l) . This has an impact on the derived parameters 
Qa and Hq which differ by 3.2 cr and 2.7 cr, respectively. This 
suggests that there is a systematic difference in the power spec- 
tra/likelihoods leading to a bias in one or other analysis. The 
theoretical temperature power spectra for the two fits are shown 
in Fig. A.l. The two models are almost identical at multipoles 
I > 1000. However, the WMAP-7+SPT fit shows a substantial 
excess over the Planck fit at multipoles I % 600. The S12 spec- 
trum spans the multipole range 650 < ( < 3000, so the region 
of the WMAP-7 +SPT spectrum that is most discrepant with the 
Planck spectrum is in the multipole range that is covered exclu- 
sively by WMAP. This suggests that the problem is related in 
some way to the matching of the SPT and WMAP spectra, i.e., 
with their relative calibration, which is determined by compar- 
ing the spectra in the multipole range 600 < I < 1000 (where 
the spectra from both experiments are noisy). 

The differences between the WMAP-7+S12 and Planck base 
ACDM cosmologies lead to different types of tension with non- 
CMB data. Whereas Planck is consistent to high precision with 
the BAO data (see Fig. 15) and shows some tension with the 
Riess et al. (201 1) measurement of H , the WMAP-7+SPT best- 
fit cosmology is consistent with the Hq measurement but in se- 
rious tension with the BAO measurements. The latter point is 
illustrated by Fig. A. 2, which is equivalent to Fig. 15, but uses 
the WMAP-7 +SPT cosmology as a reference. All of the BAO 
measurements lie systematically low compared to the WMAP- 
7+SPT cosmology 44 . If the WMAP-7+SPT parameters were cor- 
rect, assuming the base ACDM cosmology, the BAO experi- 
ments must have systematically underestimated the acoustic dis- 
tance ratio by about 5%. This seems implausible for what is es- 
sentially a simple geometric measurement. H12 thus studied ex- 



43 Note that the parameter combination Cl m h 3 , which is orthogonal to 
the acoustic-scale degeneracy in flat models, is nearly identical in the 
two fits. 

44 H12 quote a 2.3% probability of compatibility between the BOSS 
measurement and the WMAP-7+SPT data 
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Table A.l. Comparison of the base ACDM parameters from Planck+WP+highL compared to the parameters determined by S12 
from a joint likelihood analysis of SPT and WMAP-1 spectra. 



Planck+WV+highL WMAP-1+SVT (S12) 



Parameter Best fit 68% limit Best fit 68% limits 



100Q b /! 2 2.207 2.207 ±0.027 2.223 2.229 ± 0.037 

Clji 2 0.1203 0.1198 ±0.0026 0.1097 0.1093 ± 0.0040 

10 9 A S 2.211 2.198 ±0.056 2.143 2.142 ±0.061 

« s 0.958 0.959 ±0.007 0.963 0.962 ±0.010 

t 0.093 0.091 ±0.014 0.083 0.083 ±0.014 

1000, 1.0415 1.0415 ±0.0006 1.0425 1.0429 ± 0.0010 

fi A 0.683 0.685 ±0.017 0.747 0.750 ± 0.020 

H 67.2 67.3 ± 1.2 72.3 72.5 ±1.9 



tensions to the standard cosmology which reconcile the WMAP- 
7+SPT power spectra with the BAO measurements. 

Since Planck and the SPT S12 spectra have a large over- 
lap range at high multipoles, where both experiments have high 
signal-to-noise, there is no need to use WMAP as an intermediary 
to establish a relative calibration. We can compare the spectra di- 
rectly via a joint likelihood analysis using the same foreground 
model that is used in the main body of this paper. Since the S12 
spectrum is measured at a frequency of 150 GHz, we first present 
results using only the Planck 143 x 143 GHz spectrum in the 
Planck likelihood. This reduces sensitivity to the details of the 
foreground modelling. Apart from small colour corrections, the 
foregrounds are identical except for differences in the Poisson 
point source amplitudes. 

Absolute calibration of the SPT spectra is determined by 
comparing with the WMAP-1 spectrum in the multipole range 
600 < { < 1000. Since the spectra from both experiments are 
noisy in this range, there is a large (roughly 3%) uncertainty in 
the absolute calibration of the S12 data. Here we use a version 
of the SPT S12 likelihood which does not include marginaliza- 
tion over calibration uncertainties and self-consistently solve for 
a map calibration factor y??J between SPT and Planck. (This dif- 
fers from the analysis of Calabrese et al. 2013 who use an SPT 
covariance matrix that includes marginalization over calibration 
errors and combine with other experiments without solving for a 
relative calibration factor.) 

The results are shown in Fig. A. 3(a) 45 . The best-fit ACDM 
cosmology is very close to the best-fit base ACDM cosmology 
of the P/anc£+ WP+highL combination used in the main body 
of the paper (shown by the blue lines in the figure). The Planck 
spectra dominate the solution and the effect of the SPT points is 
to pull the best-fit model solution slightly upwards in the multi- 
pole range 650 S> I £ 1500. This is caused by the SPT points at 
t & 1 100, which sit high relative to Planck. 

We find similar results when we combine the S12 likelihood 
with the full Planck+ WP+highL likelihood. This is illustrated 
in Fig. A. 3(b). Relative to Planck, the SPT spectra lie system- 
atically high at t S= 1100. (The Planck spectrum sits high com- 
pared to the best-fit spectrum at £ ^ 2300, but in this region 
of the spectrum, foreground and beam errors become significant 
and introduce large correlations between the data points.) The 
parameter values for the Planck+S12 fits are listed in Table A. 2. 

To relate these results to the discrepancy in Fig. A. 1 we have 
also run MCMC chains with a joint WMAP-9 (V+W band)+S12 



45 In Figs. A.3 and A. 4 we use the window functions provided by S 12 
to band-average the Planck and theory data points at high multipoles. 
In Fig. A.4, we band-average the Planck and theory data points using 
the WMAP-9 binning scheme at C < 650. 



likelihood, self-consistently solving for a relative calibration fac- 
tor, )4pj/ WMA p, between WMAP-9 and S12 and using the fore- 
ground model adopted in this paper. The parameter values from 
this analysis are listed in the second column of Table A.2. (Note 
that, compared to the WMAP-1 +S12 values listed in Table A.l, 
the WMAP-9+S12 parameters move a little closer to the Planck 
parameters.) 

Figure A. 4(a) gives the "Planck view" of the problem. 
Here, the Planck and S12 data points are exactly as plotted 
in Fig. A. 3(b). The blue curve in the upper panel shows the 
WMAP-9+S12 best-fit ACDM model. In the lower panel, the 
blue curve shows the residual between this model and the best-fit 
base ACDM Planck model. The green points show the WMAP- 
9 data points lowered in amplitude by 2.49% to account for 
the scaling difference between WMAP and Planck, discussed in 
Planck Collaboration XI (2013). The blue curve provides a good 
fit to the S12 spectrum, but fails to match the Planck spectrum 
at multipoles below about 1100. The blue curve also fails to 
match the rescaled WMAP-9 points at i ^ 700. This suggests 
that the SPT spectrum has a small multipole-dependent bias of 
about 35 pK 2 over the multipole range 650-1 100. 

The "SPT view" of the problem is illustrated in Fig. A. 4(b). 
Here, the best-fit foreground model from the WMAP-9+S12 fit 
has been subtracted from the S12 spectrum and corrected for 
the relative calibration factor Jspt/wmap ( wn i cn 1S ver Y close to 
unity). The WMAP-9 points are plotted with no scaling correc- 
tion. The blue curve now provides a good fit to the WMAP-9 
and S12 data points (as expected). But now both WMAP and 
S12 are discrepant with the Planck data. If this view were cor- 
rect, one would have to argue that the Planck spectra suffer from 
a scale-dependent bias, varying by about 2% from low multi- 
poles to I ~ 2000, that is repeatable between detectors and 
frequencies, and is common to both the HFI and LFI instru- 
ments. There is no evidence for any such effect in the Planck 
data (Planck Collaboration XI 2013; Planck Collaboration XV 
2013). 

To reinforce this point and to illustrate further the con- 
sistency between the shapes of the Planck and WMAP spec- 
tra, Fig. A. 5 shows an estimate of the combined V+W band 
WMAP power spectrum computed on the same mask used for 
the 100 x 100 Planck spectrum. Here we have used a combined 
WMAP+Planck mask for point source holes. The pink points 
show the Planck 100 x 100 GHz spectrum. The WMAP points 
have been rescaled by a multiplicative factor of 0.975 and agree 
to remarkable precision point-by-point with the Planck spec- 
trum. The rms scatter between the Planck and WMAP points 
over the multipole range 50 < t < 400 is only 22 pK 2 , i.e., after 
a multiplicative scaling the two spectra are consistent to within 
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Fig. A.3. CosmoMC fits to the joint likelihoods for Planck and SPT S 12 spectra. Left: Fits using only the 143 x 143 GHz spectrum in 
the Planck likelihood. The blue points show the SPT data after recalibration and foreground subtraction, using the best-fit solution 
from the joint likelihood analysis. The pink points show the foreground-subtracted Planck 143 x 143 GHz spectrum. The lower 
panel shows the residuals with respect to the best-fit ACDM model to the Planck+SPT combined likelihoods (shown by the red 
line in the top panel) . The blue line in the bottom panel shows the residual of the best-fit P/ancfc+ WP+highL base ACDM model. 
Right: Foreground-subtracted recalibrated SPT spectra using the best-fit parameters from the likelihood analysis of the full Planck 
likelihood combined with the SPT S12 likelihood. The pink points show the best-fit Planck ACDM spectrum from Fig. 10, and the 
red line shows the best-fit Planck+WP +highL base ACDM model from the full Planck likelihood. The residuals with respect to this 
model are plotted in the lower panel. 
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Fig. A.4. Left: Planck best-fit spectrum and SPT S 12 spectrum calibrated and foreground corrected from the combined Planck+SPT 
likelihood analysis, as plotted in Fig. A. 3(b). The lower panel shows the residuals. The blue line shows the residual of the combined 
WMAP-9+SPTS12 cosmology (which is similar to the WMAP-1 +SPT S12 model shown in Fig. A.l). The green points show the 
WMAP-9 V+W foreground corrected power spectrum lowered in amplitude by 2.49%. Right: S12 spectrum corrected for fore- 
grounds and relative calibration using the best-fit solution from the combined analysis of WMAP-9 and SPT. The blue line shows 
the best-fit WMAP-9+SPT S12 cosmology. The Planck points are as plotted in Fig. A.4. The green points show the WMAP-9 V+W 
foreground corrected spectrum as computed by the WMAP team (i.e., with no multiplicative correction relative to Planck). 



61 



Planck Collaboration: Cosmological parameters 



Table A.2. Parameter constraints in ACDM models for various likelihood combinations as described in the text. The WMAP nine- 
year polarization likelihood is used in all of these fits. We use the standard foreground model, as described in Sect. 4. The last row 
lists the SPT^ 2 value for the best-fit parameters (47 data points). 



WMAP-9+S12 Planck 143+S12 Planck+S12 
Parameter 68% limit 68% limit 68% limit 



100fi b /j 2 2.239 ±0.035 2.232 + 0.031 2.203 ± 0.026 

fi c /z 2 0.1126±0.0037 0.1170±0.0027 0.1192 ± 0.0024 

10"A S 2.167 ±0.056 2.167 ±0.054 2.177 ±0.053 

n s 0.968 ±0.009 0.971 ±0.008 0.961 ± 0.007 

t 0.083 ±0.013 0.085 ±0.013 0.085 ±0.013 

100e» 1.0426 ±0.0010 1.0422 ±0.0006 1.0417 ± 0.0006 

£1 A 0.727 ±0.020 0.704 ±0.016 0.689 ± 0.015 

H 70.7 ± 1.7 68.8 ± 1.2 67.6 ±1.1 

y 1 ^. ... 0.995 ±0.004 0.994 ± 0.003 

^spt/wmap 0.999 ±0.006 

OTspt)"™ 53.0 55.7 56.3 
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Fig.A.5. Comparison of the Planck and WMAP power spec- 
tra. The green points show the combined WMAP V+W spec- 
trum computed on the same mask used for the 100 x 100 GHz 
Planck spectrum (with a combined WMAP+Planck mask for 
point source holes) after rescaling by a multiplicative factor of 
0.976. The pink points show the Planck 100 x 100 GHz spec- 
trum computed on the same mask. The red line shows the best-fit 
Planck+ WP+highL base ACDM model. The lower panel shows 
the residuals with respect to this model. The error bars on the 
WMAP points show errors from instrumental noise alone. 

about 0.5% of the primary CMB spectrum. Similar tests are de- 
scribed in greater detail in Planck Collaboration XI (2013). 

The results presented in this appendix illustrate the high ac- 
curacy to which systematic effects need to be controlled to re- 
alize high-precision cosmology from CMB observations. Even 
small systematic errors can masquerade as evidence for new 
physics. Using Planck as a calibrator, the SPT S12 spectrum 
shows some evidence of an excess relative to Planck at multi- 
poles I ^ 1 100. We can see from the x 1 values in Table A. 2 
that the Planck+Sl2 fit is disfavoured by S12 data compared 
to the WMAP-9+SPT fit, but not strongly so. However, for 
the Planck+ WP+highL base ACDM cosmology, x\yy f° r the 



WMAP-9+S12 foreground and calibration solution is 64.8. This 
is a significant shift and reflects the difference between the 
WMAP+SPT and Planck base ACDM cosmologies. This ap- 
pendix also illustrates the advantages of having a single exper- 
iment covering both low and high multipoles. The main reason 
for the discrepancy between the S12 and Planck cosmologies 
comes from the difficulties in matching SPT to WMAP in a lim- 
ited range of multipoles where both experiments are most vul- 
nerable to small systematic errors. 

Finally, it is worth mentioning that at £ <: 1800 where the 
S12 constraints become tighter than those from Planck they are 
in excellent agreement with the best-fit base ACDM cosmology 
from Planck. 



Appendix B: Dependence of cosmological 

parameters in extended models on foreground 
modelling and likelihood choices 

A large number of likelihood comparison tests on pa- 
rameters in the base ACDM cosmology are discussed 
in Planck Collaboration XV (2013). In the main body of this pa- 
per, we report constraints on a wide variety of extended mod- 
els. In many of these models the cosmological parameters are 
strongly degenerate with each other and are therefore more 
sensitive to the detailed modelling of foregrounds, frequency 
choices, and likelihood methodology. In this Appendix we dis- 
cuss briefly how one-parameter extensions of the ACDM model 
are affected by various choices. 

B. 1 . Impact of foreground priors 

Throughout this paper we have used a particular parameteriza- 
tion of the foreground model, and marginalized over the free pa- 
rameters using relatively wide priors. As discussed in Sect. 4, 
the choice of these priors is subjective and was guided by theo- 
retical expectations and by other data, particularly results from 
high-resolution CMB experiments and the early Planck analysis 
of the CIB power spectrum (Planck Collaboration XVIII 201 1). 
As discussed in Sect. 4, for Planck the dominant foregrounds are 
the Poisson contributions from unresolved point sources and the 
clustered CIB component at 217 GHz. The other components are 
of much lower amplitude and poorly constrained by Planck data 
alone. 
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Fig. B.l. Comparison of parameter constraints from 
Planck+ WP+highL for three CIB foreground models with 
different restrictions on the CIB spectral index y CIB (assumed 
to be the same in the 143 and 217 GHz channels). The top six 
panels show cosmological parameter constraints on « s (top left) 
in the base ACDM model and on single-parameter extensions of 
the ACDM model. These are very stable to the modelling of the 
CIB. Each sub-plot is obtained from an independent analysis 
of that model with CosmoMC. The lower six panels show the 
constraints on a subset of the foreground parameters in the base 
ACDM model, some of which change significantly. 



For the thermal and kinetic SZ amplitudes we have imposed 
uniform priors of < A tSZ < 10 and < A ksz < 10. These priors 
have little impact on the parameters derived for the base ACDM 
model, or on the parameters of extended cosmologies if Planck is 
combined with ACT and SPT data at high multipoles. However, 
for extended cosmologies the priors on these minority compo- 
nents do have a small impact on the cosmological parameters. 
Table B. 1 gives results obtained from doubling the width of the 
SZ priors. The constraints on the extended parameters change by 
small amounts compared to the Planck+WP entries in Table 10, 
giving an impression of the sensitivity of Planck+WP numbers 
to minority foregrounds . 

The use of additional high-f CMB data to constrain the fore- 
ground parameters depends on having a foreground model that 
can reliably extrapolate between the scales relevant for Planck 
and the smaller scales where the high-resolution experiments 



46 The constraint on A L for Planck+WP is not given in Table B.l; the 
result is A L = 1.22+°;|f (95% CL). 



have the tightest constraints. As a simple test of the model used 
in the main body of the paper, we relax here our assumption 
that the CIB spectral index is constant with a Gaussian prior 
y CIB = 0.7 ± 0.2. Any change in CIB index between small and 
larger scales could lead to a bias in the foreground model sub- 
tracted from the Planck spectra, particularly in the 217x217 GHz 
spectrum where the CIB is the dominant foreground compo- 
nent. However, as shown in Fig. B. 1, and Table B.l, the inferred 
cosmological parameters are actually remarkably insensitive to 
the details of the model, with very similar results obtained with 
no y CIB prior and allowing a free running of the spectral index 
through the parameter dy cm /din (. (Note that we have assumed 
here that the ^-dependence of the CIB is the same, up to an am- 
plitude, at the different frequencies.) The interpretation of the 
foregrounds does change significantly, and indeed there is mild 
evidence for running of the CIB spectral index, but it has almost 
no impact on the cosmology. This should not be too surprising, 
since the CIB signal is frequency-dependent unlike the cosmo- 
logical signal, but nonetheless it is reassuring that degeneracies 
with, for example, the SZ amplitudes do not indirectly cause bi- 
ases in cosmological parameters. 

We have not investigated extensively the impact of varying 
the tSZ and kSZ templates. A variety of different approaches (an- 
alytic, semi-analytic and numerical) have been used to estimate 
tSZ templates (e.g., Komatsu & Seljak 2002; Shaw et al. 2010; 
Sehgal et al. 2010; Trac et al. 201 1; Battaglia et al. 2010, 2012). 
These have similar shapes at multipoles I ^ 3000, relevant to 
Planck and to the tSZ template used here. The shape of our tem- 
plate is also a good match to the power spectrum of the Planck 
Compton-y map over the multipole range 100 $ I ^ 1000 
(Planck Collaboration XXI 2013). The normalization of the tSZ 
templates (i.e., their dependence on <x 8 ) and their shapes at mul- 
tipoles I ^ 3000, depend on uncertain gas physics (including 
energy injection from AGN). For this reason, we have not at- 
tempted to link the amplitude of the tSZ template to the ampli- 
tude of the matter power spectrum in the parameter analyses. The 
tSZ template used here is similar in shape to the Battaglia et al. 
(2010) template that has been used extensively in the analysis of 
ACT and SPT data. The effects of varying tSZ and kSZ templates 
on high-resolution CMB experiments have been investigated by 
Dunkley et al. (201 1), Reichardt et al. (2012), and Dunkley et al. 
(2013) who find very little effect on cosmological parameters. 

B.2. WMAP low-£ polarization likelihood 

The large-scale polarization from nine years of WMAP obser- 
vations (Bennett et al. 2012) provides our most powerful con- 
straint on the reionization optical depth t. As shown in Sect. 3 
it is not essential to use WMAP polarization information to ob- 
tain tight constraints on cosmological parameters from Planck. 
Using WMAP does, however, improve the constraints on the am- 
plitude of the power spectrum and, via the partial parameter de- 
generacies sensitive to the relative amplitude of large and small- 
scale power and the amount of lensing, WMAP polarization also 
has an impact on other cosmological parameters. Most directly, 
since the small-scale CMB power scales roughly with e~ lT <Tg, 
the inferred value of <x 8 is approximately proportional to e T as 
discussed in Sect. 3.4. 

The polarization measurement at low multipoles is challeng- 
ing because of the high level of polarized Galactic foregrounds, 
so it is important to assess the impact if the assumed constraint 
were slightly wrong. Figure B.2 shows how cosmological pa- 
rameters shift if instead of using WMAP polarization we impose 
a prior t = 0.07 ± 0.013 (lc), which is about the same width 
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CamSpec Plik /,„» = 2000 / min = 1200 no 217x217 t = 0.07 ±0.013 Running CIB < A isz < 20 

Parameter 95% limits 95% limits 95% limits 95% limits 95% limits 95% limits 95% limits 95% limits 

Im,. [eV] < 0.663 < 0.691 < 0.398 < 0.600 < 0.485 < 0.768 < 0.581 < 0.999 

X* 3 - 36 3S 33632 2.89!»;g 3.30^ 2.99^ 3.23** 3.32*g 3.67*g 

* 0.2663;°*° 0.2543;°« 0233^ 0.7&«™ 0232*»« 0359*™ 0.2643™ 

dn^ldtok -0.0153°» -0.013!°;°« -0.0073™ -0.017^ -0.0053°;? -0.0123™ -0.0143°" -0.0163°» 

^ 1-233S 1-2631 1-3831 1-313^ l - 26 -o% l - 2 ^S l-M^JS 



Table B.l. Constraints on one-parameter extensions of the ACDM model from Planck with various likelihood variations. 
P/a«c£+ WP+highL is used for these entries except for the column listing result for t = 0.07 + 0.013, where the WMAP polar- 
ization likelihood is replaced by the prior, and the final column, which does not include the high-£ experiments and doubles the 
default width of the priors of the fiat priors on the two SZ amplitudes A tSZ and A ksz . The running CIB model has no prior on y cm 
and allows for spectral curvature through the parameter dy clB /dln£. 



Planck+WP P/andr+WP+highL 

- ■ Planck+ir = 0.07 ± 0.013) - ■ Planck+(r = 0.07 ± 0.013)+highL 
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Fig. B.2. Effect on cosmological parameter constraints of re- 
placing the WMAP low-£ polarization likelihood with a prior of 
t = 0.07+0.013, which prefers lower values of the optical depth. 
The top-left sub-plot is n s in the base ACDM model, while the 
others are for one-parameter extensions. Each sub-plot shows re- 
sults from independent CosmoMC analyses of the corresponding 
model. 

but about 0.02 lower than the posterior obtained from WMAP 
(as for example might be obtained if there were some residual 
foreground contamination or instrument systematic). Table. B.l 
shows the corresponding impact on parameter constraints with 
high-/" CMB data added to Planck. The shifts are consistent with 
the known parameter degeneracies, and the ACDM constraints 
on cr 8 would shift downwards by a factor of approximately e 02 , 
or about 2%. For this reason, in Sect. 3.4 we have quoted a con- 
straint on e~ T cr%, which is insensitive to possible errors in the 
large-scale polarization likelihood. A change in the Planck cal- 
ibration would also have a similar direct effect on the inferred 
physical amplitudes. 

B.3. Planck likelihood 

The results of this paper are based on the CamSpec likeli- 
hood, which includes information from 100, 143 and 217 GHz 
channels, with a range of multipole cuts as reviewed briefly 
in Sect. 2.3 and summarized in Table 6. The combination of 
channels used allows a number of foreground parameters to be 
partially determined from Planck data alone, and as discussed 



in Sect. 4 the foreground parameters can be determined more 
precisely by including additional information from other high-^ 
data sets. Planck Collaboration XV (2013) discusses and com- 
pares the CamSpec likelihood with an alternative cross- and auto- 
spectrum likelihood, Plik. The Plik likelihood uses identical 
masks over the frequency range 100-217 GHz (retaining 48% of 
the sky), ignores correlations between multipoles and uses differ- 
ent assumptions to estimate instrument noise and to correct for 
Galactic dust, as described in Planck Collaboration XV (2013). 
For the base ACDM model, these two likelihood codes give al- 
most identical results. In this section, we investigate briefly how 
cosmological parameters in extended ACDM models vary be- 
tween the two likelihood codes and for some data cuts in the 
CamSpec likelihood. 

The results are summarized in Fig. B.3 and Table B.l. Shifts 
in parameter values are expected, since different combinations of 
data are being used and hence have differing amounts of noise 
and cosmic variance. With fewer frequencies overlapping at any 
angular scale, the foreground parameters are less well deter- 
mined, and any degeneracies with foreground parameters are ex- 
pected to open up. Fig. B.3 shows that there are noticeable shifts 
in parameter values, but when also including additional high-^ 
CMB data, which constrain the foreground parameters to higher 
precision, the relatively small differences between CamSpec and 
Plik are significantly reduced. The two likelihood codes agree 
well, even for extended models when the high-£ CMB data are 
added to the likelihoods. 

However, the preference from the temperature power spec- 
trum for Al > 1 actually becomes stronger on adding high- 
£ data. Reducing £ m . dx to 2000 in the CamSpec likelihood also 
shifts Al to higher values, particularly with the addition of the 
high-^ data. We do not, at this stage, have a full understanding 
of these shifts in A L . We note that the high-^ ACT data itself 
favours high Al (Sievers et al. 2013), and this may be part of the 
reason behind the shift to high values when high-^ data are in- 
cluded. (The truncation of Planck spectra at Z max = 2000 also 
limits the accuracy of matching high-f data to Planck. There is 
then no overlapping multipole range with the SPT data and a 
significantly narrower overlap range for ACT.) 

The analysis of extended ACDM models with strong param- 
eter degeneracies is complex and sensitive to small systematic 
errors in the CMB data and to errors in the foreground model. In 
the analysis of extended models we have usually quoted results 
from the Planck+ WP+highL data combination, using the full 
CamSpec likelihood. With this combination, we utilize the high 
signal-to-noise ratio of the Planck spectra at 143 and 217 GHz, 
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Fig. B.3. Upper, comparison of the ACDM constraints on n s (top-left) and single-parameter extensions of the ACDM model for a 
variety of data cuts for Planck+WP. Each sub-plot is obtained from a separate CosmoMC analysis of the corresponding model. The 
dashed lines show the results from Plik, an alternative likelihood discussed in Planck Collaboration XV (2013), run here with the 
same SZ and CIB foreground priors as for the CamSpec results. For the extended models, the value of the additional parameter in 
the base ACDM model is shown with the vertical dashed lines. Lower, same as the upper set of panels, but for Planck+WP+highL. 
Additional data from the high-^ CMB experiments significantly reduce the foreground degeneracies. 
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which have a wide multipole range with which to match high- 
£ experiments to Planck, and gain better control of foregrounds 
via the inclusion of high-^ data. The general "rule of thumb" 
adopted in this paper has been to use the differences between pa- 
rameter constraints from Planck+WP and Planck+ WP+highL 
as a guide to whether parameters are sensitive to errors in the 
foreground model (or other sources of error). We do sometimes 
see shifts of up to around 1 <x between these likelihoods in the 
parameter values for extended models and this needs to be borne 
in mind when interpreting our results. In the absence of any ad- 
ditional information, we take the cosmological parameters from 
Planck+ WP+highL as our best estimates for extended models. 
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